Ejercicio: comparativa numérica de métodos iterativos

Newton, Halley, Ostrowski, Traub, punto medio, Jarratt y Newton doble sobre funciones de prueba con tolerancia 10⁻¹⁰⁰: iteraciones, residuos y ACOC confirmando los órdenes teóricos.

Condiciones del experimento

Se comparan siete métodos bajo condiciones idénticas: criterio de parada xk+1xk<10100|x_{k+1}-x_k|<10^{-100} o 60 iteraciones, con aritmética de precisión extendida (400 dígitos), imprescindible para tolerancias tan extremas. Las funciones de prueba y sus raíces son f1(x)=sinxexf_1(x)=\sin x-e^{-x} con α0.58853274\alpha\approx 0.58853274 y f3(x)=(x1)31f_3(x)=(x-1)^3-1 con α=2\alpha=2.

Resultados para f₁ (x₀=0.1)

Métodoiterf(xk+1)|f(x_{k+1})|xk+1xk|x_{k+1}-x_k|ACOC
Newton86.5531·10⁻²⁰⁵1.0865·10⁻¹⁰²2.0000
Halley605.3661·10⁻¹⁸⁷3.0000
Ostrowski51.2716·10⁻⁴⁰⁸6.7766·10⁻¹⁹⁹4.0000
Traub61.2716·10⁻⁴⁰⁸9.3924·10⁻¹⁶⁶3.0000
Punto medio61.2716·10⁻⁴⁰⁸2.9422·10⁻¹⁹²3.0000
Jarratt51.2716·10⁻⁴⁰⁸5.1327·10⁻¹⁹⁸4.0000
Newton doble51.2716·10⁻⁴⁰⁸4.7250·10⁻²⁰⁵4.0000
Resultados para f1(x)=sinxexf_1(x)=\sin x-e^{-x} con x0=0.1x_0=0.1.

Resultados para f₃ (x₀=1.5)

Métodoiterf(xk+1)|f(x_{k+1})|xk+1xk|x_{k+1}-x_k|ACOC
Newton112.8174·10⁻³⁵⁹3.0646·10⁻¹⁸⁰2.0000
Halley701.7850·10⁻²¹⁴3.0000
Ostrowski607.3471·10⁻²³⁹4.0000
Traub581.2798·10⁻³⁹³5.9750·10⁻¹³²3.0000
Punto medio72.1994·10⁻³⁹⁹9.2824·10⁻¹³⁴3.0000
Jarratt607.3471·10⁻²³⁹4.0000
Newton doble603.0646·10⁻¹⁸⁰4.0000
Resultados para f3(x)=(x1)31f_3(x)=(x-1)^3-1 con x0=1.5x_0=1.5.