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

[Илья Кибель]

Дана процедура умножения матрицы на вектор:

procedure MultMatVecString(a:aType; b,c:bType);
begin
for i:=1 to m do begin
c[i]:=0;
for j:=1 to m do begin
c[i]:=c[i]+a[i,j]*b[j];
end;
end;
end;

Вызов процедуры:

MultMatVecString(n,u,a);

При вызове процедуры на выходе получается нулевой вектор. Помогите, пожалуйста, найти ошибку.

[Maslan]

А что там под m прячется?

[Илья Кибель]

m=4 - константа, размерность квадратной матрицы.

[Virtson]

procedure MultMatVecString(a:aType; b:bType; Var c: bType);

// OFF : с фантаизией туговато ..

[Илья Кибель]

Я попробовал. Ничего не получилось. Вектор на выходе по-прежнему нулевой.

[Maslan]

Всё работает.... Может ты не инициализировал матрицы с первым вектором?

program Project2;
{$APPTYPE CONSOLE}
uses
  SysUtils;
Type
 TMatrix=array of array of real;
 TVector=array of real;
procedure MultMatVecString(a:tMatrix; b:tVector; var c:TVector);
Var i,j:integer;
begin
SetLEngth(c,Length(b));
for i:=0 to Length(b)-1 do begin
c[i]:=0;
for j:=0 to Length(b)-1 do begin
c[i]:=c[i]+a[i,j]*b[j];
end;
end;
end;

Var Matr:TMatrix;Vec1,Vec2:TVector;
Var i,j:integer;
begin
Randomize;
 SetLength(Matr,4);
 SetLength(Matr[0],4);
 SetLength(Matr[1],4);
 SetLength(Matr[2],4);
 SetLength(Matr[3],4);
 SetLength(Vec1,4);
 for i:=0 to 3 do for j:=0 to 3 do
 Matr[i,j]:=Random;
 for i:=0 to 3 do Vec1[i]:=Random;
 MultMatVecString(Matr,Vec1,Vec2);
 WriteLn('[',Vec2[0]:3:4,'][',Vec2[1]:3:4,'][',Vec2[2]:3:4,'][',Vec2[3]:3:4,']');
 Readln;
end.