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') 各位大神帮帮忙,,,告诉我错在哪里应该怎么解决