Equalisation

function out=equagamma(img)
 maxcount=sum(img)
 maxcntniv=length(img)/256.0
 inlevel=0    // continue input  level
 outlevel=0   // discrete output level
 routlevel=0  // continue output level
 [rows,cols] = size(img)
 out=zeros(rows,cols)

 while inlevel<256
    for index=find(img==inlevel)
        out(index)=outlevel
        cntniv=cntniv+1
        while cntniv>=maxcntniv
            cntniv=cntniv-maxcntniv
            routlevel=routlevel+1
        end

    end
    outlevel=routlevel
    inlevel=inlevel+1
 end
endfunction