procedure vsprod(var s: scalar; var v1,v2: vector);
var
  n: integer;
  st: scalar;
begin { vsprod }
  for n:=0 to vdim do sprod(s,v1[n],v2[n]);
  sabs0(v2[vdim],st);
  v2[vdim]:=st;
end { vsprod };