procedure pexp(d2: integer; var s: scalar; var v: vector);
var
  i: integer;
begin { pexp }
  v[0]:=sone;
  v[1]:=s;
  for i:=2 to d2 do begin
    sprod(s,v[i-1],v[i]);
    isdiv(i,v[i]);
  end;
  for i:=d2+1 to vdim do szero(v[i]);
end { pexp };