code 34212;
    procedure DAVUPD(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:= V[K] * C1; WK:= W[K] * C2;
            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 DAVUPD;
        eop