rsolve (y(n+2)-y(n+1)-y(n),y(n))

2**(-n)*(C0*(1 + sqrt(5))**n + C1*(-sqrt(5) + 1)**n)

Want to compute something more complicated? Try a full Python/SymPy console at SymPy Live.