cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A118239 Engel expansion of cosh(1).

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