%函数
function dxdt=differential2(t,x,flag,X,Y,Z,T,M,N,O,P)
dxdt=[x(1)*(1-x(1))*[X+Y-Z-T];x(2)*(1-x(2))*[M+N-O-P]];

%y-x
E(1)=2000;r(1)=0.03;e(1)=600;
d(1)=0.08;d(2)=0.02;t(1)=0.8;k(1)=0.06;c(1)=0.4;
m=0.05;n=0.02;
E(2)=800;r(2)=0.01;e(2)=350;t(2)=0.7;k(2)=0.04;c2=0.2;
p=25;
X=x(2)*(E(1)+r(1)*e(1)*p+d(1)*e(1)+d(1)*e(2)-0.5*c(1)*e(1)*e(1)-(1-t(1))*k(1)*r(1)*e(1)*p)
Y=(1-x(2))*(E(1)+r(1)*e(1)*p-0.5*c(1)*e(1)*e(1)-m*E(1))
Z=x(2)*(E(1)+d(1)*e(2)-n*E(1))
T=(1-x(2))*(E(1)-m*E(1)-n*E(1))
M=x(1)*(E(2)+r(2)*e(2)*p+d(2)*e(2)+d(2)*e(1)-0.5*c(2)*e(2)*e(2)-(1-t(2))*k(2)*r(2)*e(2)*p)
N=(1-x(1))*(E(2)+r(2)*e(2)*p-0.5*c(2)*e(2)*e(2)-m*E(2))
O=x(1)*(E(2)+d(2)*e(1)-n*E(2))


for i=0:0.1:1
    for j=0:0.1:1

P=(1-x(1))*(E(2)-m*E(2)-n*E(2))
        [T,Y]=ode45('differential2',[0 5],[i j]);
        figure
        grid on
        plot(Y(:,1),Y(:,2));
        hold on
    end
end