OFFSET
0,2
COMMENTS
LINKS
FORMULA
a(2*n) = A007696(n+1)^2 = ( Product {k = 0..n} 4*k + 1 )^2.
a(2*n-1) = 4*n*A007696(n)^2 = 4*n * ( Product {k = 0..n-1} 4*k + 1 )^2.
a(n) = 4*a(n-1) + (2*n - 1)^2*a(n-2) with a(0) = 1, a(1) = 4.
a(2*n+1) = 4*(n + 1)*a(2*n); a(2*n) = (4*n + 2)*a(2*n-1) + a(2*n-2).
Empirical e.g.f.: ((-Q(1/2, -3)-Q(-1/2, -3))*P(1/2, (2*x+3)/(2*x-1))+Q(1/2, (2*x+3)/(2*x-1))*(P(1/2, -3)+P(-1/2, -3)))/((1-2*x)^(3/2)*(-Q(-1/2, -3)*P(1/2, -3)+Q(1/2, -3)*P(-1/2, -3))) where P and Q are Legendre functions of the first and second kinds. - Robert Israel, Feb 24 2015
EXAMPLE
54/25 = 2.16, 441/200 = 2.205 etc approach 2.188..
MAPLE
a[0] := 1: a[1] := 4:
for n from 2 to 18 do a[n] := 4*a[n-1] + (2*n-1)^2*a[n-2] end do:
seq(a[n], n = 0 .. 18);
MATHEMATICA
RecurrenceTable[{a[0] == 1, a[1] == 4, a[n] == 4 a[n - 1] + (2 n - 1)^2 a[n - 2]}, a, {n, 20}] (* Vincenzo Librandi, Feb 24 2015 *)
PROG
(Magma) I:=[1, 4]; [n le 2 select I[n] else 4*Self(n-1)+(2*n-3)^2*Self(n-2): n in [1..20]]; // Vincenzo Librandi, Feb 24 2015
CROSSREFS
KEYWORD
nonn,frac,easy
AUTHOR
Peter Bala, Feb 23 2015
STATUS
approved