%函数functiondxdt=differential(t,x)dxdt=[x(1)*(1-x(1))*(5-6*x(2));x(2)*(1-x(2))*(1-6*x(1))];end%主程序clear%y-xfori=0:0.1:1forj=0:0.1:1[T,Y]=ode45('differential',[05],[ij]);figure(1)gridonplot(Y(:,1),Y(:,2));holdonendend%dx/dt-tfori=0:0.1:1forj=0:0.1:1[T,Y]=ode45('differential',[05],[ij]);figure(2)gridonplot(T,Y(:,1));holdonendend%dy/dt-tfori=0:0.1:1forj=0:0.1:1[T,Y]=ode45('differential',[05],[ij]);figure(3)gridonplot(T,Y(:,2));holdonendendY-X图像dx/dt-t图像dy/dt-t图像