figure;

N = length(y(:,1));
init = floor(5/6*N);

% Copy over the steady state
y_ss = zeros(N-init+1,3);
y_ss(:,1) = y([init:N],1);
y_ss(:,2) = y([init:N],2);
y_ss(:,3) = y([init:N],3);

t_ss = t([init:N]);

plot3(y_ss(:,1),y_ss(:,2),y_ss(:,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;

