login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A122057 a(n) = (n+1)! * (H(n+1) - H(2)), where H(n) are the harmonic numbers. 1

%I #31 Sep 08 2022 08:45:28

%S 0,2,14,94,684,5508,49104,482256,5185440,60668640,767940480,

%T 10462227840,152698210560,2377651449600,39350097561600,

%U 689874448435200,12773427499929600,249097496204390400,5103595024496640000,109608397522606080000

%N a(n) = (n+1)! * (H(n+1) - H(2)), where H(n) are the harmonic numbers.

%C Former title (corrected): A Legendre-based recurrence sequence. Let b(n) = ((4*n+2)*x -(2*n+1) )/(n+1)*b(n-1) - (n/(n+1))*b(n-2), where x=1, then a(n) = (n+1)!*b(n)/6. - _G. C. Greubel_, Oct 03 2019

%D Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964, 9th Printing (1970), pp. 782

%H G. C. Greubel, <a href="/A122057/b122057.txt">Table of n, a(n) for n = 1..445</a>

%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].

%F Let b(n) = ((-2*n-1) +(4*n+2)*x)/(n+1)*b(n-1) - (n/(n+1))*b(n-2) with x=1, then a(n) = b(n)*(n+1)!/6.

%F a(n) = (n+1)! * Sum_{k=3..n+1} 1/k. - _Gary Detlefs_, Jul 15 2010

%p a:=n-> (n+1)!*add(1/k,k=3..n+1): seq(a(n),n=1..30); # _Gary Detlefs_, Jul 15 2010

%t x=1; b[1]:=0; b[2]:=2; b[n_]:= b[n]= ((-2*n-1) +(4*n+2)*x)/(n+1)*b[n-1] - (n/(n+1))*b[n-2]; Table[b[n]*(n+1)!/6, {n,30}]

%t Table[(n+1)!*(HarmonicNumber[n+1] - 3/2), {n,30}] (* _G. C. Greubel_, Oct 03 2019 *)

%o (PARI) vector(30, n, (n+1)!*(sum(k=1,n+1, 1/k) - 3/2) ) \\ _G. C. Greubel_, Oct 03 2019

%o (Magma) [Factorial(n+1)*(HarmonicNumber(n+1) - 3/2): n in [1..30]]; // _G. C. Greubel_, Oct 03 2019

%o (Sage) [factorial(n+1)*(harmonic_number(n+1) - 3/2) for n in (1..30)] # _G. C. Greubel_, Oct 03 2019

%o (GAP) List([1..30], n-> Factorial(n+1)*(Sum([1..n+1], k-> 1/k) - 3/2) ); # _G. C. Greubel_, Oct 03 2019

%Y Cf. A001008, A002805.

%K nonn

%O 1,2

%A _Roger L. Bagula_, Sep 14 2006

%E If all terms are really negative, sequence should probably be negated. - _N. J. A. Sloane_, Oct 01 2006

%E Negated terms and edited by _G. C. Greubel_, Oct 03 2019

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 07:27 EDT 2024. Contains 371265 sequences. (Running on oeis4.)