martes, 30 de marzo de 2010

Oblongos y pitagóricos (3)

La ecuación x2+(x+1)2 =y2 se puede desarrollar de esta forma: x2+(x+1)2 =y2; 2x2+2x+1=y2; (2x+1)2 + 1 = 2y2 ; (2x+1)2 - 2y2 = -1, por lo que llamando z=2x+1 desembocamos en una ecuación de Pell con segundo miembro igual a -1

z2-2y2 = -1

Utilizamos las hojas de cálculo que presentamos en una entrada anterior

http://www.hojamat.es/sindecimales/aritmetica/herramientas/hoja/pell.ods para Calc
http://www.hojamat.es/sindecimales/aritmetica/herramientas/hoja/pell.xls para Excel

con el resultado que indica la imagen:



en la que valdrán las soluciones correspondientes a -1

Z=1; Y=1; Imposible, pues X sería negativo
Z=7; Y=5 X=3; X+1=4; Y=5
Z=41; Y=29 X=20; X+1=21;Y=29
Z=239;Y=169 X=119; X+1=120; Y=169
Z=1393;Y=985 X=696;X+1=697; Y=985

que coinciden con las propuestas en la anterior entrada dedicada a este tema.

Este método tiene el inconveniente de que depende de la precisión que tenga la hoja de cálculo en los datos con coma flotante, lo que hará que se rompa en algún momento la periodicidad de los cocientes, en este caso el 2. Por ello se puede completar con una fórmula recursiva que obtenga soluciones exactas conociendo las primeras.

En este ejemplo cada elemento de las distintas celdas cumple la fórmula

an+2 = 2an+1 + an

pero como las soluciones aparecen de forma alternada, deberemos reiterar dos veces, y nos quedará:

an+4 = 2an+3 + an+2 = 2(2an+2 + an+1)+ 2an+1 + an = 4an+2 + 4an+1+ an = 6an+2 - an

Con esta fórmula recursiva se van obteniendo las soluciones sin errores a partir de las dos primeras:

Z0 = 1; Z2 = 7; Z2 = 6*7-1 = 41; Z2 = 6*41-7 =239; …
Y0 = 1; Y2 = 5; Y4 = 6*5-1 = 29; Y6 = 6*29-5 =169; …

Pero no olvidemos que Z es una variable auxiliar Z=2X+1 y que después debemos despejar X

La siguiente lista de ternas, que coincide con la primera que propuso Girard, se ha obtenido mediante esta técnica

1 0 1
5 3 4
29 20 21
169 119 120
985 696 697
5741 4059 4060
33461 23660 23661
195025 137903 137904
1136689 803760 803761
6625109 4684659 4684660
38613965 27304196 27304197
225058681 159140519 159140520
1311738121 927538920 927538921
7645370045 5406093003 5406093004
44560482149 31509019100 31509019101
259717522849 183648021599 183648021600
1513744654945 1070379110496 1070379110497
8822750406821 6238626641379 6238626641380
51422757785981 36361380737780 36361380737781
299713796309065 211929657785303 211929657785304


Un reto: Fermat propuso una fórmula de recurrencia para generar ternas de este tipo a partir de otras similares. Dada la terna (x,x+1,y), se puede generar otra similar (x’,x’+1,y’)
mediante las fórmulas x’=2x+3y+1 y’=4x+3y+2.

¿Sabrías demostrarlo? ¿Engendra todas las ternas posibles a partir de 3,4,5?

viernes, 26 de marzo de 2010

Oblongos y pitagóricos (2)

¿Cómo organizar una búsqueda de soluciones para x2+(x+1)2 =y2 con hoja de cálculo?

Si no lo pensaste al leer la entrada anterior, ahí van dos propuestas:

Elemental

Rellena una columna con los primeros números naturales consecutivos, y en la columna de su derecha auméntalos en una unidad. Supongamos que has comenzado en las celdas B4 y C4 respectivamente. En ese caso puedes rellenar la celda D4 con la fórmula B4^2+C4^2, y en la E4 una condición que nos devuelva la palabra “Vale” si es cuadrado perfecto:

SI(D4=(ENTERO(RAÍZ(D4))^2; “Vale”;””).

Si usas Excel suprime la tilde de la palabra RAIZ. De esta forma descubriremos las soluciones, con algo de paciencia, tiempo y muchas filas de hoja de cálculo:



Con Basic

La misma idea de construir una lista para X, otra para X+1 y una tercera en la que buscamos los cuadrados perfectos se puede construir en Basic. X lo almacenamos en la variable i, X+1 en la j, y la hipotenusa en k. Una sentencia IF nos presenta las soluciones en las que k es un entero.

Con este código se buscan las soluciones para números inferiores a 1000000.

Sub busquedas
Dim i,j,k
for i=1 to 1000000
j=i+1
k=i^2+j^2
if k=Int(sqr(k))^2 then
msgbox(i)  
msgbox(j)  
msgbox(sqr(k))  
end if
next i  
End Sub

Otro día nos pondremos más serios e intentaremos un estudio algebraico. Si te quieres adelantar…

lunes, 22 de marzo de 2010

Oblongos y pitagóricos (1)

Una cuestión que ha dado juego desde los tiempos de Girard y Fermat y que permite recorrer alternativas de cálculo es la siguiente:

De todos los triángulos rectángulos de lados enteros ¿Cuáles cumplen que la diferencia entre los catetos es la unidad?

El primero que la cumple es el popular de lados 3, 4 y 5, pues la diferencia entre 3 y 4 es una unidad. ¿Habrá más? ¿Cómo abordamos el cálculo?

Casi todos los caminos nos llevan a una ecuación diofántica de segundo grado, pero hay que ver cuál y cómo resolverla. También podemos intentar una búsqueda con hoja de cálculo.

Quizás fuera prudente comenzar con esta última posibilidad. Si lo intentas descubrirás al menos estas soluciones:

3, 4 y 5
20, 21 y 29
119, 120 y 169
696, 697 y 985
4059, 4060 y 5741
23660,23661 y 33461
137903, 137904 y 195025

En la siguiente entrada propondremos la organización de una posible búsqueda por si no se te ocurre nada.