code 31055;
procedure GRNNEW(N,X,A);
value N; integer N; array X,A;
begin
procedure ELMCEV(L,U,SHIFT,A,B,X);
value L,U,SHIFT,X; integer L,U,SHIFT;
real X; array A,B;
for L:=L step -1 until U do A[L]:=A[L]+B[L+SHIFT]*X;
integer K;
for K:=N-1 step -1 until 0 do
ELMCEV(N-1,N-1-K,1,A,A,X[N-1-K]);
end GRNNEW
eop