%
% Calculate and plot Finbonacci numbers
%

clear 
clf

n = 50;
alpha = log((1 + sqrt(5))/2);

tic()
fff = myfibonacci(n);
toc()

x = 1:n;

semilogy(x, fff, '.-')

z = exp(alpha*x)/sqrt(5);

hold on

title('Fibonacci numbers')
xlabel('n')
ylabel('F(n)')
grid on
semilogy(x, z, 'k--')
legend('numerics', 'analytics')

hold off