code 34183;
    comment MCA 2413;
    procedure REASCL(A, N, N1, N2); value N, N1, N2;
    integer N, N1, N2; array A;
    begin integer I, J; real S;
        for J:= N1 step 1 until N2 do 
        begin S:= 0;
            for I:= 1 step 1 until N do 
            if ABS(A[I,J]) > ABS(S) then S:= A[I,J];
            if S ^= 0 then 
            for I:= 1 step 1 until N do A[I,J]:= A[I,J] / S
        end 
    end REASCL;
        eop