Ya que nadie se ha manifestado interesado en el método, muy interesante por cierto, aquí se los detallo para que lo prueben una o dos veces ... sólo como curiosidad intelectual, claro

Se quiere calcular, digamos, SQR(311553)
La fórmula general es: Raíz(n) = (311553 / Raiz(n-1) + Raiz(n-2) ) / 2
La raiz exacta es 558.169...
Lo primero es aproximar una primera raíz como la mitad (311553 / 2 = 155776). Esto es un buen punto de partida, ya que es uno de los peores casos. Sería mucho mejor (más rápido) aproximar la primera vez con algo más parecido a la raiz, digamos 500.
Se hace el primer ciclo de cálculo: 311553 / 155776 = 2 ; ( 2 + 155776 ) / 2 = 77889 --> MC y M+
El segundo: 311553 / MR(77889) = 3.99996148 ; ( 3.9996148 + MR(77889) ) / 2 = 38946.4999 --> MC y M+
El tercero: 311553 / MR(38946.4999) = 7.999512 ; ( 7.999512 + MR(38946.4999) ) / 2 = 19477.2497 --> MC y M+
El cuarto:311553 / MR(19477.2497) = 15.995738 ; ( 15.995738 + MR(19477.2497) ) / 2 = 9746.6227 --> MC y M+
El quinto: 311553 / MR = 31.965 ; ( 31.965 + MR ) / 2 = 4889.293 --> MC y M+
El sexto da 2476.515
El septimo 1301.159
El octavo 770.30009
El noveno 587.3785
El decimo 558.8956
El onceavo 558.16980 ... este resultado ya tiene una diferencia con la raiz real de 5 diezmilésimas
El doceavo 558.169329 ... y en este, la diferencia es de 1x10 a la -9
Si en lugar de comenzar con la mitad como primera aproximación hubiéramos comenzado con 500:
1) 623.106 diferencia con la raiz real de 64.9
2) 561.553 3.38
3) 558.1795234 0.0102
4) 558.169329235 0.000000089
Dicho de otra forma. el par de sucesionas monótonas de valores convergentes (en otros paises llamadas series telescópicas) creado por el algoritmo tiene una tendencia mucho más rápida al valor real si la primera aproximación es mejor.
En fin, esto no pretende ser una lección de mates, sino simplemente una curiosidad antediluviana necesaria cuando había que calcular una frecuencia de resonancia con calculadoraas muuuy antiguas.
Algorítmicos saludos