function rho_M = rhoMAX(n,h,m)

AL_b = 16;
rho = zeros(m);

for l=1:m
    
    AJ_b = 0;
    
    
    for j=1:m
        %Define the Al & Aj
        Al = bitget(uint8(AL_b), n:-1:1);
        Aj = bitget(uint8(AJ_b), n:-1:1);
        % Define the Rho value for each and every Al & Aj Compination 
        rho(l,j) = Rho(Al,Aj,n,h);
       
        AJ_b = AJ_b + 1;
    end

    AL_b = AL_b + 1;
end

rho_M = max(max(rho));