Даны веществ квадратная матрица A порядка n, n-мерный веществ вектор x

uses crt;
const n=5;
Type Matrix=array [1..n,1..n] of real;
Vector=array [1..n] of real;
var a: Matrix;
x,ax,b: Vector;
i,j: integer;
max: real;
procedure MultAx(a: Matrix; x: Vector; var ax: Vector);
begin
for i:=1 to n do
begin
ax[i]:=0;
for j:=1 to n do
ax[i]:=ax[i]+a[i,j]*x[j];
end;
end;
function Modul(x: Vector): real;
var s: real;
begin
s:=0;
for i:=1 to n do
s:=s+sqr(zx[i]);
Modul:=sqrt(s);
end;
begin
writeln('Matrix:');
for i:=1 to n do
begin
for j:=1 to n do 
begin
a[i,j]:=random(9)+1;
write(a[i,j]:6:2);
end;
writeln;
end;
writeln('Vector:');
for i:=1 to n do
begin
x[i]:=random(9)+1;
write(x[i]);
end;
writeln;
for i:=1 to n do
for j:=1 to n do
MultAx(a[i,j],x[i],ax[i]);
writeln('Vector Ax:');
for i:=1 to n do
writeln(ax[i]);
readln;
b[i]:=ax[i]-x[i];
writeln('Vector Ax-x');
writeln(b[i]);
readln;
if Modul(ax[i])>Modul(x[i]) then
max:=Modul(ax[i]);
else  max:=Modul(x[i]);
if Modul(b[i])>Modul(ax[i]) then
max:=Modul(b[i]);
else max:=Modul(ax[i]);
if Modul(b[i])>Modul(x[i]) then
max:=Modul(b[i]);
else max:=Modul(x[i]);
writeln('max:=', max);
readln;
end.