. mata : function fib(n) { return((((1+sqrt(5))/2):^n-((1-sqrt(5))/2):^n)/sqrt(5)) } : fib(0..10) 1 2 3 4 5 6 7 8 9 10 11 +--------------------------------------------------------+ 1 | 0 1 1 2 3 5 8 13 21 34 55 | +--------------------------------------------------------+ : end