% Simulate LV Model
% Requires Stiff Solver ode15s
[t1,y]=ode15s(@lv_ode,[0,10],[1000,600]);

figure(101)
plot(t1,y);grid;
xlim([0,10]);
fontlabs = 'Times New Roman';   
xlabel('$t$','FontSize',14,'FontName',fontlabs,...
    'interpreter','latex');  
ylabel('Molecules','FontSize',14,'FontName',fontlabs);  
set(gca,'FontSize',12);
print -depsc lv_ode_ts.eps

figure(102)
plot(y(:,1),y(:,2),'b-');grid;
xlabel('$Y_1$','FontSize',14,'FontName',fontlabs,...
    'interpreter','latex');  
ylabel('$Y_2$','FontSize',14,'FontName',fontlabs,...
    'interpreter','latex');
print -depsc lv_ode_pp.eps