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.

Showing 1-3 of 3 results.

A174503 Continued fraction expansion for exp( Sum_{n>=1} 1/(n*A087799(n)) ), where A087799(n) = (5+sqrt(24))^n + (5-sqrt(24))^n.

Original entry on oeis.org

1, 8, 1, 96, 1, 968, 1, 9600, 1, 95048, 1, 940896, 1, 9313928, 1, 92198400, 1, 912670088, 1, 9034502496, 1, 89432354888, 1, 885289046400, 1, 8763458109128, 1, 86749292044896, 1, 858729462339848, 1, 8500545331353600, 1
Offset: 0

Views

Author

Paul D. Hanna, Mar 21 2010

Keywords

Examples

			Let L = Sum_{n>=1} 1/(n*A087799(n)) or, more explicitly,
L = 1/10 + 1/(2*98) + 1/(3*970) + 1/(4*9602) + 1/(5*95050) +...
so that L = 0.1054740177896236251618898675297390156061405857647...
then exp(L) = 1.1112372317482311056432125938345153306039099019639...
equals the continued fraction given by this sequence:
exp(L) = [1;8,1,96,1,968,1,9600,1,95048,1,940896,1,...]; i.e.,
exp(L) = 1 + 1/(8 + 1/(1 + 1/(96 + 1/(1 + 1/(968 + 1/(1 +...)))))).
Compare these partial quotients to A087799(n), n=1,2,3,...:
[10,98,970,9602,95050,940898,9313930,92198402,912670090,9034502498,...].
		

Crossrefs

Programs

  • PARI
    {a(n)=local(L=sum(m=1,2*n+1000,1./(m*round((5+sqrt(24))^m+(5-sqrt(24))^m))));contfrac(exp(L))[n]}

Formula

a(2n-2) = 1, a(2n-1) = A087799(n) - 2, for n>=1 [conjecture].
The above conjectures are correct. See the Bala link for details. - Peter Bala, Jan 08 2013
a(n) = 11*a(n-2)-11*a(n-4)+a(n-6). G.f.: -(x^4+8*x^3-10*x^2+8*x+1) / ((x-1)*(x+1)*(x^4-10*x^2+1)). [Colin Barker, Jan 20 2013]

A174508 Continued fraction expansion for exp( Sum_{n>=1} 1/(n*A086594(n)) ), where A086594(n) = (4+sqrt(17))^n + (4-sqrt(17))^n.

Original entry on oeis.org

1, 7, 65, 1, 535, 4353, 1, 35367, 287297, 1, 2333751, 18957313, 1, 153992263, 1250895425, 1, 10161155671, 82540140801, 1, 670482282087, 5446398397505, 1, 44241669462135, 359379754094593, 1, 2919279702218887, 23713617371845697, 1
Offset: 0

Views

Author

Paul D. Hanna, Mar 21 2010

Keywords

Examples

			Let L = Sum_{n>=1} 1/(n*A086594(n)) or, more explicitly,
L = 1/8 + 1/(2*66) + 1/(3*536) + 1/(4*4354) + 1/(5*35368) +...
so that L = 0.1332613701545977545822925541573311424901819508933...
then exp(L) = 1.1425485874089841897117810754210805471767735522069...
equals the continued fraction given by this sequence:
exp(L) = [1;7,65,1,535,4353,1,35367,287297,1,2333751,...]; i.e.,
exp(L) = 1 + 1/(7 + 1/(65 + 1/(1 + 1/(535 + 1/(4353 + 1/(1 +...)))))).
Compare these partial quotients to A086594(n), n=1,2,3,...:
[8,66,536,4354,35368,287298,2333752,18957314,153992264,...].
		

Crossrefs

Programs

  • PARI
    {a(n)=local(L=sum(m=1,2*n+1000,1./(m*round((4+sqrt(17))^m+(4-sqrt(17))^m))));contfrac(exp(L))[n]}

Formula

a(3n-3) = 1, a(3n-2) = A086594(2n-1) - 1, a(3n-1) = A086594(2n) - 1, for n>=1 [conjecture].
a(n) = 67*a(n-3)-67*a(n-6)+a(n-9). G.f.: -(x^2-x+1)*(x^6-8*x^5-8*x^4-2*x^3+72*x^2+8*x+1) / ((x-1)*(x^2+x+1)*(x^6-66*x^3+1)). [Colin Barker, Jan 20 2013]
From Peter Bala, Jan 25 2013: (Start)
The above conjectures are correct. The real number exp( Sum {n>=1} 1/(n*A086594(n)) ) is equal to the infinite product F(x) := product {n >= 0} (1 + x^(4*n+3))/(1 - x^(4*n+1)) evaluated at x = sqrt(17) - 4. Ramanujan has given a continued fraction expansion for the product F(x). Using this we can find the simple continued fraction expansion of the numbers F(1/2*(sqrt(N^2 + 4) - N)), N a positive integer. The present case is when N = 8. See the Bala link for details.
The theory also provides the simple continued fraction expansion of the numbers F({sqrt(17) - 4}^(2*k+1)), k = 1, 2, 3, ...: if [1; c(1), c(2), 1, c(3), c(4), 1, ...] denotes the present sequence then the simple continued fraction expansion of F({sqrt(17) - 4}^(2*k+1)) is given by [1; c(2*k+1), c(2*(2*k+1)), 1, c(3*(2*k+1)), c(4*(2*k+1)), 1, ...].
(End)

A174510 Continued fraction expansion for exp( Sum_{n>=1} 1/(n*A080040(n)) ), where A080040(n) = (1+sqrt(3))^n + (1-sqrt(3))^n.

Original entry on oeis.org

1, 1, 3, 1, 9, 13, 1, 37, 51, 1, 141, 193, 1, 529, 723, 1, 1977, 2701, 1, 7381, 10083, 1, 27549, 37633, 1, 102817, 140451, 1, 383721, 524173, 1, 1432069, 1956243, 1, 5344557, 7300801, 1, 19946161, 27246963, 1, 74440089, 101687053, 1, 277814197
Offset: 0

Views

Author

Paul D. Hanna, Mar 21 2010

Keywords

Examples

			Let L = Sum_{n>=1} 1/(n*A080040(n)) or, more explicitly,
L = 1/2 + 1/(2*8) + 1/(3*20) + 1/(4*56) + 1/(5*152) + 1/(6*416) +...
so that L = 0.5855329921665857283309456463364081071245363598803...
then exp(L) = 1.7959479567807442397990076546690432122217738278933...
equals the continued fraction given by this sequence:
exp(L) = [1;1,3,1,9,13,1,37,51,1,141,193,1,529,723,1,...]; i.e.,
exp(L) = 1 + 1/(1 + 1/(3 + 1/(1 + 1/(9 + 1/(13 + 1/(1 +...)))))).
Compare these partial quotients to A080040(n)/2^[n/2], n=1,2,3,...:
[2,4,10,14,38,52,142,194,530,724,1978,2702,7382,10084,27550,...],
where A080040 begins:
[2,8,20,56,152,416,1136,3104,8480,23168,63296,172928,472448,...].
		

Crossrefs

Programs

  • PARI
    {a(n)=local(L=sum(m=1,2*n+1000,1./(m*round((1+sqrt(3))^m+(1-sqrt(3))^m))));contfrac(exp(L))[n]}

Formula

a(3n-3) = 1, a(3n-2) = A080040(2n-1)/2^(n-1) - 1, a(3n-1) = A080040(2n)/2^n - 1, for n>=1 [conjecture].
a(n) = 5*a(n-3)-5*a(n-6)+a(n-9). G.f.: -(x^2-x+1)*(x^6-2*x^5-2*x^4-2*x^3+4*x^2+2*x+1) / ((x-1)*(x^2+x+1)*(x^6-4*x^3+1)). [Colin Barker, Jan 20 2013]
Showing 1-3 of 3 results.