code 34021;
    procedure ELMVECCOL(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]:= A[L] + B[L,I] * X;
        eop