Составить процедуру и функцию для вычисления значений элементов матрицы по заданной формуле

uses crt;
const n=10;
      m=10;
      xn=0;
      xk=1;
      yn=0;
      yk=1;
var a:array[1..n,1..m] of real;
    i,j:byte;
    x,y,dx,dy:real;
begin
clrscr;
dx:=(xk-xn)/n;
dy:=(yk-yn)/m;
writeln('Њ*ваЁж* A:');
for i:=1 to n do
 begin
  for j:=1 to m do
   begin
    x:=xn+(i-1)*dx;
    y:=yn+(j-1)*dy;
    a[i,j]:=abs(exp(x/(1+abs(x*y))));
    write(a[i,j]:6:3);
   end;
  writeln;
 end;
readln
end.

uses crt;
const n=10;
      m=10;
      xn=0;
      xk=1;
      yn=0;
      yk=1;
function Z(x,y:real):real;
begin
Z:=exp(x)/(1+abs(x*y));
end;
type mtr=array[1..n,1..m] of real;
procedure MatricaA(var a:mtr);
var i,j:byte;
    x,y,dx,dy:real;
begin
dx:=(xk-xn)/n;
dy:=(yk-yn)/m;
writeln('Матрица A:');
for i:=1 to n do
 begin
  for j:=1 to m do
   begin
    x:=xn+(i-1)*dx;
    y:=yn+(j-1)*dy;
    a[i,j]:=abs(Z(x,y));
    write(a[i,j]:6:3);
   end;
  writeln;
 end;
end;
var a:mtr;
    i,j:byte;
begin
clrscr;
MatricaA(a);
readln
end.

uses crt;
const n=10;
      m=10;
      xn=0;
      xk=1;
      yn=0;
      yk=1;
function Z(x,y:real):real;
begin
Z:=exp(x)/(1+abs(x*y));
end;
type mtr=array[1..n,1..m] of real;
procedure MatricaA(var a:mtr);
var i,j:byte;
    x,y,dx,dy:real;
begin
dx:=(xk-xn)/n;
dy:=(yk-yn)/m;
writeln('Матрица A:');
for i:=1 to n do
 begin
  for j:=1 to m do
   begin
    x:=xn+(i-1)*dx;
    y:=yn+(j-1)*dy;
    a[i,j]:=Z(x,y);
    write(a[i,j]:6:3);
   end;
  writeln;
 end;
end;
procedure MatricaB(a:mtr;var b:mtr);
var x:array[1..m] of real;
    i,j,k:byte;
    z:real;
begin
writeln('Матрица B:');
for i:=1 to n do
 begin
  for j:=1 to m do
  x[j]:=a[i,j];
  for j:=1 to m-1 do
  for k:=j+1 to m do
  if x[j]>x[k] then
   begin
    z:=x[j];
    x[j]:=x[k];
    x[k]:=z;
   end;
  for j:=1 to m do
   begin
    b[i,j]:=x[j];
    write(b[i,j]:6:3);
   end;
  writeln;
 end;
end;
 
var a,b:mtr;
begin
clrscr;
MatricaA(a);
MatricaB(a,b);
readln
end.

uses crt;
const n=10;
      m=10;
      xn=0;
      xk=1;
      yn=0;
      yk=1;
function Z(x,y:real):real;
begin
Z:=exp(x)/(1+abs(x*y));
end;
type mtr=array[1..n,1..m] of real;
     vc1=array[1..1,1..m] of real;
     vc2=array[1..m,1..1] of real;
procedure MatricaA(var a:mtr);
var i,j:byte;
    x,y,dx,dy:real;
begin
dx:=(xk-xn)/n;
dy:=(yk-yn)/m;
writeln('Матрица A:');
for i:=1 to n do
 begin
  for j:=1 to m do
   begin
    x:=xn+(i-1)*dx;
    y:=yn+(j-1)*dy;
    a[i,j]:=Z(x,y);
    write(a[i,j]:6:3);
   end;
  writeln;
 end;
end;
procedure MatricaB(a:mtr;var b:mtr);
var x:array[1..m] of real;
    i,j,k:byte;
    z:real;
begin
writeln('Матрица B:');
for i:=1 to n do
 begin
  for j:=1 to m do
  x[j]:=a[i,j];
  for j:=1 to m-1 do
  for k:=j+1 to m do
  if x[j]>x[k] then
   begin
    z:=x[j];
    x[j]:=x[k];
    x[k]:=z;
   end;
  for j:=1 to m do
   begin
    b[i,j]:=x[j];
    write(b[i,j]:6:3);
   end;
  writeln;
 end;
end;
procedure VectorX(b:mtr;var x:vc1);
var i:byte;
begin
writeln('Вектор-строка X:');
for i:=1 to m do
 begin
  x[1,i]:=(b[1,1]+b[1,m])/2;
  write(x[1,i]:6:2);
 end;
writeln;
end;
procedure VectorX1(x:vc1;var x1:vc2);
var i:byte;
begin
writeln('Вектор-столбец X1:');
for i:=1 to m do
 begin
  x1[i,1]:=x[1,i];
  writeln(x1[1,i]:6:2);
 end;
writeln;
end;
var a,b:mtr;
    x,xa:vc1;
    x1:vc2;
    res:real;
    i,j:byte;
begin
clrscr;
MatricaA(a);
MatricaB(a,b);
write('Press Enter');
readln;
VectorX(b,x);
VectorX1(x,x1);
{Вычисление М, ходя так ее назвать нельзя, поскольку m уже занято}
{умножаем x*A}
for i:=1 to m do
 begin
  xa[1,i]:=0;
  for j:=1 to n do
  xa[1,i]:=xa[1,i]+x[1,i]*a[j,i];
 end;
{умножаем a*x1}
res:=0;
for i:=1 to m do
res:=res+a[1,i]*x1[i,1];
write('M=',res:0:3);
readln
end.