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