procedure vvmtprod(var m: matrix; var vv1, vv2: vtuple);
var
  i,j: integer;
  v0,v1,v2: vector;
begin { vvmtprod }
  for i:=0 to lmax do begin
    vzero(v0);
    for j:=0 to lmax do begin
      vsprod(m[j,i],vv1[j],v1);
      vsum(v0,v1,v2);
      v0:=v2;
    end;
    vv2[i]:=v0;
  end;
end { vvmtprod };

