%I #33 Mar 07 2020 08:50:49
%S 1,2,12,30,56,90,132,182,240,306,380,462,552,650,756,870,992,1122,
%T 1260,1406,1560,1722,1892,2070,2256,2450,2652,2862,3080,3306,3540,
%U 3782,4032,4290,4556,4830,5112,5402,5700,6006,6320,6642,6972,7310,7656,8010,8372
%N Engel expansion of cosh(1).
%C Differs from A002939 only in first term.
%C This sequence is also the Pierce expansion of cos(1). - _G. C. Greubel_, Nov 14 2016
%H G. C. Greubel, <a href="/A118239/b118239.txt">Table of n, a(n) for n = 1..1000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/EngelExpansion.html">Engel Expansion</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PierceExpansion.html">Pierce Expansion</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,1).
%F a(n) = A002939(n-1) = 2*(n-1)*(2*n-3) for n>1.
%F From _Colin Barker_, Apr 13 2012: (Start)
%F a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
%F G.f.: x*(1 - x + 9*x^2 - x^3)/(1-x)^3. (End)
%F E.g.f.: -6 + x + 2*(3 - 3*x + 2*x^2)*exp(x). - _G. C. Greubel_, Oct 27 2016
%t Join[{1}, Table[(2 n - 2) (2 n - 3), {n, 2, 50}]] (* _Bruno Berselli_, Aug 04 2015 *)
%t Join[{1}, LinearRecurrence[{3,-3,1},{2,12,30},25]] (* _G. C. Greubel_, Oct 27 2016 *)
%t PierceExp[A_, n_] := Join[Array[1 &, Floor[A]], First@Transpose@ NestList[{Floor[1/Expand[1 - #[[1]] #[[2]]]], Expand[1 - #[[1]] #[[2]]]} &, {Floor[1/(A - Floor[A])], A - Floor[A]}, n - 1]]; PierceExp[N[Cos[1] , 7!], 50] (* _G. C. Greubel_, Nov 14 2016 *)
%o (PARI) a(n)=max(4*n^2-10*n+6, 1) \\ _Charles R Greathouse IV_, Oct 22 2014
%o (Sage)
%o A118239 = lambda n: falling_factorial(n*2,2) if n>0 else 1
%o print([A118239(n) for n in (0..46)]) # _Peter Luschny_, Aug 04 2015
%Y Cf. A006784, A002939, A068377.
%K nonn,easy
%O 1,2
%A _Eric W. Weisstein_, Apr 17 2006