matlab code for water-filling using binary searching method

N=1024;
NoB=1/N;
g2n=chi2rnd（1[1 N]）;
g2n=sort（g2n‘descend‘）;
delta=10^-5;
% 后面要用到的一些求和中间变量
sum=0;
for i=1:N    
    b（i）=NoB/g2n（i）;
    sum=sum+NoB/g2n（i）;
end
N_bar=N;
for Pmax=0:0.1:2
   c=0;
    d=max（g2n）/NoB;
    u_bar=（c+d）/2;
    for n=1:N
        Pn（n）=max（[1/u_bar-NoB/g2n（n） 0]） ;
    end
    sum_pn=0;
    for n=1:N
        sum_pn=sum_pn+Pn（n）;
    end