code 34357;
    procedure ROTCOMCOL(L, U, I, J, AR, AI, CR, CI, S);
    value L, U, I, J, CR, CI, S; integer L, U, I, J;
    real CR, CI, S; array AR, AI;
    begin real ARLI, AILI, ARLJ, AILJ;
        for L:= L step 1 until U do 
        begin ARLI:= AR[L,I]; AILI:= AI[L,I]; ARLJ:= AR[L,J];
            AILJ:= AI[L,J];
            AR[L,I]:= CR * ARLI + CI * AILI - S * ARLJ;
            AI[L,I]:= CR * AILI - CI * ARLI - S * AILJ;
            AR[L,J]:= CR * ARLJ - CI * AILJ + S * ARLI;
            AI[L,J]:= CR * AILJ + CI * ARLJ + S * AILI;
        end 
    end ROTCOMCOL;
        eop