
xlab = 'Population ($P$)';             % X-label
ylab = '$f(P)$';                     % Y-label
mytitle = 'Phase Portrait';   % Title

xx = linspace(-500,2500,100);

ms=7;
yy = 0.05*xx.*(1-xx/2000);

plot(xx,yy,'b-','LineWidth',2.5);
hold on
plot([0 0],[-10 30],'k-','LineWidth',1.5);
plot([-500 2500],[0 0],'k-','LineWidth',1.5);
plot([0 2000],[0 0],'ro','MarkerSize',10);
plot([-500 2500],[0 0],'r-','LineWidth',3);
plot([0],[0],'wo','MarkerSize',8);
plot([2000],[0],'ro','MarkerFaceColor','r','MarkerSize',9);

fs_big=30;
dx=0.05;
text(250,dx,'$>$','FontSize',fs_big,'Color','r','interpreter','latex');
text(700,dx,'$>$','FontSize',fs_big,'Color','r','interpreter','latex');
text(1150,dx,'$>$','FontSize',fs_big,'Color','r','interpreter','latex');
text(1600,dx,'$>$','FontSize',fs_big,'Color','r','interpreter','latex');
text(2220,dx,'$<$','FontSize',fs_big,'Color','r','interpreter','latex');
text(-280,dx,'$<$','FontSize',fs_big,'Color','r','interpreter','latex');
text(750,5,'$\frac{dP}{dt} = f(P)$','FontSize',20,'interpreter','latex');
text(15,-2,'$0$','FontSize',20,'interpreter','latex');
text(1850,-2,'$M$','FontSize',20,'interpreter','latex');

xlim([-500 2500]);
ylim([-10 30]);

fontlabs = 'Times New Roman';   % Font type used in labels
xlabel(xlab,'FontSize',14,'FontName',fontlabs,'interpreter','latex');  
ylabel(ylab,'FontSize',14,'FontName',fontlabs,'interpreter','latex'); 
title(mytitle,'FontSize',16,'FontName',fontlabs,'interpreter','latex'); 

set(gca,'ytick',[]);
set(gca,'xtick',[]);

print -depsc popphase.eps   

