function f = fib(n) f = [1 1 ; 1 0]^(n-1); f = f(1,1); end