code 34213;
procedure FLEUPD(H, N, V, W, C1, C2); value N, C1, C2;
integer N; real C1, C2; array H, V, W;
begin integer I, J, K;
real VK, WK;
K:= 0; for J:= 1, J + K while K < N do
begin K:= K +1; VK:= - W[K] * C1 + V[K] * C2; WK:= V[K] * C1;
for I:= 0 step 1 until K - 1 do
H[I + J]:= H[I + J] + V[I + 1] * VK -W[I + 1] * WK
end
end FLEUPD;
eop