Умножение матрицы на вектор

uses crt;
const n=16;
var a:array[1..n]of real;
i,k,j:integer;
x:real;
begin
clrscr;
randomize;
writeln('ish mass:');
for i:=1 to n do
begin
a[i]:=random(10)+random-4;
write(a[i]:6:2);
end;
writeln;
k:=0;
for i:=1 to n do
if a[i]>=0 then {если элемент положительный}
begin
k:=k+1; {увеличиваем счетчик 1,2,3...}
x:=a[i]; {запоминаем элемент}
for j:=i downto k+1 do {сдвигаем все элементы от найденоого до к=1,2,3.. вправо}
begin
a[j]:=a[j-1]; 
end;
a[k]:=x; {на 1,2,3.. место вставляем наш элемент}
end;
writeln('rez:');
for i:=1 to n do
write(a[i]:6:2);
readln
end.

uses crt;
var A:array[1..50,1..50] of real;
b,Ab:array[1..50]of real;
n,i,j:integer;
begin
clrscr;
randomize;
write('Vvedite razmernost n=');
readln(n);
writeln('Matrica A:');
for i:=1 to n do
begin
for j:=1 to n do
begin
A[i,j]:=random(10)+random-3;
write(A[i,j]:6:2);
end;
writeln;
end;
writeln('Vector b:');
for i:=1 to n do
begin
if odd(i) then b[i]:=1/(i*i+2) {если нечетное}
else if not odd(i) then b[i]:=1/i; {если четное}
write(b[i]:6:2);
end;
writeln;
writeln('Vector Ab:');
for i:=1 to n do
begin
Ab[i]:=0;
for j:=1 to n do
Ab[i]:=Ab[i]+A[i,j]*b[j]; {произведение матрицы на вектор}
write(Ab[i]:6:2);
end;
readln
end.