不好意思!请教中

yanyun2929

for N1=1:1:100
    s(N1)=0;
end
r1=0.43;r2=0.5;r3=0.5;
a1=0.01;a2=0.02;a3=0.04;
N2=1;N3=4;
for  N1=1:100
    for T=1:2000
        n1(N1)=binornd(N1,r1);
        n2=binornd(N2,r2);
        n3=binornd(N3,r3);
        if r1*a1*N1<a1*n1
            B1=r1*a1*N1;
        else B1=0;end
        if r2*a2*N2<a2*n2
            B2=r2*a2*N2;
        else B2=0;end
        if r3*a3*N3<a3*n3
            B3= r3*a3*N3;
        else B3=0;end
        B=B1+B2+B3;
        if B~=0
            d1=(r1*a1*N1)/B;我想使这边的N1也从１变到１００，怎么编？谢谢
            d2=(r2*a2*N2)/B;
            d3=(r2*a3*N3)/B;
            wlr1=max(r1*a1*N1,a1*n1)-a1*n1;
            wlr2=max(r2*a2*N2,a2*n2)-a2*n2;
            wlr3=max(r3*a3*N3,a3*n3)-a3*n3;
            wlr=wlr1+wlr2+wlr3+1-0.01-(r1*a1*N1+r2*a2*N2+r3*a3*N3);
            t1=(wlr*d1+r1*a1*N1)/a1;
            t2=(wlr*d2+r2*a2*N2)/a2;
            t3=(wlr*d3+r3*a3*N3)/a3;
            if n1(N1)>0
                p1=t1/n1(N1);
            else p1=1;
            end
            if n2>0
                p2=t2/n2;
            else p2=1;
            end
            if n3>0
                p3=t3/n3;
            else p3=1;
            end
            nt1=binornd(n1(N1),p1);
            nt2=binornd(n2,p2);
            nt3=binornd(n3,p3);
        else nt1=n1(N1);
             nt2=n2;
             nt3=n3;
          if nt1*a1+nt2*a2+nt3*a3<1-0.01;
              s(N1)=nt1+nt2+nt3+s(N1);
          end
      end
      s(N1)=s(N1)/2000;
end
plot(s)
end