请教一下，多谢！！急！！

yanyun2929

这个程序的结果出不来，请各位高手帮帮忙！！谢谢！！r1=0.43;r2=0.5;r3=0.5;
a1=0.01;a2=0.02;a3=0.04;
for N1=1:100;
    N2=4;N3=1;
    s(N1)=0;
   for T=1:2000;
   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=(a1*r1*N1)/B;
    d2=(a2*r2*N2)/B;
    d3=(a3*r3*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-(a1*r1*N1+a2*r2*N2+a3*r3*N3);
    t1=(wlr*d1+a1*r1*N1)/a1;
    t2=(wlr*d2+a2*r2*N2)/a2;
    t3=(wlr*d3+a3*r3*N3)/a3;
    p1=t1/n1;p2=t2/n2;p3=t3/n3;
    nt1=binornd(n1,p1);
    nt2=binornd(n2,p2);
    nt3=binornd(n3,p3);
  else   t1=r1*N1;t2=r2*N2;t3=r3*N3;
    p1=t1/n1;p2=t2/n2;p3=t3/n3;
    nt1=binornd(n1,p1);
    nt2=binornd(n2,p2);
    nt3=binornd(n3,p3);
    end
    if nt1*a1+nt2*a2+nt3*a3<1-0.01;
        s(N1)=nt1+nt2+nt3+s(N1);
        end
        s(N1)=s(N1)/2000;
    end
end   
plot(s)