%函数 
function dxdt=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-x 
for i=0:0.1:1     for j=0:0.1:1 
        [T,Y]=ode45('differential',[0 5],[i j]);          figure(1)         grid on 
        plot(Y(:,1),Y(:,2));          hold on     end end