clear
clf

n = 50;
F = myfibonacci(n);

alpha = log((1 + sqrt(5))/2);
np = 100;
x = linspace(1, n, np);
y = exp(alpha*x)/sqrt(5);

semilogy(1:n, F, '.-')
hold on
semilogy(x, y)
grid on
title('Fibonacci numbers')
xlabel('n')
ylabel('F_n')
legend('recurrence', 'analytic', 'location', 'best')