OFFSET
0,2
COMMENTS
a(n) = (a(n-1) + a(n+1))/(2n). - Emeric Deutsch, Feb 03 2006
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
FORMULA
a(-n) = (-1)^n * a(n). - Michael Somos, Jan 25 2014
a(n+1) ~ 1.0903447... * 2^n * n! * f(1/(4*n)) where f(x) = 1 + x + x^2/2! + 13*x^3/3! + 241*x^4/4! + 6201*x^5/5! + ...). - Michael Somos, Jan 26 2014
0 = a(n) * (a(n+2)) + a(n+1) * (-a(n+1) + 2*a(n+2) - a(n+3)) + a(n+2) * (a(n+2)) for all n in Z. - Michael Somos, Jan 25 2014
0 = u0 * u3 - u1 * (u2 + 2*u3) + u2 * (4*u2) where u0 = e.g.f. a(n), u1=u0', u2=u1', u3=u2'. - Michael Somos, Jan 25 2014
EXAMPLE
a(3) = 2*2*a(2)-a(1) = 12-2 = 10.
G.f. = 1 + 2*x + 3*x^2 + 10*x^3 + 57*x^4 + 446*x^5 + 4403*x^6 + 52390*x^7 + ...
MAPLE
a[0]:=1: a[1]:=2: for n from 1 to 20 do a[n+1]:=2*n*a[n]-a[n-1] od: seq(a[n], n=0..20); # Emeric Deutsch, Feb 03 2006
MATHEMATICA
RecurrenceTable[{a[0]==1, a[1]==2, a[n]==2(n-1)*a[n-1]-a[n-2]}, a, {n, 20}] (* Harvey P. Dale, May 31 2016 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Amarnath Murthy, Apr 19 2004
EXTENSIONS
More terms from Emeric Deutsch, Feb 03 2006
STATUS
approved