
clear		
xlab = 'Year'; 		
ylab = 'Population ($\times 1000$)'; 			
mytitle = 'Skeena Sockeye Salmon - Simulation'; 
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);
popl = [1096.8];
popr = [1103.7];
for i = 2:n
    popl(i) = 1.32765*popl(i-1) - 0.0006146*popl(i-1).^2;
    popr(i) = 1.5344*popr(i-1).*exp(-0.0007816*popr(i-1));
end

plot(yr,salmon,'ro','MarkerSize',7);
hold on
plot(yr,popl,'bo-','LineWidth',1.5);
plot(yr,popr,'o-','color',[0 0.5 0],'LineWidth',1.5);

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

xlim([1908 1948]);
ylim([0 1200]);
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 salmonx.eps
