
plot3(y(:,1),y(:,2),y(:,3));
hold on;

xlabel('N_1');ylabel('N_2');zlabel('P');
title('Phase space in N_1, N_2, P. (a=1.5; b=0.01; c=0.5; d=1.0; r=1.0; eps=0.009)');
grid on;

N = length(y(:,1));
% Starting position marked with green;
% Ending position marked with red;
plot3(y(1,1),y(1,2),y(1,3),'g*');
plot3(y(N,1),y(N,2),y(N,3),'r*');

