OFFSET
0,2
LINKS
Gheorghe Coserea, Table of n, a(n) for n = 0..201
Robert S. Maier, On Rationally Parametrized Modular Equations, arXiv:math/0611041 [math.NT], 2006.
FORMULA
n^2 * a(n) = 3*(39*n^2-52*n+20) * a(n-1) - 441*(3*n-4)^2 * a(n-2), with a(0)=1, a(1)=21.
0 = 9*x*(x^2+13*x+49)*y'' + (21*x^2 + 195*x + 441)*y' + (4*x+21)*y, where y(x) = A(x/-441).
EXAMPLE
A(x) = 1 + 21*x + 693*x^2 + 23940*x^3 + ... is the g.f.
MATHEMATICA
a[0] = 1; a[1] = 21; a[n_] := a[n] = (3(39n^2 - 52n + 20) a[n-1] - 441(3n - 4)^2 a[n-2])/n^2;
Table[a[n], {n, 0, 16}] (* Jean-François Alcover, Oct 19 2018 *)
PROG
(PARI)
seq(N) = {
my(a = vector(N)); a[1] = 21; a[2] = 693;
for (n=3, N,
a[n] = (3*(39*n^2 - 52*n + 20) * a[n-1] - 441*(3*n-4)^2 * a[n-2])/n^2);
concat(1, a);
};
seq(17)
(Magma) I:=[21, 693]; [1] cat [n le 2 select I[n] else (3*(39*n^2-52*n+20)*Self(n-1)-441*(3*n-4)^2*Self(n-2)) div n^2: n in [1..30]]; // Vincenzo Librandi, Aug 25 2016
CROSSREFS
KEYWORD
sign
AUTHOR
Gheorghe Coserea, Aug 23 2016
STATUS
approved