Happy Pi Day!

[DM48]

Gerson, below is part of a python spigot algo my son uses to calculate PI. The range 10,000,000 would need to be reduced a lot for the HP48G. Any thoughts on the best way to program this up in RPL?


q, r, t, k, m, x = 1, 0, 1, 1, 3, 3
for j in range(10,000,000):
if 4 * q + r - t < m * t:
yield m
q, r, t, k, m, x = 10*q, 10*(r-m*t), t, k, (10*(3*q+r))//t - 10*m, x
else:
q, r, t, k, m, x = q*k, (2*q+r)*x, t*x, k+1, (q*(7*k+2)+r*x)//(t*x), x+2