procedure lmprod(var m: matrix; var l1,l2: ltuple);
var
  i,j: integer;
  s0,s1,s2: scalar;
begin { lmprod }
  for i:=0 to lmax do begin
    szero(s0);
    for j:=0 to lmax do begin
      sprod(m[i,j],l1[j],s1);
      ssum(s0,s1,s2);
      s0:=s2;
    end;
    l2[i]:=s0;
  end;
end { lmprod };