Матрицы

[Родик]

Помогите с Delphi

Соседями элемента Aij в матрице N на N назовем элементы Akl с i-1<=k<=i+1, j-1<=l<=j+1, (k,l)<>(i,j). Операция сглаживания матрицы дает новую матрицу того же размера, каждый элемент которой получается как среднее арифметическое имеющихся соседей соответствующего элемента исходной матрицы.

Матрица генерируется с помощью датчика псевдослучайных чисел.
Матрица квадратная, поэтому вводим только один размер.

Вот начало кода где формируется матрица:

unit Unit1;

interface

uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls, Math;

type
TForm1 = class(TForm)
Label1: TLabel;
Edit1: TEdit;
Button1: TButton;
Memo1: TMemo;
procedure Button1Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;

var
Form1: TForm1;

implementation

{$R *.dfm}

procedure TForm1.Button1Click(Sender: TObject);
var
i,j,n,k,l:integer;
aij:array of array of double;
akl:array of array of double;
s:string;

begin
n:=strtoint(edit1.Text);
memo1.Clear;
setlength(aij,n);
setlength(akl,n);

for i:=0 to n-1 do
setlength(aij,n);
for i:=0 to n-1 do
setlength(akl,n);

//////////////////////////////////////////////////////////
for i:=0 to n-1 do
for j:=0 to n-1 do
begin
aij[i,j]:=-15+30*random;
aij[i,j]:=RoundTo(aij[i,j],-2);
end;
//////////////////////////////////////////////////////////
memo1.Lines.Add('matrix');
for i:=0 to n-1 do
begin
s:='';
for j:=0 to n-1 do
s:=s+floattostr(aij[i,j])+' ';
memo1.Lines.Add(s);
end;

//////////////////////////////////////////////////////////
for i:=0 to n-1 do
for j:=0 to n-1 do

aij:=nil;
akl:=nil;

end;
initialization
randomize;

end.

[ex.cluz]

А где сам вопрос?