I'll grant you that we both started accidentally started taking on water, but that's where the similarities end!

The death toll from me wetting myself has been remarkably low, thank God 🤭

I got one today too😭

Ah shite...there is another typo. The square was on the wrong bracket.

Looks like it's actually approaching pi. Oops!

tn - pn*(an - an_1(an, bn))**2

n = 0
func_iter = []
while n < 1000:
n += 1
func_iter.append(f(an, bn, tn)) # tn should come first
tn = tn_1(tn, pn, an, bn)
a0 = an # so that bn can be updated with the current "a_n" value, not a_{n+1}
an = an_1(an, bn)
bn = bn_1(a0, bn)
pn = pn_1(pn)

def f(an, bn, tn):
return ((an + bn)**2)/(4*tn)
def an_1(an, bn):
return (an + bn)/2
def bn_1(an, bn):
return np.sqrt(an*bn)
def tn_1(tn, pn, an, bn):
return tn - pn*(an - an_1(an, bn)**2)
def pn_1(pn):
return 2*pn

import numpy as np
import matplotlib.pyplot as plt
a0 = 1
b0 = 1/np.sqrt(2)
t0 = 1/4
p0 = 1
an = a0
bn = b0
tn = t0
pn = p0

Oops, posted this and immediately noticed I had a typo in my code... Maybe there's more, I wrote this pretty quick. Someone should check my work :3

Result is the same, just changed the convergence speed. Seems like it's approaching 0 though!