三方进化博弈数值分析

请问以下程序有问题，怎么改呀？？？求大神指点！！

%函数

function dxdt=differential(t,x)

dxdt=[x(1)*(1-x(1))*(2*x(3)+1;x(2)*(1-x(2))*(3*x(3)+1;x(3)*(1-x(3))*(-1*x(1)-x(2)-x(1)*x(2)+4)];

end

%主程序

clear

%z-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

plot3(Y(:,1),Y(:,2),Y(:,3));

hold on

end

end

%dx/dt-t

for i=0:0.1:1

for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]);

figure(2)

grid on

plot3(T,Y(:,1));

hold on

end

end

%dy/dt-t

for i=0:0.1:1

for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]);

figure(3)

grid on

plot3(T,Y(:,2));

hold on

end

end

%dz/dt-t

for i=0:0.1:1

for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]);

figure(4)

grid on

plot3(T,Y(:,3));

hold on

end

end