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