Определите u=Ax+By-Cz+Bx

uses crt;
const nmax=10;
type mtr=array[1..nmax,1..nmax] of real;
     vkt=array[1..nmax] of real;
var n:byte;
procedure vv_mtr(var a:mtr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
a[i,j]:=9*random
end;
procedure vv_vek(var x:vkt);
var i:byte;
begin
for i:=1 to n do
x[i]:=9*random
end;
procedure um_mt_vk(a:mtr;x:vkt;var ax:vkt;cm,cv:string);
var i,j:byte;
begin
writeln('Матрица ',cm);
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j]:5:2);
  writeln;
 end;
gotoXY(n*5+2,1+n div 2);
write('x');
gotoXY(n*5+3,1);
writeln('Вектор ',cv);
for i:=1 to n do
 begin
  gotoXY(n*5+4,1+i);
  write(x[i]:5:2);
 end;
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;
gotoXY(n*5+11,1+n div 2);
write('=');
gotoXY(n*5+14,1);
writeln(cm,'*',cv);
for i:=1 to n do
 begin
  gotoXY(n*5+13,1+i);
  write(ax[i]:5:2);
 end;
gotoXY(1,n+2);
writeln('Нажмите Enter');
readln;
clrscr
end;
var a,b,c:mtr;
    x,y,z,ax,by,cz,bx,u:vkt;
    i:byte;
begin
clrscr;
randomize;
repeat
write('Размер матриц и векторов до ',nmax,' n=');
readln(n);
until n in [1..nmax];
vv_mtr(a);
vv_mtr(b);
vv_mtr(c);
vv_vek(x);
vv_vek(y);
vv_vek(z);
clrscr;
um_mt_vk(a,x,ax,'A','x');
um_mt_vk(b,y,by,'B','y');
um_mt_vk(c,z,cz,'C','z');
um_mt_vk(b,x,bx,'B','x');
writeln('  u  ');
for i:=1 to n do
 begin
  u[i]:=ax[i]+by[i]-cz[i]+bx[i];
  writeln(u[i]:5:2);
 end;
readln
end.