code 34022;
    procedure ELMCOLVEC(L, U, I, A, B, X); value L,U,I,X;
    integer L,U,I; real X; array A,B;
    for L:= L step 1 until U do A[L,I]:= A[L,I] + B[L] * X;
        eop