rsolve (y(n+2)-y(n+1)-y(n),y(n),{y(0):0,y(1):1})

2**(-n)*sqrt(5)*((1 + sqrt(5))**n - (-sqrt(5) + 1)**n)/5

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