1 | initial version |
In your example all f(n) are equal to 1, so the limit is 1 :)
2 | No.2 Revision |
In your example all f(n) are equal to 1, so the limit is 1 :)
In the case of linear recursive definitions you can use rsolve from sympy or solve_rec from Maxima, and then use limit function. In the general nonlinear case you can experiment numerically: http://ask.sagemath.org/question/1247/how-to-define-a-recursive-sequence-on-sage.
Using Maxima you can do:
sage: maxima('f[n]:=1/5*(f[n-1]^2+f[n-1]+3)')
f[n]:=(3+f[n-1]+f[n-1]^2)/5
sage: maxima('f[0]:1')
1
sage: maxima('f[100]')
1