典型的虫口混沌问题。状态的递推关系为Xn+1=a*(Xn-Xn^2) 

Xn+1和Xna分别为次年和当年的虫口数。 x1=0.5 , n=50。用图形的方法研究系数a的变化为a=[2.0 3.2 3.5 3.8]对xn+1变化的影响。图形横、纵坐标分别为n和Xn+1。 

我的程序如下: 

n=50; 
a=[2.0 3.2 3.5 3.8]; 
x(1)=0.5; 
[p,k]=size(a); 
for j=1:k 
for i=1:n 
x(i+1)=a(j)*(x(i)-x(i)^2) 
plot(i,x) 
end 
end 

此程序只能求出n=50时的情况,望高手指点。