p10_13.m

来自「Programs for the book Advanced Engineeri」· M 代码 · 共 13 行

% P10_13.M Plot the rabbit population
clear
k=[0:12];    % 12 months
y= 1/5^0.5*((1+5^0.5)/2).^k - 1/5^0.5*((1-5^0.5)/2).^k;
fprintf('\n     Months    Population\n')
for i=2:1:13
fprintf('        %g      %g\n',i-1,y(i))
end 
%
plot(k,y)
title('Growth of rabbit population')
xlabel('Months')
ylabel('Population')