code 31045;
procedure ALLORTPOL(N,X,B,C)RESULTS:(P);
value N,X; integer N;real X;array B,C,P;
if N=0 then P[0]:=1else 
begin integer K,K1; real R,S,H;
  R:=P[1]:=X-B[0]; S:=P[0]:=1;K:=1;
  for K1:=2 step 1 until N do 
begin H:=R;  P[K1]:=R:=(X-B[K])*R-C[K]*S;
  S:=H; K:=K1;
end;
end ALLORTPOL

        eop