
clear		
xlab = '$P_n$'; 		% X-label
ylab = '$P_{n+1}$'; 			% Y-label
mytitle = 'Skeena Sockeye Salmon - Ricker''s Model'; 	% Title
yr = [1908 1912 1916 1920 1924 1928 1932 1936 1940 1944 1948];
salmon = [1098 740 714 615 706 510 278 448 528 639 523];
n = length(yr);
salmon1 = salmon(1:n-1);
salmon2 = salmon(2:n);
pn = linspace(0,1200,100);
pn1 = 1.5344*pn.*exp(-0.0007816*pn);
pn2 = 1.32765*pn - 0.0006146*pn.^2;

plot(salmon1,salmon2,'ro','MarkerSize',7);
hold on
plot([0 800],[0 800],'r--','LineWidth',1.5);
plot(pn,pn2,'b-','LineWidth',1.5);
plot(pn,pn1,'-','color',[0 0.5 0],'LineWidth',1.5);

grid
legend('Data','Identity','Logistic','Ricker''s','Location','Southeast');
%text(5,80,'h(a) = 6.45a + 73.9','FontSize',fs_text);

axis equal
xlim([0 1200]);
ylim([0 800]);
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,'FontSize',12);

print -depsc salmonupx.eps
