Ecuaciones No-Lineales

Una ecuación lineal es del tipo ax + b = 0 donde los coeficientes a, b  R . Este tipo de
ecuación puede
b
1. Tener solución única, si a  0 , x = − .
a
2. Tener infinitas soluciones, si a = 0 y b = 0 .
3. No tener solución, si a = 0 y b  0 .
De lo comentado arriba, se observa que el estudio de las raíces o soluciones de una
ecuación lineal es un tema muy simple.
Una ecuación no lineal es aquella que tiene un formato diferente al de la lineal, por ejemplo
ax 2 + bx + c = 0 , sen( x 2 ) + x cos(x) + c = 0 ,…etc.

A diferencia de las ecuaciones lineales, las no-lineales presentan, en general, mucha más
dificultad en el análisis de sus raíces.
Proceso de cálculo de las raíces reales de ecuaciones no-lineales en una variable
f ( x) = 0

1. Acotación de las raíces

Se trata de encontrar un intervalo Ci, Cs  R , llamado de acotación de las raíces,
de modo que todas las raíces (si existen) se encuentren dentro de este intervalo.
Esto puede lograrse construyendo el grafico de f en un intervalo suficientemente
amplio.

2. Separación de las raíces

En esta etapa, deben encontrarse pequeños intervalos, cada uno conteniendo en su
interior una única raíz.
Para lograr este cometido, recurrimos a los siguientes teoremas del cálculo:
 a. Teorema de Bolzano
Sea g una función continua en el intervalo real a, b
Si g (a) g (b)  0 , entonces en a, b existe por lo menos una raíz de g (x) .

b. Teorema de Rolle
Sea g una función derivable en el intervalo real a, b .

, entonces en a, b existe c tal que

dg
Si (c ) = 0 .
dx

3. Cálculo aproximado de las raíces de la ecuación

Una vez separadas las raíces, el cálculo de las mismas es efectuado utilizando
métodos iterativos. Estos métodos son de dos tipos:
 1. Métodos de intervalo o cerrados
Mantienen a la raíz dentro de un intervalo, que se cierra en torno de ella, desde
el inicio de los cálculos hasta la obtención final de la raíz.
Son métodos seguros (no pierden a la raíz), pero son lentos, ya que hay muchas
operaciones involucradas en la actualización de los intervalos más pequeños.
2. Métodos abiertos o de punto(s) inicial(es)
No necesitan intervalo en torno a la raíz, solo requieren punto(s) cercano(s) a la
raíz. Son métodos inseguros (pueden perder a la raíz), pero son rápidos ya que
no se realizan operaciones de actualización de intervalo.
Ambas familias de métodos utilizan una formula iterativa, propia del método, que
funciona como una fábrica de aproximaciones a la raíz buscada.
La fórmula iterativa de cada método se obtiene cambiando la curva f por una recta, y
calculando la aproximación de la raíz, como la intersección de dicha recta con el eje x .

En la practica la secuencia de aproximaciones a la raíz, generada por un metodo, puede no

ser convergente y el método falla.
El procedimiento se detiene cundo el error, calculado entre dos iteraciones consecutivas
x k , x k +1 resulta menor que la tolerancia fijada en el problema.

Métodos de intervalo
1. Método de la bisección

ak + bk
Utiliza la siguiente formula iterativa: x k +1 =
2
 
, donde a k , b k es el intervalo

 
actualizado en la iteración k . El procedimiento inicia con un intervalo a 0 ,b 0 .
Calculada la primer aproximación x1 , se actualiza el intervalo que tendrá como
uno de sus extremos al nuevo valor x1 , de modo que a 0 o b 0 deberá cambiarse
por x1 . Esa decisión depende de los signos de f (a 0 ), f ( x1 ) y f (b 0 ) .
Si f (a 0 ) f ( x1 )  0 entonces b1 = x1
Si f (b 0 ) f ( x1 )  0 entonces a1 = x1 .
 
Cuando se trate de actualizar el intervalo a k , b k , usando la nueva
aproximación x k +1 , se tiene que realizar el siguiente análisis:
Si f (a k ) f ( x k +1 )  0 entonces b k +1 = x k +1
Si f (b k ) f ( x k +1 )  0 entonces a k +1 = x k +1 .
El número de iteraciones n que se tendrán que ejecutar en este método, para
alcanzar la raíz, depende de la tolerancia tol y la longitud del intervalo inicial,
valores fijados al inicio del problema.
Note que cada iteración reduce la longitud del intervalo a la mitad, luego de las
b0 − a0
n iteraciones, la longitud el intervalo inicial se reduce a y el
2n
b0 − a0
procedimiento concluye cuando  tol , entonces considerando la
2n
igualdad en la expresión anterior y usando logaritmos para despejar el número
ln(b 0 − a 0 ) − ln(tol )
total de iteraciones n , se obtiene n = .
ln( 2)
 Ejemplo dada la ecuación x 2 − 2 = 0 , calcule su raíz en 1,2 . Utilice el método
de la bisección con tol = 0.0001.

2. Método de la Falsa Posición

( )( )
Utiliza la recta secante a través de los puntos a k , f (a k ) , b k , f (b k ) del intervalo
k +1
actualizado en la iteración k para calcular la nueva iteración x .

El grafico de arriba muestra el cálculo de la primera iteración usando la recta

( )( )
secante que pasa por los puntos a 0 , f (a 0 ) , b 0 , f (b 0 ) .
  b0 − a0
x = a +  , donde
1 0
=
− f (a 0 ) f (b 0 ) − f (a 0 )
 b0 − a0 
 = − f (a 0 )  0 
 f (b ) − f (a ) 
0

 b0 − a0 
x1 = a 0 − f (a 0 )  0 
 f (b ) − f (a ) 
0

Para la iteración k + 1, resulta la formula iterativa general del método

k +1  bk − ak 
x = a − f (a ) 
k k
k 
 f (b ) − f (a ) 
k

Ejemplo dada la ecuación x 2 − 2 = 0 , calcule su raíz en 1,2 . Utilice el método de la Falsa

Posición con tol = 0.0001.

a x b Error
1 1.3333333 2
1.3333333 1.4 2
1.4 1.4117647 2
1.4117647 1.4137931 2
1.4137931 1.4141414 2
1.4141414 1.4142012 2 <tol

Modificación a la Falsa Posición

Para algunas ecuaciones, el método de la falsa posición mantiene fijo uno de los extremos
del intervalo actualizado, durante todas o un número grande de iteraciones. Este fenómeno
contribuye a disminuir la velocidad del método.
Para resolver este problema se modifica ligeramente el método, de la siguiente manera:
Sea e el extremo del intervalo actualizado que se mantiene fijo durante dos iteraciones
consecutivas, entonces donde aparece f (e) en la forrmula iterativa de la falsa posición,
f (e)
esta expresión se sustituye por y se aplica la iteración, si el mismo extremo continua
2
f (e)
fijo, entonces cambiamos f (e) por 2 y si continua fijo seguimos aplicando la
2
modificación aumentando la potencia de 2 hasta conseguir que el extremo fijo se mueva.
Una vez se moviliza el extremo e se aplica la formula iterativa del método sin
modificación. El procedimiento vuelve a repetirse si nuevamente uno de los dos extremos
se mantiene fijo durante dos iteraciones consecutivas.

Ejemplo dada la ecuación x 2 − 2 = 0 , calcule su raíz en 1,2 . Utilice el método de la Falsa

Posición Modificada con tol = 0.0001.
a x b Error
1(N) 1.33333333 2(N)
1.33333333333(N) 1.4 2(S)
1.4(N) 1.42307692 2(S)
1.4(S) 1.41416894 1.42307692307692(N)
1.4141689373297(N) 1.41421342 1.42307692307692(S) <tol

Métodos de punto(s) inicial(es)

1. Método de la secante

Utiliza dos puntos iniciales a, b y aproxima la raíz con la recta secante pasando por
los puntos (a, f (a)) y (b, f (b)) .
Los puntos iniciales a, b no necesitan formar intervalo en torno de la raíz, solo es
requerido que sean cercanos a la raíz.

El primer paso consiste en identificar al punto más cercano x 2 y al más alejado x1

de la raíz. Se utiliza el siguiente criterio:
𝑆𝑖|𝑓(𝑎)| ≥ |𝑓 (𝑏)𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠|→𝑥 1 = 𝑎 𝑦 𝑥 2 = 𝑏
𝑆𝑖|𝑓(𝑎)| ≤ |𝑓 (𝑏)𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠|→𝑥 1 = 𝑏 𝑦 𝑥 2 = 𝑎
La fórmula iterativa se deduce del grafico mostrado:

 x 2 − x1 
x 3 = x 2 +  , pero  = − f ( x 2 )  1 
Luego
 f (x ) − f (x ) 
2

 x 2 − x1 
x = x − f ( x )
3 2 2
1 
 f (x ) − f (x ) 
2
  x k − x k −1 
Generalizando se obtiene x k +1 = x k − f ( x k )  k −1 
 f (x ) − f (x ) 
k

Ejemplo dada la ecuación x 2 − 2 = 0 , calcule su raíz en 1,2 . Utilice el método de la

Secante con tol = 0.0001.
2 1 1.33333333
1 1.33333333 1.42857143
1.33333333 1.42857143 1.4137931
1.42857143 1.4137931 1.41421144
1.4137931 1.41421144 1.41421356
1.41463415 1.41421144 1.41421356 <tol

2. Metodo de Newton-Raphson

Utiliza un punto inicial x 0 , el cual debe estar cerca de la raíz, y aproxima la raíz con
la recta tangente al gráfico de f por el punto ( x 0 , f ( x 0 ) .
 − f (x0 ) f (x0 )
x = x +,  =
1 0
por lo tanto x = x −
1 0
.
f ( x 0 ) f ( x 0 )
f (x k )
Generalizando x k +1 = x k −
f ( x k )

Cuando tenemos un intervalo conteniendo a la raíz, como dato inicial, entonces podemos
elegir entre los dos puntos del intervalo, el punto e que cumple f (e) f (e)  0 .

Ejemplo dada la ecuación x 2 − 2 = 0 , calcule su raíz en 1,2 . Utilice el método de Newton-

Raphson con tol = 0.0001.
x(k-1) x(k)
1 1.5
1.5 1.41666667
1.41666667 1.41421569
1.41421569 1.41421356
Problema del Cable

Ecuaciones de equilibrio

F x = 0 → H1 = H 2 = H

M Apoyo2
z
= 0 → V1 Lx + HLy + P1 ( Lx − x1 ) + P2 ( Lx − x2 ) = 0

V1 = −
1
Lx

HL y + P1 ( Lx − x1 ) + P2 ( Lx − x 2 ) 

F y = 0 → V1 + V2 + P1 + P2 = 0

V2 = −V1 + P1 + P2 

Falta una ecuación para determinar las reacciones.

Ahora se utiliza la extrema flexibilidad del cable, que permite asumir momento flector nulo
en todos sus puntos.

M P1
z , izq
= 0 → V1 x1 + Hy1 = 0

V1
y1 = − x1
H

M P2
z , izq
= 0 → V1 x2 + Hy2 + P1 ( x2 − x1 ) = 0

y2 = −
1
V1 x2 + P1 ( x2 − x1 )
H
Se incorporaron dos ecuaciones pero también dos incógnitas, continua faltando una
ecuación para salvar la indeterminación.
Se trata de una condición de compatibilidad L = L0 siendo L la longitud del cableen la
configuración de equilibrio.

L = x12 + y12 + ( x2 − x1 ) 2 + ( y 2 − y1 ) 2 + ( Lx − x2 ) 2 + ( L y − y 2 ) 2

Finalmente, La ecuación no-lineal que debe resolverse es: f ( H ) = 0 siendo f ( H ) = L − L0 .

Ejemplo del cable

Datos : L0 = 12 m , Lx = 4 m , L y = 3 m , x1 = 2 m , x 2 = 3 m , P1 = 10 tnf , P2 = 5 m .

Calcule H usando bisección.

Primero se determina el intervalo que contiene a la raíz.

Del grafico se deduce que la raíz de la ecuación esta entre 2 tnf. y 4tnf,
A continuación, utilizamos bisección para el cálculo de la raíz, con tol=0.0001.
Hi H Hf Error
2 3 4
2 2.5 3
2 2.25 2.5
2 2.125 2.25
2.125 2.1875 2.25
2.1875 2.21875 2.25
2.21875 2.234375 2.25
2.234375 2.2421875 2.25
2.2421875 2.24609375 2.25
2.24609375 2.248046875 2.25
2.24609375 2.247070313 2.248046875
2.24609375 2.246582031 2.247070313
2.24609375 2.246337891 2.246582031
2.246337891 2.246459961 2.246582031
2.246459961 2.246520996 2.246582031 <tol

H=2.246520996 tnf.