Algoritmo babilónico

En este artículo o archivo se presentan y explican varios métodos que se pueden utilizar para calcular la raíz cuadrada de un número real positivo, siendo el más conocido el método de resolución.

En la imagen adjunta podemos ver cinco partes esenciales de la raíz cuadrada en este método de resolución:

Los pasos a seguir son estos:

La raíz cuadrada de 5836,369 es 76,39, con un residuo de 0,9369. Recordemos que el cero es sólo un auxiliar. Es también que la operación anterior utilizada como ejemplo no está completa. Si la continuáramos daría como resultado 76,396132101 (con nueve decimales).

Los pasos se pueden resumir en ciclos de cuatro después de separar en grupos de dos cifras y teniendo en cuenta cuando se coloca la coma decimal en la raíz:

Las calculadoras de bolsillo típicamente implementan buenas rutinas para calcular la función exponencial y el logaritmo natural; entonces calculan la raíz cuadrada de ${displaystyle x}$ utilizando la identidad

La misma identidad es usada cuando se calculan las raíces cuadradas con tablas de logaritmos o reglas de cálculo.

Muchos de los métodos de cálculo para raíces cuadradas requieren un valor inicial. Si el valor inicial está muy lejos de la raíz cuadrada real, el cálculo será muy lento. Por lo tanto es útil tener un cálculo aproximado, que puede ser muy inexacto pero fácil de calcular. Una forma de obtener tal estimación para ${displaystyle {sqrt {x}}}$ está calculando ${displaystyle 3^{D}}$, donde ${displaystyle D}$ es el número de dígitos (a la izquierda del punto decimal) de ${displaystyle x}$. Si ${displaystyle x<1}$, ${displaystyle D}$ es el negativo del número de ceros a la derecha inmediata del punto decimal.

Un mejor método de estimación es éste:

Al trabajar en el sistema de numeración binario (como lo hacen las computadoras internamente), un método alternativo es utilizar ${displaystyle 2^{leftlfloor D/2 ight floor }}$ (aquí D es el número de dígitos binarios).

El algoritmo babilónico^[2]​ se centra en el hecho de que cada lado de un cuadrado es la raíz cuadrada del área. Fue usado durante muchos años para calcular raíces cuadradas a mano debido a su gran eficacia y rapidez. Para calcular una raíz, dibuje un rectángulo cuya área sea el número al que se le busca raíz y luego aproxime la base y la altura del rectángulo hasta formar o por lo menos aproximar un cuadrado.

El algoritmo se puede enunciar sin el uso de dibujos como sigue:

Raíz(x):

Este algoritmo aproxima la raíz cuadrada de cualquier número real tanto como se desee. Es claro que no se necesita conocer el valor de ${displaystyle h}$, puesto que depende directamente de ${displaystyle x}$ y que el área del rectángulo siempre se aproxima a la raíz cuadrada de ${displaystyle x}$ sin importar el valor de ${displaystyle b}$ siempre y cuando ${displaystyle b>0}$. De esta manera surge la función recursiva

de manera tal que ${displaystyle n}$ es la ${displaystyle n}$-ésima aproximación a ${displaystyle {sqrt {x}}}$. Esto implica que

Puesto que algunas raíces son números irracionales es necesario definir qué tanto es "aproximadamente". Afortunadamente nadie es capaz de escribir un número con una infinita cantidad de dígitos, por lo que el umbral de aproximación se limita a la cantidad de dígitos que se es capaz de escribir. Entonces podemos definir que el algoritmo termine en el momento que la última aproximación es la misma que la anterior (es decir, ya no se puede aproximar más).

De manera formal, se expresa el algoritmo babilónico usando pseudocódigo de la siguiente manera:

función ${displaystyle mathrm {raiz} (x),}$

donde ${displaystyle xleftarrow y}$ significa "sustituya el valor de ${displaystyle x}$ por del de ${displaystyle y}$", y devolver expresa el resultado del algoritmo y su terminación.

Los irracionales cuadráticos (números de la forma ${displaystyle {frac {a+{sqrt {b}}}{c}}}$, donde a, b y c son enteros), y en particular, las raíces cuadradas de números enteros, tienen fracciones continuas periódicas. Podemos estar interesados a veces no en encontrar el valor numérico de una raíz cuadrada, sino por algo en su expansión como fracción continua. El algoritmo iterativo siguiente se puede utilizar para este propósito (S es cualquier número natural que no sea un cuadrado perfecto):

Hay que notar que m_n, d_n, y a_n son siempre enteros. El algoritmo termina cuando en este trío el resultado nuevo que obtenemos ya empieza a ser igual al anterior. La expansión se repetirá entonces. La secuencia [a₀; a₁, a₂, a₃, …] es la expansión fracción continua:

Comenzamos con m₀=0; d

Ahora de enlaza de nuevo con la segunda ecuación de arriba.

Por lo tanto, la fracción continua para la raíz cuadrada de 114 es: ${displaystyle {sqrt {114}}=[10;1,2,10,2,1,20,1,2,10,2,1,20,1,2,10,2,1,20,...].}$

Este método para encontrar una aproximación a la raíz cuadrada fue descrito en un manuscrito antiguo llamado manuscrito de Bakhshali. Equivale a dos iteraciones del método babilónico comenzando con el número ${displaystyle n}$ tal que ${displaystyle n^{2}}$ es el cuadrado más cercano a ${displaystyle x}$.

Si queremos calcular ${displaystyle {sqrt {10,5}}}$ con este método lo primero que hacemos es asignarle el número cuadrado perfecto cuyo cuadrado se acerque más a 10,5, ese número va a ser 3, ya que al dar ${displaystyle 3^{2},!}$ como resultado 9 se acerca más a 10,5 que ${displaystyle 4^{2},!}$ que da 16, con lo que ahora en la igualdad sustituimos:

Siendo las cifras 384615 periódicas.

Este método da un valor bastante aproximado de la raíz cuadrada del número, se puede observar también que este método al dar el resultado mediante una fracción da un número racional, mientras que la raíz cuadrada real de un número es irracional siempre que este no sea un cuadrado perfecto (o el cociente de dos cuadrados perfectos).

Si N es una aproximación a ${displaystyle {sqrt {S}}}$, una aproximación mejor puede ser encontrada usando la serie de Taylor de la función de la raíz cuadrada:

Como método iterativo, el orden de convergencia es igual al número de los términos usados. Con 2 términos, es idéntica al método babilónico; con 3 términos, cada iteración toma casi tantas operaciones como la aproximación de Bakhshali, pero converge más lentamente. Por lo tanto, esta no es una manera particularmente eficiente de la operación.

---