code 34040;
    procedure ROTCOL(L, U, I, J, A, C, S); value L,U,I,J,C,S;
    integer L,U,I,J; real C,S; array A;
    begin real X, Y;
        for L:= L step 1 until U do 
        begin X:= A[L,I]; Y:= A[L,J]; A[L,I]:= X * C + Y * S;
            A[L,J]:= Y * C - X * S
        end 
    end ROTCOL;
        eop