function f = fibonacci(limit)
%% Fibonacci sequence
f = [1 1]; pos = 1;
while f(pos) + f(pos+1) < limit
   f(pos+2) = f(pos) + f(pos+1);
   pos = pos + 1;
end