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