be(1)=10;ne(1)=10;bf(1)=20;nf(1)=20;

N=20;Qf=0.5063; r(1:20)=0.1;

for k=1:20;

    rand('state',4800);

    A=rand(1,power(2,(k+1)));

    Q(k)=mean(A(:));

    for t=1:300;

        for i=1:20;

        p(1)=0;

        p(i+1)=p(i)+0.2

        sbe(t)=power(Qf,p(i))*power(ne(t),r(1:20)*p(i))/(power(20,power(ne(t),p(i)*r(1:20)))+19*power(20,Qf*power(Qf,p(i))*power(nf(t),p(i)*r(1:20))));

        sbf(t)=power(Qf,p(i))*power(nf(t),r(1:20)*p(i))/(power(20,power(ne(t),p(i)*r(1:20)))+19*power(20,Qf*power(Qf,p(i))*power(nf(t),p(i)*r(1:20))));

        end

        dbe(t)=5*power((1+0.05),t-1)*sbe(t)-0.001*be(t)

        dbf(t)=5*power((1+0.05),t-1)*sbf(t)-0.001*bf(t)

        dne(t)=0.1001*be(t)/1/(1-power(0.1,t))-0.001*ne(t);

        dnf(t)=0.1001*bf(t)/1.2/(1-power(0.1,t))-0.001*nf(t);

        be(t+1)=dbe(t)+ be(t);

        bf(t+1)=dbf(t)+ bf(t);

        ne(t+1)=dne(t)+ne(t);

        nf(t+1)=dnf(t)+nf(t);

end

      I(k)=k/N;

      be300(k)=be(301);

      bf300(k)=bf(301);

      ne300(k)=ne(301);

      nf300(k)=nf(301);

    x(k)=be300(k)/(be300(k)+19*bf300(k));

    y(k)=ne300(k)/(ne300(k)+19*nf300(k));

end

I(21)=[]

plot(I,x,'-.',I,y, '-o')      各位大神帮帮忙，，，告诉我错在哪里应该怎么解决