function fib,n
  q=1/( p=(1+sqrt(5))/2 )
  return,round((p^n+q^n)/sqrt(5))
end