semilogy(fibonacci2(18), '-o')
% We can either get the slope as the average (i.e., mean) of all of the
% pairwise differences (i.e., diff) in the sequence of logs of Fibonacci
% numbers. Note that differences of y-values are sufficient since the
% difference in the x-values is always 1 here.
slope = mean(diff(log10(fibonacci2(18))));
disp(sprintf('Approximate slope of line from average is %f',slope))
% or just as the last difference (which seems to be very stable)
slopes = diff(log10(fibonacci2(18)));
disp(sprintf('Approximate slope of line from last value is %f',slopes(end)))

disp(sprintf('For comparison purposes, the base 10 log of the golden ratio is: %f',log10((1+sqrt(5))/2)))
