rsolve(y(n+2)-y(n+1)-y(n), y(n), {y(0): 0, y(1): 1})
Want to compute something more complicated? Try a full Python/SymPy console at SymPy Live.