ENSAYOS Y TEMAS ESPECIALES

Ecuaciones no lineales en física y su resolución mediante el uso de métodos iterativos multipaso de orden alto

Nonlinear equations in Physics, and their resolution using high order multi-step iterative methods

 

Santiago Quinga¹

 

1 departamento de Ciencias Exactas, Universidad de las Fuerzas Armadas ESPE, S/N y Ambato, Av. Gral. Rumiñahui, Sangolquí, Ecuador.

E-mail: sant.david@hotmail.com

Recibido el 30 de Abril de 2021 | Aceptado el 6 de Agosto de 2021

 

Resumen

El presente trabajo proporciona al docente de física fundamentos teóricos y computacionales para resolver ecuaciones no lineales, muy comunes en la solución de problemas físicos. En el presente trabajo de investigación se resuelven tres problemas de física, los cuales son: una esfera flotando en agua, caída no libre de un paracaidista, compresión de un resorte real; haciendo uso de principios referentes a fluidos, cinemática y dinámica. Se obtienen ecuaciones no lineales difíciles y en algunos casos imposibles de ser resueltas mediante métodos analíticos. Para encontrar una solución aproximada a dichas ecuaciones se hace uso de métodos iterativos partiendo desde los métodos tradicionales como son Newton, Secante, Steffensen hasta la introducción de métodos multipaso con alto orden de convergencia como son Traub, Ostrowski y métodos de orden ocho diseñados a partir del método de Ostrowski. Finalmente, se realiza un análisis de los resultados obtenidos al aplicar todos estos métodos a cada uno de los problemas físicos seleccionados y de esta formar establecer qué método iterativo es más adecuado ante cada situación.

Palabras clave: Ecuaciones no lineales; Métodos iterativos multipaso; Newton; Métodos de Ostrowski; Física.

 

Abstract

This work provides the physics teacher with theoretical and computational foundations to solve nonlinear equations, very common in solving physical problems. In the present research three physics problems are solved, which are: a sphere floating in water, non-free fall of a parachutist, compression of a real spring; making use of principles related to fluids, kinematics and dynamics. Nonlinear equations are obtained which are difficult and, in some cases, impossible to be solved by means of analytical methods. To find an approxi-mate solution to these equations we use iterative methods starting from traditional methods such as Newton, Secant, Steffensen to the introduction of multi-step methods with high order of convergence such as Traub, Ostrowski and methods of order eight designed from Ostrowski's method. Finally, an analysis of the results obtained by applying all these methods to each of the selected physical problems is carried out and, in this way, establish which iterative method is more appropriate in each situation.

Keywords: Nonlinear equations; Iterative multi-step methods; Newton; Ostrowski's method; Physics.

 

I. INTRODUCCIÓN

El presente trabajo proporciona al docente de física fundamentos teóricos y computacionales para resolver ecuaciones no lineales, muy comunes en la solución de problemas físicos; haciendo uso de principios referentes a fluidos, cinemática y dinámica se han planteado tres problemas cuyo modelamiento matemático lleva a la obtención de una ecuación no lineal, dichos ejemplos son: una esfera flotando en agua, caída no libre de un paracaidista, compresión de un resorte real. Al momento de encontrar la solución de una ecuación no lineal, el método que, por excelencia es utilizado por el docente, es el de Newton, sin embargo, este método no siempre es el más óptimo ante cada situación, hecho que se va a demostrar mediante un análisis de diferentes métodos iterativos. En forma general, una ecuación no lineal se puede escribir de la forma f(_x) = 0, y el objetivo será encontrar una solución de la misma, donde f: I c R ^ R es una función real que está definida en un intervalo I abierto.

El desarrollo de los métodos numéricos ha permitido la aparición de técnicas computacionales que permiten resolver varios tipos de modelos matemáticos que no tienen una solución exacta, tal es el caso de los métodos iterativos los cuales nos van a permitir obtener la solución aproximada de una ecuación no lineal.

II. MARCO TEORICO

A. Métodos iterativos

Como mencionan Mathews & Fink (2000), una iteración consiste en realizar un proceso por repetidas ocasiones hasta que se obtiene un resultado con el menor error posible. En este sentido los métodos iterativos se utilizan para hallar raíces de ecuaciones y sistemas de ecuaciones sean estos lineales o no lineales, y en la solución numérica de ecuaciones diferenciales.

Para resolver una ecuación no lineal mediante el uso de métodos iterativos se va a realizar un proceso en el cual mediante una fórmula (función de punto fijo) se va a sustituir una primera aproximación de la raíz, con ello se obtendrá una nueva estimación de la raíz, este proceso se realizará varias veces obteniéndose una sucesión de valores (iterados) que se aproximarán cada vez más a la raíz de la ecuación estudiada (García, 2017).

García (2017) indica que un método iterativo de punto fijo se representa de la siguiente forma:

^{Xn₊1 = G(Xn X    n = 0¿2,...,    (1)}

donde G es la función de punto fijo y a al ser la raíz de la ecuación debe cumplir que a = G(a) convirtiéndose de esta forma en un punto fijo de la función G ya que los puntos fijos son aquellos valores de x que satisfacen la ecuación

G( x) = x.

B. Clasificación de métodos iterativos

Según Chicharro (2017), los métodos iterativos podrían clasificarse como métodos con memoria o sin memoria en función de los datos necesarios para obtener el valor del siguiente iterado; un método sin memoria solo depende del valor del último iterado:

^{Xn_{+1 =^(Xn);    n = 0,U,...,    (2)}}

mientras que un método con memoria depende de varios iterados anteriores:

^Xn+1 =t^(Xn > ^Xn-1> ^Xn-2^,'-^{)    n} = ⁰>U,---,    ⁽³⁾

Amorós (2020) menciona que los métodos iterativos se pueden clasificar según el número de pasos en cada iteración, obteniéndose de este modo métodos de un paso o multipaso; un método de un paso tiene la forma de la ecuación (1), mientras que un método multipaso viene descrito por:

y_n = w^(x„⁾

^Xn+1 =^^{(Xn > yn X n = °X2,---    (4)}

El ultimo criterio que se puede emplear para clasificar un método iterativo es dependiendo de la presencia o ausencia de derivadas, obteniéndose métodos con derivadas o libres de derivadas (Chicharro, 2017).

C. Orden de convergencia

Sea {_Xjt} _> una sucesión de números reales tales que n_{m x} = g. La convergencia se llama:

a)    Lineal, si existe M, 0 <M < 1 y k₀ e N tales que _^<M\x_k -g¡, para todo k > k_{0 ;}

b)    De orden p, si existe M > 0 y k₀ e N tales que    _^<M\x_k-g¡^p, para todo k > k₀.

Denotamos por _e<_ = _Xk _% el error de la k - ésima iteración. La ecuación ^ +1 = M_e^p + o(e^p+1), se llama ecuación del error, donde p es el orden de la convergencia (Artidiello, 2014).

Del orden de convergencia de un método iterativo depende el número de iteraciones que se debe realizar para obtener una solución aproximada de la ecuación no lineal, siendo los métodos de orden de convergencia alto quienes se aproximan más rápido a la solución, necesitando de este modo menos iteraciones (Amorós, 2020).

Para realizar una comparación entre métodos iterativos se pude utilizar parámetros relacionados al orden de convergencia tales como el índice de eficiencia:

I = P^{d    (5)}

donde p es el orden de convergencia del método y d es el número de evaluaciones de la función f y sus derivadas, realizadas en cada iteración (Ostrowski, 1966).

Por otro lado, Kung y Traub (1974) establecen una conjetura que permite saber si un método es óptimo o no, dicha conjetura establece que el orden de convergencia p de un método sin memoria con d evaluaciones funcionales por iteración cumple con:

p < 2^{d-1    (6)}

siendo un método óptimo aquel que cumple la igualdad en la expresión anterior.

D. Método de Newton

Como indica Ypma (1995) el método de Newton es uno de los más conocidos y utilizados al momento de resolver una ecuación no lineal, en donde la función de punto fijo tiene la siguiente forma:

^X„₊1 = ^Xn

^{f ( Xn ) .}

f( Xn )’

n = 0,1,2,...,

(7)

en ocasiones también es nombrado como método de Newton-Raphson.

Además, este método tiene una representación geométrica en donde, en cada etapa n, el valor _X corresponde

a la abscisa del punto de corte con el eje OX de la recta tangente a la curva y = f (_X) en el punto X; f (_Xn)), debido

a ello este método también es conocido como el método de la tangente (Artidiello, 2014).

Como menciona Amorós (2020) en base a la teoría planteada anteriormente, el método de Newton tiene las siguientes características:

a)    Es un método sin memoria.

b)    Es un método de un paso.

c)    Es un método con derivada

d)    Tiene un orden de convergencia cuadrática

e)    Según la conjetura de Kung y Traub es un método óptimo

E. Método de la secante

Cordero, Hueso, y Torregrosa (2006) indican que en el método de Newton se necesitan realizar evaluaciones funcionales tanto de la función como de su derivada en cada paso, en ocasiones podría resultar un poco complicada la

obtención de la derivada de la función. El método de la secante no necesita la derivada para el cálculo de la estimación, pero a diferencia del método de Newton su orden de convergencia es menor, la función de punto fijo a utilizarse es:

x

H+1

= x

^{f (x}n ⁾

^{f (X}„ ) - ^{f (X}„-l)

^Xn ^{- X}n-1

n = 0,1,2,...

(8)

El método de la secante tiene las siguientes características:

a)    Necesita dos estimaciones iniciales.

b)    Es un método con memoria.

c)    Es un método libre de derivadas.

d)    Es un método de un paso.

e)    El orden de convergencia aproximado es 1,618033... (Mathews & Fink, 2000).

f)    Al ser un método con memoria no se puede analizar si es óptimo con la conjetura de Kung y Traub.

F. Método de Steffensen

Steffensen (1933) menciona que partiendo del método de Newton descrito en la ecuación (7) y sustituyendo la derivada f '(_x ) por su diferencia finita progresiva se obtiene un método con orden de convergencia cuadrática y cuya

expresión es la siguiente:

^xn+1 = ^xn

f (^x„ )²    .

^f i^x„ + ^{f (x}n ))^{- f ( x}n )’

n = 0,1,2,...

(9)

García (2017) menciona algunas de las características del método de Steffensen como son:

a)    Es un método sin memoria.

b)    Es un método de un paso.

c)    Es un método libre de derivadas.

d)    Tiene un orden de convergencia cuadrática

e)    Según la conjetura de Kung y Traub es un método óptimo

G. Métodos Multipaso de orden alto

Con el objetivo de encontrar métodos iterativos que tengan un mayor orden de convergencia, surgen métodos que requieren, en forma general, un aumento en el número de evaluaciones funcionales, así como del número de pasos, obteniéndose de esta forma métodos multipaso cuya forma es la expresada en la ecuación (4). A continuación, se presentan dos métodos iterativos multipaso que se van a utilizar en esta investigación.

El método de Traub (1964) también conocido como el método de Potra-Pták es un método de tercer orden que se obtiene al realizar una composición entre el método de Newton cuya expresión se encuentra en la ecuación (7) consigo mismo y realizando una sustitución en el segundo paso por la derivada empleada en el primer paso, su esquema es el siguiente:

yn = ^xn -

^xn+1 = y„

^{f ( xn )}

f(x„)

_- ^{f (y»)}

f '( x„)’

n = 0,1,2,.

(10)

Como menciona García (2017) el método de Traub no es óptimo.

Ostrowski (1960) Diseñó el primer método iterativo cuyo orden de convergencia es cuatro, para ello hizo uso de un esquema multipaso y, además, utilizó técnicas de interpolación, la expresión obtenida es la siguiente: ^{f (Xn )}

y = X

n    n

f( Xn )

x

n+1

^{f (X}n ^{)    f} (y«^{) .}

f (Xn ) - 2f (yn ) f'( Xn )’

n = 0,1,2,...
(11)

Como indica (Chicharro, 2017), el método de Ostrowski es un método óptimo.

Como indican Dzunic y Petkovic (2012), partiendo del esquema de dos pasos de Ostrowski se puede llegar a obtener métodos iterativos de tres pasos de orden ocho. El método óptimo de orden ocho propuesto por los mismos es el siguiente:

y = X

n    n

z = y

nn

X

n+1

=z

n

^{f (X}n ⁾

f( Xn )

^{f (X}n ^{)    f (}y„⁾

f (Xn) - 2f (yn) f'(Xn)

f (z„) 1 + f (z_n)/f (y_n)1 + 2f (z_n)/f (X_n). f(X„) 1 -2f (y„)/f (X_n)-(f (y„)/f (X_n))²â€™

n = 0,1,2,.
(12)

Cordero, Torregrosa y Vassileva (2011) diseñan un método iterativo de orden ocho óptimo, partiendo del esquema de Ostrowski y cuyo esquema es el siguiente:

y = X

n    n

^{f (X}n ⁾

f'(Xn )

^{f ( X}n ^{) f (X}n ^{) - f (}yn ⁾

^Xn+1 = ^Un ^{- 3}

f \x« ) f (Xn ) - 2f (yn ) f (z ) u - z

nnn

(13)

^f' ^(Xn ⁾ yn ^{- X}n

n = 0,1,2,.,

donde

 

+ 1    j2 f (yn ) - 2f ( z« ) J
III. METODOLOGIA

Referentes a física, en el presente trabajo de investigación se plantean los tres problemas siguientes:

a)    Esfera flotando en agua

b)    Caída no libre de un paracaidista

c)    Compresión de un resorte real

Cuyo análisis lleva a la obtención de una ecuación no lineal, y cuya solución se va a obtener mediante diferentes métodos iterativos para luego realizar un análisis entre los resultados obtenidos por cada uno de estos métodos, y de esta forma poder concluir cual método responde de mejor manera ante cada uno los problemas planteados, los métodos iterativos que se van a emplear son:

a)    Método de la Secante

b)    Método de Newton

c)    Método de Steffensen

d)    Método de Traub

e)    Método de Ostrowski

f)    Método de Dzunic y Petkovic

g)    Método de Cordero, Torregrosa y Vassileva

cuyas expresiones iterativas fueron descritas en la sección anterior, y cuya programación se realizó en MATLAB, los códigos se pueden revisar en anexos.

A.    Proceso

A continuación, se detalla el proceso que se va a realizar para resolver cada uno de los problemas planteados:

1)    Aplicación de los fundamentos físicos y matemáticos para obtener una ecuación no lineal de la forma f (x) = 0, donde x será la variable física en estudio.

2)    Obtención de una gráfica de la función f (_x), para de este modo observar los cortes con el eje x, y de esta forma poder elegir una estimación inicial lo más cercana a dicho corte, esta estimación es un parámetro de entrada para los métodos iterativos programados.

3)    Selección de valores para los parámetros de entrada de los métodos programados como son: la tolerancia, el número máximo de iteraciones que junto a la estimación inicial seleccionada en el paso anterior nos permitirán obtener la solución aproximada a la ecuación no lineal.

B.    Esfera flotando en agua

Una esfera de densidad p y de radio r está sumergida parcialmente en agua hasta una profundidad d, calcular dicha profundidad.

FIGURA 1. Se muestra una esfera de radio r, sumergida en agua una profundidad d .

Para calcular el volumen del casquete esférico sumergido en el agua necesitamos resolver la siguiente integral:

V =

d

JV(r² - (r - y)²)dy

(14)

0

la cual se obtiene mediante la utilización del método de discos, donde y es la distancia desde un disco de posición arbitraria hasta el extremo del casquete esférico, y dy representa el ancho infinitesimal de dicho disco. Además, la fórmula _r ² _-(_{r -} y)² obtenida mediante el uso del teorema de Pitágoras, nos permite calcular el radio al cuadrado del disco.

Al resolver la integral planteada en la ecuación (14), se obtiene el siguiente volumen del casquete esférico sumergido en el agua:

V = ^f (3r - d),

(15)

de esta forma el peso del agua desplazada es:

^Pa = gPa    ^{(3r - d}^

(16)

donde p_a representa la densidad del fluido, en este caso agua.

Por otro lado, la masa de la esfera se puede calcular con la siguiente formula:

m = p

e re

(17)

Como la esfera se encuentra parcialmente sumergida en el agua, sobre la superficie de la misma actúan los efectos de presión tanto del aire como del agua, dichos efectos aparecen en forma de fuerzas perpendiculares a su superficie. Si se analiza la flotación de la esfera en el agua se obtiene que el empuje experimentado por la esfera tiene que ser igual al peso de la misma, ya que esta se encuentra en equilibrio. De esta forma se puede calcular el empuje con la siguiente fórmula:

^E = m_eg = gp_e

4mr³ 3

(18)

Aplicando el principio de Arquímedes el cual establece que "Un cuerpo total o parcialmente sumergido en un fluido en reposo experimenta un empuje vertical hacia arriba igual al peso del fluido desalojado". De esta forma igualando las ecuaciones (16) y (18), obtenemos la siguiente expresión:

(19)

4mr³    md² ¡    \

gPe    = gPa    l^3- - ^d)^,

igualando a cero la expresión anterior:

p_ad³ - 3rp_ad² + 4p_er³ = 0^{,    (20)}

se obtiene una ecuación no lineal de la forma f (d) = 0, para continuar con la solución se va a considerar los siguientes datos:

Densidad del agua: p = \ g / _cm³

Radio de la esfera: r = 10 cm

Esfera de madera de densidad: p = 0.638 g / _cm ³

Utilizando la función f(d) = p_ad³ _-3-pd² + 4p _r³ y los respectivos datos considerados anteriormente, se obtiene la figura 2:

FIGURA 2. Representación gráfica de la función f (d) = p_ad³ _-3_rp d² + 4p_rr³, para los valores de los parámetros

p_a= 1 g / cm³, r = 10 cm, p_e = 0.638 g / cm³.

Como se puede apreciar en la figura 2, existe dos ceros de la función en el intervalo mostrado, tomando en consideración que la esfera está parcialmente sumergida, se va a calcular solo el valor de la raíz que está del lado izquierdo, y cuyo valor está en torno a d = 11 cm. Esta figura es de mucha utilidad para establecer el valor de la estimación inicial necesaria en los diferentes métodos iterativos.

Para resolver la ecuación no lineal (20) escogiendo una estimación inicial de la profundidad d₀ = 11 y una tolerancia 10^-10, para los valores de los parámetros p_a= 1 g/_cm³,r = 10 cm, p = 0.638 g/_cm³ con el método de Ostrowski obtenemos la siguiente profundidad aproximada d «11.8615 cm, luego de realizar 3 iteraciones.

B. Caída no libre de un paracaidista

Calcular la masa m de un paracaidista que tiene una caída no libre y tiene una rapidez v₀ cuando t = 0. Para simplificar

la modelización matemática se va a despreciar el efecto de la fuerza de rozamiento vorticoso misma que es proporcional al cuadrado de la rapidez del paracaidista, para solo considerar como fuerza de resistencia la referente al rozamiento viscoso que ofrece el aire al paracaidista la cual es proporcional a la rapidez del mismo. Otros datos necesarios para el cálculo de la masa del paracaidista son el instante de tiempo t > 0 en donde adquiere la rapidez v misma que debe ser menor a la rapidez límite del paracaidista que es cuando el mismo desciende con fuerzas en equilibrio.

i Fr r	a
i	r mg

FIGURA 3. Fuerzas que actúan sobre el paracaidista durante su caída, las mismas que son fuerza de resistencia Fr dirigida hacia arriba y el peso W = _mg dirigido hacia abajo.

Como la fuerza de rozamiento es proporcional a la rapidez del paracaidista entonces, su módulo es:

Fr = kv,    (21)

Considerando la dirección del movimiento positiva hacia abajo haciendo uso de la figura 3 y aplicando la segunda ley de Newton, obtenemos:

(22)

dv

W - Fr = m ^, dt

en donde al reemplazar la ecuación (21) y el equivalente del peso W = mg , nos encontramos con la siguiente ecuación diferencial:

(23)

dv

mg - kv = m .

dt

Al resolver por varias separables la ecuación (23), se plantea:

dt = m-

dv

(24)

mg - kv

como tiene una rapidez v₀ cuando t = 0, entonces la condición inicial es que v(t = 0) = v₀, así que obtenemos:

(25)

f,    } dv

I dt =1 m    ^,

0    Ã­ ^mg^{- kv}

al integrar los dos miembros de la ecuación:

m,

t =--ln

k

mg - kv

^mg ^{- kv}0

(26)

si despejamos v, de la ecuación (26), encontramos la fórmula que relaciona la rapidez en función del tiempo, así:

v = mg ₊

k

k \

mg —

v - m^{te m}

v ^k y

(27)

Como el objetivo es encontrar la masa del paracaidista, supuesto conocido la rapidez en un instante, igualando a cero la expresión anterior:

( - (28)

mg

k

0

mg

--e

k

m

- v = 0^,

se obtiene una ecuación no lineal de la forma f (_m) = 0, para continuar con la solución se va a considerar los siguientes datos:

Rapidez inicial: v₀ = 1 m / 5

Gravedad: g = 9.8 m / s²

Constante de proporcionalidad: k = 14 kg / s

Rapidez: v = 9.8 m / s

Tiempo: t = 8 s

Utilizando la función

^{f (m})

mg_ ₊

k

(

^v0

- ms_e

k

-v

y los respectivos datos considerados anteriormente, se ob

y

FIGURA 4. Representación gráfica de la función

- v’

tiene la figura 4:

f (m) = ^ ₊ k

0

k

para los valores de los parámetros

v₀ = 1 m /s, g = 9.8 m/s², k = 14 kg/s, v = 9.8 m/s, t = 8 s .

Como se puede apreciar en la figura anterior, existe un cero de la función en el intervalo mostrado, y cuyo valor está en torno al valor de _m = 55 kg, esta figura es de mucha utilidad para establecer el valor de la estimación inicial necesaria en los diferentes métodos iterativos.

Para resolver la ecuación no lineal (28) escogiendo una estimación inicial de la masa _m = 60 y una tolerancia 10" ¹⁰, para los valores de los parámetros v₀ = 1 m/s , g = 9.8m/s², k = 14 kg/s, v = 9.8m/s, t = 8 s con el método de Traub obtenemos la siguiente masa aproximada m « 56.1698 kg, luego de realizar 15 iteraciones.

C. Compresión de un resorte real

En los sistemas mecánicos reales se utiliza resortes no lineales, como ejemplo de aplicación se presenta el siguiente problema: un objeto de masa m se deja caer desde una altura h sobre un resorte real cuya fuerza elástica es Fe = -(k_ld + k₂d^3/2) donde d es la compresión del mismo; calcular la máxima compresión del resorte.

FIGURA 5. Se muestra un objeto de masa m que cae desde una altura h sobre un resorte real, a la izquierda el resorte está en posición normal, a la derecha el resorte está comprimido al máximo, d es la máxima compresión del resorte.

Utilizando el principio del Trabajo y Energía:

(29)

EC +    = EC₂,

como el objeto se deja caer del reposo y además se busca la máxima compresión del resorte, entonces la energía cinética del objeto es nula en las posiciones 1 y 2, además reemplazando el trabajo del peso y del resorte, tenemos:

(30)

d

WAy - J(kd + k₂d^3/2)dd = 0,

0

la variación de y se calcula como Ay = h + d además integrando el segundo término de la ecuación (27), obtenemos:

mgh + mgd — 1 k_xd² — 2 k₂d^5/2 = 0,    ⁽³¹⁾

se obtiene una ecuación no lineal de la forma f (d) = 0, para continuar con la solución se va a considerar los siguientes datos:

Gravedad: g = 9.81m/s²

Constantes de proporcionalidad: k = 40000g / s², k₂ = 40g /{s²m⁰'⁵) Masa del objeto : _m = 95 g Altura: h = 0.43 m

Utilizando la función

2.

f (d) = mgh + mgd — k_xd — k₂d

5/2

2

5

y los respectivos datos considerados anteriormente, se

obtiene la siguiente figura:

FIGURA 6. Representación gráfica de la función f (d) = g = 9.81 m/s², k = 40000 g/s², k₂ = 40 g/(s²m⁰⁵),

: mgh + mgd — — k_xd² — — k₂d⁵

, para los valores de los parámetros

m = 95 g, h = 0.43 m.

Como se puede apreciar en la figura anterior, existe un cero de la función en el intervalo mostrado, y cuyo valor está en torno al valor de d = 0.2 m, esta figura es de mucha utilidad para establecer el valor de la estimación inicial necesaria en los diferentes métodos iterativos.

Para resolver la ecuación no lineal (31) escogiendo dos estimaciones iniciales de la máxima compresión d = 0.2, d = 0.1 y una tolerancia 10^-10, para los valores de los parámetros g = 9.81 m / s², k = 40000 g / s², k₂ = 40 g l{s²m°), _m = 95 g, con el método de la Secante obtenemos la siguiente máxima compresión aproximada d « 0.1667 m, luego de realizar 6 iteraciones.

IV. RESULTADOS

Para las pruebas respectivas se ha utilizado MATLAB R2015a, el sistema computacional utilizado para realizar todas las pruebas tiene las siguientes especificaciones: Intel(R) Í5-420U CPU @ 1.70GHz 2.40GHz, memoria RAM de 6 GB.

Los parámetros de entrada necesarios para los programas de los métodos iterativos utilizados en la presente investigación son: la función respectiva de cada uno de los problemas físicos analizados, la estimación inicial * de la raíz

buscada, la tolerancia que además sirve como criterio de parada y cuya representación es _x | < 10_¹⁰, es decir,

estamos trabajando con una tolerancia de 10 _¹⁰ , otro parámetro de entrada es el número máximo de iteraciones, ya que se podría producir un desborde de iteraciones al no alcanzar la tolerancia establecida; dichos parámetros de entrada se los puede apreciar en los códigos de cada programa presentados en anexos.

TABLA I. Resultados obtenidos al utilizar los métodos de la Secante, Newton, Steffensen, Traub y Ostrowski para resolver tres problemas físicos la esfera flotando en agua, la caída no libre de un paracaidista y la compresión de un resorte real.

			Método	P	z	K+1 _ Xn\	Iter	Tiempo
			Secante	1.6180	11.8615	1.776E-15	6	0.00552200
			Newton	2	11.8615	6.573E-14	4	0.00534150
Esfera flotando en agua Pad3 _ 3rPad2 + 4Per' = 0			Steffensen	2	---	---	---	---
	d„ = 11		Traub	3	11.8615	1.066E-14	3	0.00527680
	d * = 11		Ostrowski	4	11.8615	1.776E-15	3	0,00534590
	II		Dzunic y Petkovic	8	11.8615	8.882E-15	3	0.00447070
			Cordero, Torregrosa y Vassileva	8	11.8615	1.776E-15	3	0.00441745
			Secante	1.6180	56.1698	3.2188E-12	5	0.00483015
Caída no libre de un paracaidista			Newton	2	56.1698	6.1036E-11	28	0.00478300
mg	1 _ e k	_ v = 0	Steffensen	2	---	---	---	---
k	^    J		Traub	3	56.1698	3.4710E-11	15	0.00476565
	O o II o s		Ostrowski	4
	3 ll k/\ L/t
	3 ll Q\ O		Dzunic y Petkovic	8	---	---	---	---
			Cordero, Torregrosa	8
			y Vassileva
			Secante	1.6180	0.1667	1.002E-11	6	0.00536435
Compresión de un resorte real			Newton	2	0.1667	2.776E-17	5	0.00510065
mgh + mgd	_1 k,d2 _ -kd5/2 = 0 2 1    5 2 d0 = 0.2		Steffensen	2	0.1667	6.500E-12	682	0.01932860
			Traub	3	0.1667	0.000E+00	4	0.00524890
	II O		Ostrowski	4	0.1667	2.776E-17	3	0.00491215
	d2 = 0.2		Dzunic y Petkovic	8	---	---	---	---
			Cordero, Torregrosa	8
			y Vassileva

* Estimaciones iniciales para el caso del método de la Secante.

En la tabla I se muestran los resultados obtenidos al emplear los métodos iterativos de la Secante, Newton, Stef-fensen, Traub, Ostrowski, Dzunic y Petkovic y Cordero, Torregrosa y Vassileva; para resolver tres problemas físicos; encontramos la estimación inicial * utilizada para alcanzar la raíz £, el orden de convergencia p, el error cometido

en la última iteración |^_i __x |, el número de iteraciones (Iter) necesarias para cumplir con la tolerancia del criterio de

parada, así como el tiempo necesario en cómputo (tiempo) el cual se calculó mediante el promedio de 20 ejecuciones consecutivas de cada programa.

Para los parámetros de entrada como son la tolerancia menor a 10 _¹⁰ y las respectivas estimaciones iniciales de la raíz para cada problema físico, se producen problemas con la utilización de algunos métodos ya que no son capaces de resolver las diferentes ecuaciones no lineales utilizadas; como ejemplo el método de Steffensen no consigue la solución del problema de la esfera flotando en agua, los dos métodos de orden ocho no consiguen resolver los problemas de la caída no libre de un paracaidista y el de la compresión real de un resorte.

Para el caso de la esfera flotando en agua, tomando como criterios de análisis el tiempo computacional requerido, el número de iteraciones y el error, se observa que el método de orden 8 de Cordero, Torregrosa y Vassileva es el mejor de todos los métodos utilizados, recordando además que es un método óptimo.

Para el caso de la caída no libre de un paracaidista, tomando como criterios de análisis el tiempo computacional requerido, el número de iteraciones y el error, se observa que el método de orden 3 de Traub es el mejor de todos los métodos utilizados, recordando además que es un método no óptimo.

Finalmente, para el caso de la compresión de un resorte real, tomando como criterios de análisis el tiempo computacional requerido, el número de iteraciones y el error, se observa que el método de orden 3 de Traub es el mejor de todos los métodos utilizados, recordando además que es un método no óptimo.

V. CONCLUSIONES

Si bien al momento de resolver una ecuación no lineal el método por excelencia a ser utilizado es Newton, en el presente trabajo se ha demostrado que existen otros métodos que ofrecen mejores características ya sea por la estructura de los mismos como es el caso de aquellos métodos libres de derivadas, así como por el orden de convergencia alto.

Al utilizar métodos iterativos con mayor orden de convergencia se reduce el número de iteraciones requeridas en la obtención de la solución aproximada de la ecuación no lineal estudiada, reduciendo de esta forma el tiempo de ejecución siendo más eficientes que los métodos tradicionales de segundo orden como son Newton o Steffensen.

Para los problemas planteados el método de Newton siempre converge a la solución, convirtiéndose de esta forma en el método más estable empleado en nuestra investigación.

Al utilizar funciones un tanto más complejas como en el caso de esta investigación la caída no libre de un paracaidista que tiene una naturaleza exponencial, y la compresión real de un resorte que tiene términos irracionales los métodos de orden ocho utilizados en esta investigación no siempre convergen a la solución, por lo tanto, son métodos inestables.

Para entender de mejor forma los criterios de estabilidad de un método que hacen que algunos métodos no convergen a la solución, es necesario realizar un estudio dinámico de los mismos, para ello existe diferentes investigaciones que pueden ser revisadas con este propósito.

Los métodos iterativos de punto fijo tienen una amplia gama de aplicación en diferentes campos del conocimiento, en todos aquellos en donde la modelación matemática lleve a la obtención de ecuaciones no lineales.

En próximas investigaciones se podría mejorar los fenómenos físicos estudiados, como ejemplo para la caída no libre del paracaidista se podría utilizar un modelo más real para la fuerza de rozamiento, en donde la misma sea proporcional a la rapidez al cuadrado.

 

REFERENCIAS

Amorós, C. (2020). Estudio sobre convergencia y dinámica de los métodos de Newton, Stirling y alto orden. Madrid: Universidad Internacional de la Rioja.

Artidiello, S. (2014). Diseño, implementación y convergencia de métodos iterativos para resolver ecuaciones y sistemas no lineales utilizando funciones peso. Valencia: Universitat Politécnica de Valencia.

Chicharro, F. (2017). Análisis dinámico y aplicaciones de métodos iterativos de resolución de ecuaciones no lineales. Valencia: Universitat Politécnica de Valencia.

Cordero, A., Hueso, J., & Torregrosa, J. (2006). Problemas resueltos de métodos numéricos. Madrid: Paraninfo.

Cordero, A., Torregrosa, J. R., & Vassileva, M. P. (2011). Three-step iterative methods with optimal eighth-order convergence. Journal of Computational and Applied Mathematics, 235, 3189-3194.

Dzunic, J., & Petkovic, M. S. (2012). A Family of Three-Point Methods of Ostrowski's Type for Solving Nonlinear Equations. Journal of Applied Mathematics.

García, J. (2017). Análisis dinámico. Valencia: Universidad Politécnica de Valencia.

Kung, H. T., & Traub, J. F. (1974). Optimal order of one-point and multipoint iteration. J. Assoc. Comput. Math, 21(4), 643651.

Mathews, J. H., & Fink, K. D. (2000). Métodos numéricos con Matlab. Madrid: Prentice Hall.

Ostrowski, A. M. (1960). Solution of equations ans systems of equations. New York: Academic Press.

Ostrowski, A. M. (1966). Solutions of equations and systems of equations. Academic Press.

Steffensen, J. F. (1933). Remarks on iteration. Skand. Aktuarietidskrift, 64-72.

Traub, J. F. (1964). Iterative methods for the solution of equations. Englewood Cliffs, N.J., Prentice-Hall.

Ypma, T. J. (1995). Historical development of the Newton-Raphson method. SIAM Rev., 37, 531-551.