xn1axn.m

来自「分形计算的几个程序原码」· M 代码 · 共 18 行

hold on; 
aspan=[3.4,3.99]; 
set(gcf,'DoubleBuffer','on'); 
axis([aspan,-1,2]); 
x0=0.3; 
for a=linspace(aspan(1),aspan(2),200); 
   x=x0; 
   for k=1:20; 
       x=a*x*(1-x); 
   end 
   title(['a=',num2str(a)]); 
   for k=1:100; 
       x=a*x*(1-x); 
       plot(a,x,'k.','markersize',2); 
       pause(0.02); 
   end 
end