fib(n)=my(phi=quadgen(5));(phi^n-(-1/phi)^n)/(2*phi-1)