procedure vdquot(var r: scalar; var v1,v2: vector);
var
  i: integer;
begin { vdquot }
  write('(vdquot'); flush(output);
  for i:=0 to vdim-2 do v2[i]:=v1[i+1];
  ssymm(v1[vdim],v2[vdim-1]);
  squot(v1[vdim],r,v2[vdim]);
  isdiv(vdim+1,v2[vdim]);
  write(')'); flush(output);
end { vdquot };