code 31022;
    procedure MULCOL(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]:= B[L,J]*X;
        eop