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