**HP Prime and Casio fx-CG50 MicroPython: √(x + √(x + √(x + … )))**

**Introduction**

Made with LibreOffice Math |

The program SQRADD attempts to find a calculation to:

√(x + √(x + √(x + … )))

That is:

1. Start with a number x (any positive number, we’re assuming that we are in the Real Number domain)

2. Take the square root

3. Add x to the result

4. Take the square root of the result

5. Repeat steps 3 and 4 forever. (or until a condition is satisfied, such as machine accuracy is reached)

**HP Prime Program SQRADD**

The loop repeats until the difference between two terms is less than machine accuracy, 1E-12 (since HP Prime has floating point accuracy of 12 digits). Output: {x, answer, number of iterations it took to get machine accuracy}

EXPORT SQRADD(X)

BEGIN

// 2019-04-28 EWS

// √(X+√(X+√(X+…)))

// setup

LOCAL T,U,N:=1;

T:=√X;

U:=√(X+T);

// loop

WHILE ABS(U-T)≥1E−12 DO

T:=U;

U:=√(X+T);

N:=N+1;

END;

// display results {X,U,N}

RETURN {X,U,N};

END;

**Casio fx-CG50 MicroPython Script sqradd.py**

Use the Python mode of the Casio fx-CG50 (software versions 3.2 or greater). For more information of this procedure, please see: http://edspi31415.blogspot.com/2018/10/casio-fx-cg50-version-320-python-comes.html

# 2019-04-28 EWS

import math

# sqrt(X+sqrt(X+…))

# declare variables

# as float

N=1.0

X=float(input(“X=”))

T=float(math.sqrt(X))

U=float(math.sqrt(X+T))

# loop

# Casio Micropython

# has precision of 15

while abs(U-T)>=1e-15:

T=U

U=math.sqrt(X+T)

N+=1

# results

print(“X= “+str(X))

print(“Result= “)

print(str(U))

print(“Iterations= “+str(N))

**Procedure on a Four-Function or Desktop Calculator**

1. Clear the memory (usually [ MRC ] [ MRC ] )

2. Enter x and store it in memory (x [ M+ ] )

3. Take the square root [ √ ]

4. Repeat this loop (until machine accuracy is reached, that is, the the result doesn’t “change”): [ + ] [ MR ] [ = ] [ √ ]

**Results**

HP Prime

x = 1, Result: 1.61803398875, 23 iterations

x = 2, Result: 2, 20 iterations

x = 3, Result: 2.30277563773, 18 iterations

x = 4, Result: 2.56155281281, 17 iterations

x = 5, Result: 2.79128784748, 17 iterations

x = 6, Result: 3, 16 iterations

x = 7, Result: 3.19258240357, 16 iterations

x = 8, Result: 3.37228132327, 15 iterations

x = 9, Result: 3.54138126514, 14 iterations

x = 10, Result: 3.70156211871, 14 iterations

x = 11, Result: 3.85410196624, 13 iterations

x = 12, Result: 4, 14 iterations

x = 20, Result: 5, 12 iterations

x = 30, Result: 6, 12 iterations

x = 42, Result: 7, 11 iterations

x = 56, Result: 8, 11 iterations

Casio fx-CG50 MicroPython

x = 1, Result: 1.618033988749895, 30 iterations

x = 2, Result: 2, 26 iterations

x = 3, Result: 2.302775637731994, 24 iterations

x = 4, Result: 2.56155281280883, 22 iterations

x = 5, Result: 2.79128784747792, 21 iterations

x = 6, Result: 3, 20 iterations

x = 7, Result: 3.192582403567252, 20 iterations

x = 8, Result: 3.372281323269014, 19 iterations

x = 9, Result: 3.54138126514911, 19 iterations

x = 10, Result: 3.701562118716424, 18 iterations

x = 11, Result: 3.854101966249685, 18 iterations

x = 12, Result: 4, 18 iterations

x = 20, Result: 5, 16 iterations

x = 30, Result: 6, 15 iterations

x = 42, Result: 7, 14 iterations

x = 56, Result: 8, 14 iterations

**Some Interesting Observations**

When x = 1, the result is the Golden Ratio, Φ = (1 + √5)/2.

When x = n * (n-1) where n is an integer and n ≥ 2, the result is n.

Examples:

x = 2 = 2 * 1, result: 2

x = 12 = 4 * 3, result: 4

x = 56 = 8 * 7, result: 8

x = 240 = 16 * 15, result: 16

x = 992 = 32 * 31, result: 32

Fun stuff. What other patterns can we find with this or other similar calculation? Until next time,

Eddie

All original content copyright, © 2011-2019. Edward Shore. Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited. This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author.