%
% calculate and plot Fibonacci numbers
%
clear
clf

format compact

n = 50;

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

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

x = 1:n;

semilogy(x, fff, '.-')
title('Fibonacci numbers')
xlabel('n')
ylabel('F(n)')
legend('numerics')
grid on
hold on
semilogy(x, exp(x*alpha), 'k--');
hold off