A332396 Decimal expansion of lim_{n->infinity} (b(n)/n - (n - log(n))/4) where {b(n)} is the real-valued sequence defined by b(1)=1 and b(n+1) = b(n) + sqrt(b(n)) for n > 0.

6, 7, 5, 1, 7, 7, 4, 4, 2, 4, 5, 8, 5, 5, 7, 1, 3, 9, 8, 1, 3, 2, 8, 5, 6, 2, 5, 0, 7, 5, 8, 2, 7, 6, 3, 3, 6, 8, 4, 0, 7, 3, 1, 5, 8, 9, 8, 9, 0, 5, 1, 5, 9, 1, 4, 0, 1, 1, 3, 0, 8, 1, 0, 8, 5, 9, 1, 5, 1, 9, 7, 2, 9, 4, 9, 8, 4, 0, 2, 3, 7, 9, 3, 5, 3, 3, 7

Offset

0,1

Comments

Consider the "add the square root" sequence {b(n)} starting with 1, i.e., b(1)=1, b(n+1) = b(n) + sqrt(b(n)) for n > 0: at each step, the next term is obtained simply by adding the current term to its square root, so the sequence begins {1, 2, 2 + sqrt(2), 2 + sqrt(2) + sqrt(2 + sqrt(2)), ...}, i.e., {1, 2, 3.414213..., 5.261972..., ...}.

The 1st, 10th, 100th, and 1000th terms are

b(10^0) = 1.0

b(10^1) = 25.956...

b(10^2) = 2452.850...

b(10^3) = 248949.869...

b(10^4) = 24983729.385...

b(10^5) = 2499779700.562...

Since the difference between successive terms is b(n+1) - b(n) = sqrt(b(n)) and the similar differential equation (d/dx)f(x) = sqrt(f(x)) is satisified by the function f(x) = (1/4)*x^2, it is perhaps not surprising that lim_{n->infinity} b(n)/n^2 = 1/4. More precisely, it can be shown that, as n increases, b(n) approaches

(1/4)*n^2

+ u*n

+ u^2 - u/2

+ (-(1/2)*u^2 + u/4 - 1/96)/n

+ ((1/3)*u^3 + 0*u^2 + (-5/48)*u + 7/576)/n^2

+ ...

where u = -(1/4)*log(n) + c

and c = 0.675177442458557139813285625075...

It follows that lim_{n->infinity} (b(n)/n - (n - log(n))/4) = c.

If we were to define {b(n)} instead as b(1)=1, b(n+1) = b(n) + 1/b(n) for n > 0 (i.e., "add the reciprocal" rather than "add the square root"), we would obtain the rational-valued sequence {b(n)} = {1, 2, 5/2, 29/10, 941/290, ...} (see A073833 and, for a constant arising from that sequence, A233770).

Links

Table of n, a(n) for n=0..86.

Example

0.67517744245855713981328562507582763368407315898905...

Crossrefs

Cf. A073833, A233770.

Sequence in context: A296480 A021601 A305328 * A100124 A355336 A190648

Adjacent sequences: A332393 A332394 A332395 * A332397 A332398 A332399

Keyword

nonn,cons

Author

Jon E. Schoenfield, Feb 10 2020

Status

approved