code 31047;
real procedure SUMORTPOL(N,X,B,C,A);
value N,X; integer N; real X; array B,C,A;
if N=0 then SUMORTPOL:=A[0] else 
begin integer K; real H,R,S;
  R:=A[N]; S:=0;
  for K:=N-1 step -1 until 1 do 
  begin 
    H:=R;
    R:=A[K]+(X-B[K])*R+S;
    S:=       -C[K]*H
  end;
  SUMORTPOL:=A[0]+(X-B[0])*R+S
end SUMORTPOL;
        eop