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

