procedure vscale(var s: scalar; var v1,v2: vector);
var
  n: integer;
  st: scalar;
  v0: vector;
begin { vscale }
  write('(vscale'); flush(output);
  vpowers(s,v0);
  for n:=0 to vdim do sprod(v0[n],v1[n],v2[n]);
  st:=v2[vdim];
  sabs0(st,v2[vdim]);
  write(')'); flush(output);
end { vscale };

