Получить квадратную матрицу порядка 10

uses crt;
const n=10;
var a:array[1..n,1..n] of real;
    x,p:real;
    i,j:byte;
begin
write('x=');
read(x);
p:=1;
for i:=1 to n do
 begin
  a[1,i]:=p;//заполняем периметры
  a[i,1]:=p;
  a[n,n-i+1]:=p;
  a[n-i+1,n]:=p;
  p:=p*x; //повышаем степень х
 end;
writeln('Матрица:');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
     // заполним середину
    if not( i in [1,n])and not(j in [1,n]) then a[i,j]:=0;
    write(a[i,j]:8:2);
   end;
  writeln;
 end;
end.

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, Grids, StdCtrls, Spin;
 
type
  TForm1 = class(TForm)
    Label1: TLabel;
    Button1: TButton;
    StringGrid1: TStringGrid;
    Edit1: TEdit;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
const n=10;
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
var a:array[1..n,1..n] of real;
    x,p:real;
    i,j:integer;
begin
x:=Strtofloat(Edit1.Text);
StringGrid1.ColCount:=n;
StringGrid1.RowCount:=n;
p:=1;
for i:=1 to n do
 begin
  a[1,i]:=p;//заполняем периметры
  a[i,1]:=p;
  a[n,n-i+1]:=p;
  a[n-i+1,n]:=p;
  p:=p*x; //повышаем степень х
 end;
for i:=1 to n do
for j:=1 to n do
 begin
  // заполним середину
  if not( i in [1,n])and not(j in [1,n]) then StringGrid1.Cells[j-1,i-1]:='0'
  else StringGrid1.Cells[j-1,i-1]:=FloatToStrF(a[i,j],fffixed,8,2);
 end;
end;
 
end.