OFFSET
0,4
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..54
R. W. Gosper and Richard C. Schroeppel, Somos Sequence Near-Addition Formulas and Modular Theta Functions, arXiv:math/0703470 [math.NT], 2007. See p. 18.
MAPLE
f:=proc(n) option remember;
if n <= 2 then n-1
elif n=3 then 12
elif n=4 then -432
elif n=5 then 93312
else (144*f(n-3)*f(n-1)+432*f(n-2)^2)/f(n-4); fi; end;
(144*a(n-3)*a(n-1)+432*a(n-2)^2)/a(n-4); ai; end;
[seq(f(n), n=0..15)];
MATHEMATICA
a[n_] := a[n] = Which[n <= 2, n - 1, n == 3, 12, n == 4, -432, n == 5, 93312, True, (144*a[n - 3]*a[n - 1] + 432 a[n - 2]^2)/a[n - 4]];
Table[a[n], {n, 0, 12}] (* Jean-François Alcover, Dec 05 2017 *)
nxt[{a_, b_, c_, d_}]:={b, c, d, (144b*d+432c^2)/a}; Join[{-1, 0}, NestList[nxt, {1, 12, -432, 93312}, 10][[All, 1]]] (* Harvey P. Dale, Dec 30 2017 *)
PROG
(Magma) I:=[1, 12, -432, 93312]; [-1, 0] cat [n le 4 select I[n] else (144*Self(n-3)*Self(n-1) + 432*Self(n - 2)^2)/Self(n-4): n in [1..20]]; // G. C. Greubel, Aug 08 2018
CROSSREFS
KEYWORD
sign
AUTHOR
N. J. A. Sloane, May 16 2014
STATUS
approved