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