options=odeset(‘RelTol’,1.0e-3);
[t,x]=ode23(@f,[0 20],[2 1],options);
plot(t,x(:,1),’k–‘,t,x(:,2),’r-‘)
xlabel(‘t’,’fontsize’,16);
ylabel(‘x’,’fontsize’,16);
legend(‘x_0(t)’,’x_1(t)’);
axis([0 20 0 6]);
function dxdt=f(t,x)
dxdt=[1.2*x(1)-0.6*x(1)*x(2);
-0.8*x(2)+0.3*x(1)*x(2)];
end