matlab 虫口方程特性

function test
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 研究不同R和初始值x对周期的影响，横坐标是R的变化值，纵坐标是X（t）的变化值
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
R=input（‘R=  ‘）;    %给定R的初始值，观察不同的R，对周期的变化
x1=input（‘x1=  ‘）;    %虫口方程对初始值很敏感，给定不同的初始值
x2=input（‘x2=  ‘）;    %虫口方程对初始值很敏感，给定不同的初始值
%x=0.001;
hold on;
a=zeros（1100）;
b=zeros（1100）;
   i=1;
   while x1<=1;
    if i==101   % 判断x是否等于1，是跳出while循环，否就需执行
        break;
    else
        dx1=R.*（1-x1）.*x1;
        x1=dx1;