function f = fibonacci(n)
f = zeros(n,1);
f(1) = 0;
f(2) = 1;
for k = 3:n
f(k) = f(k-1) + f(k-2);
end
