Logistic模型参数的matlab求解

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）*