function?f=diedai（NNmN0n）
r=0;
for?i=1:n
????r=r+i*log（（Nm/N（i）-1）/（Nm/N0-1））;
end
r=-6/n/（n+1）/（2*n+1）*r;
a=Nm;b=r;
e=1;
while（e>0.000001）
????y=0;z=0;
????for?i=1:n
????????z=z+1/（（N0*exp（r*i）+a-N0）/（N0*exp（r*i）-N（i）））^2;
????????y=y+（exp（r*i）-1）*