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