code 34020;
    procedure ELMVEC(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;
        eop