1
CAPÍTULO III
ECUACIONES LINEALES Y NO LINEALES
CONTENIDO
1.- Definición
2.- Método de Horner y Horner Iterado
3.- Método de Bisección
4.- Método de Newton Raphson
5.- Método de Punto Fijo
6.- Método de la Secante
7.- Ejercicios
--------------- * ---------------
1.- Definición.-
Una ecuación lineal o de primer grado, es aquella que tiene la forma de un polinomio y que
involucra solamente sumas y restas de variables elevadas a la primera potencia (elevadas a uno,
que no se escribe); no se consideran variables multiplicadas entre sí , ni en el denominador.
Pueden ser representadas como rectas en el sistema cartesiano.
Mientras que una ecuación no lineal, serán aquellas relaciones de igualdad que consideran y o más
variables con exponentes diferentes a 1.
2.- Método de Hornen y Horner Iterado.-
El método de Horner, también llamado la regla de Horner es un algoritmo que permite calcular el
resultado de un polinomio para un determinado valor de x; es para evaluar de forma eficiente las
funciones polinómicas de una forma monomial.
El procedimiento a seguir para resolver una ecuación mediante el método de Horner cuenta con
los siguientes pasos:
a) El dividendo y el divisor deben ser polinomios ordenados de manera decreciente y completa
respecto de la misma variable.
b) En caso que el polinomio no sea completo, se completará con ceros los términos que faltan en
el dividendo y divisor.
2
c) Se colocan los coeficientes del dividendo con su signo; se colocan los coeficientes del divisor
todos cambiados de signo menos el primero que se mantiene.
d) Colocar la línea vertical que separa el cociente del residuo. Esta línea se obtiene contando de
derecha a izquierda tantas columnas como nos indica el grado del divisor.
e) El primer coeficiente del cociente se obtiene dividiendo el primer coeficiente del dividendo
entre el primero del divisor. Luego este (primer coeficiente del cociente) se multiplica por cada
uno de los coeficientes del divisor que han cambiado de signo y el resultado se coloca en la
segunda fila, corriéndose un lugar a la derecha.
f) Se reduce la siguiente columna (se suman los coeficientes), y se repite el paso anterior tantas
veces hasta que la última operación efectuada caiga debajo del último coeficiente del dividendo.
g) Se suman directamente los números que están en las columnas que corresponde a los
coeficientes del residuo.
h) El grado del polinomio del cociente queda determinado por la diferencia entre los grados del
dividendo y divisor; y el grado del residuo queda determinado según la cantidad de términos.
con su mismo
signo D D I V I D E N D O
I
V
Con signo I
cambiado S
O
R COCIENTE RESIDUO
Horner Iterado
Método basado en la división sintética repetitiva para la obtención de las soluciones a partir de
una posible solución que puede ser iniciada en los rangos que presenta cambio de signo en la
función evaluada para cada valor asignado. Es decir se toma el polinomio de grado n y se divide
entre (x - t), donde t es un valor particular de x.
La primera división sintética da como resultado un cociente y un residuo R 1(x); la segunda división
sintética se realiza entre el cociente obtenido de la primera división y el mismo valor del divisor,
obteniéndose un nuevo cociente y residuo R 2(x). Con estos residuos y el valor inicial de la posible
raíz se calcula una nueva aproximación mediante:
3
xi+1 = xi – [R1(x) / R2(x)]
Con la nueva aproximación xi+1 se repiten las divisiones sintéticas para la obtención de los nuevos
residuos y con estos la nueva aproximación, proceso que se repite hasta que las aproximaciones
penúltima y última sean aproximadamente iguales.
Ventajas:
- Se puede expresar un polinomio en forma: P(x) = (x – x 0) Q(x) + dn .
- Se trabaja fácil y rápido.
Desventajas:
- Sólo es útil para funciones polinómicas no lineales.
Ejemplo
Hallar una aproximación a una de las raíces del polinomio 2x 4 – 3x2 + 3x – 4 utilizando como valor
inicial de partida x = -2.
1º Iteración
2 0 -3 3 -4
-2 -4 8 -10 14
2 -4 5 -7 10
-2 -4 16 -42
2 -8 21 -49
La nueva aproximación según xi+1 = xi – [R1(x) / R2(x)] es xi+1 = (-2) – [10 / (-49)] = -1,796
Ahora se repite el proceso pero con esta nueva aproximación:
2º Iteración
2 0 -3 3 -4
-1.746 -3.592 6.451 -6.198 5.744
2 -3.592 3.451 -3.198 1.744
-1.746 -3.592 12.90 -29.36
2 -7.184 16.35 -32.56
xi+1 = (-1.796) – [1.744 / (-32.56)] = -1,742
4
3º Iteración
2 0 -3 3 -4
-1.742 -3.484 6.069 -5.346 4.086
2 -3.484 3.069 -2.346 0.086
-1.742 -3.484 12.138 -26.490
2 -6.968 15.207 -28.836
xi+1 = (-1.742) – [0.086 / (-28.836)] = -1,739
4º Iteración
2 0 -3 3 -4
-1.739 -3.478 6.048 -5.300 3.999
2 -3.478 3.048 -2.300 -0.001
-1.739 -3.478 12.096 -26.335
2 -6.956 15.144 -28.635
xi+1 = (-1.739) – [-0.001 / (-28.635)] = -1,739
La raíz aproximada es – 1,739.
3.- Método de Bisección
El método de Bisección para la resolución de la ecuación f(x) = 0 se basa en el Teorema de Bolzano
que nos asegura la existencia de, al menos, una raíz de una función f(x) en un cierto intervalo [a,
b], bajo ciertas condiciones.
Este método busca la solución de una ecuación en el intervalo [a, b] suponiendo que el signo de
f(a) es distinto que el signo de f(b), es decir, f(a)*f(b) < 0. En tal caso, es evidente que la función
cortará al eje X en un punto o en un número impar de puntos, por lo que seguramente habrá por
lo menos una solución.
Para aplicar el algoritmo, este consiste en obtener el punto medio del intervalo x m = [(a+b)/2] y
hacer que uno de los extremos del intervalo sea este nuevo punto calculado: si f(a)*f(m) < 0,
entonces se hace que b = x m mientras a = a (se mantiene), y en caso contrario, a = x m mientras b = b
(se mantiene); de tal manera que siempre los dos extremos del intervalo cumplan que el signo de
f(a) sea distinto que el de f(b). El procedimiento se repite obteniendo sucesivamente los nuevos
puntos medios y parará si el número de iteraciones es mayor que un número prefijado o si el valor
absoluto de la función en el punto medio calculado es menor que cierta cifra o es casi un cero.
5
Ventajas:
- Es siempre convergente.
- Es óptimo para ecuaciones f(x)=0 cuando no se sabe nada de f, excepto calcular su signo.
- Es fácil de aplicar.
Desventajas:
- Converge muy lentamente.
- Sólo encuentra una raíz, aunque existan más en el intérvalo.
- Difícil generalizar para dimensiones superiores.
Ejemplo
Resolver la ecuación f(x) = ex – x4 + 1 = 0 en el intervalo [1; 2]
# ai bi xm f(ai) f(xm) f(ai)*f(xm)
1 1 2 1.5 2.718 0.419 1.138
2 1.5 2 1.75 0.419 -2.624 -1.099
3 1.5 1.75 1.625 0.419 -0.894 -0.374
4 1.5 1.625 1.562 0.419 -0.189 -0.079
5 1.5 1.562 1.531 0.419 0.126 0.052
6 1.531 1.562 1.546 0.126 -0.028 -0.003
7 1.531 1.546 1.539 0.126 0.049 0.006
8 1.539 1.546 1.542 0.049 0.010 0.0004
9 1.542 1.546 1.544 0.010 -0.009 -0.00009
10 1.542 1.544 1.543 0.010 0.0006 0.000006
La raíz es aproximadamente 1,543.
4.- Método de Newton Raphson
Este método es uno de los más utilizados para localizar raíces ya que en general es muy eficiente y
siempre converge para una función polinomial. Se requiere que las funciones sean diferenciables,
y por tanto, continuas, para poder aplicar este método.
Se debe partir de un valor inicial para la raíz: xi, este puede ser cualquier valor, el método
convergirá a la raíz más cercana. Si se extiende una tangente desde el punto b, el punto donde
esta tangente cruza al eje x representa una aproximación mejorada de la raíz.
El esquema iterativo de Newton puede derivarse del desarrollo de Taylor de la función alrededor
de la estimación inicial.
6
La fórmula de Newton-Raphson se deduce a partir de la fórmula de la pendiente de una recta.
Pendiente de una recta:
El valor absoluto de la diferencia de la m debe ser menor que la tolerancia o el resultado de alguna
fórmula de error debe ser menor que la tolerancia dada. Una de las fórmulas de error más útiles es
la del error relativo porcentual aproximado.
El método de Newton-Raphson es convergente cuadráticamente, es decir, el error es
aproximadamente al cuadrado del error anterior.
7
Esto significa que el número de cifras decimales correctas se duplica aproximadamente en cada
interacción. Cuando el método de Newton-Raphson converge, se obtienen resultados en
relativamente pocas interacciones, ya que para raíces no repetidas este método converge con
orden 2 y el error Ei+1 es proporcional al cuadrado del resultado anterior E i.
Ventajas:
- Es el más conocido.
- Es eficiente en la solución de ecuaciones no lineales.
- Converge rápidamente.
- Proporciona una buena precisión en los resultados
Desventajas:
- No existe un criterio general de convergencia.
- Las funciones deben ser derivables.
- Evaluación de las derivadas.
Ejemplo
Hallar una raíz de la función f(x) = x 3 + 4x2 -10 a partir de x0 = 1
Para lo cual se requiere de la derivada de la función, que está dada por: f’(x) = 3x 2 + 8x
i xi f(xi) f’(xi) xi+1
0 1 -5 11 1.454
1 1.454 1.540 17.983 1.368
2 1.368 0.060 16.572 1.365
3 1.365 0.0001 16.513 1.365
La solución es aproximadamente 1.365
5.- Método de Punto Fijo
Éste método se aplica para resolver ecuaciones de la forma x = g(x). Si la ecuación es f(x) = 0,
entonces puede despejarse x ó bien sumar x en ambos lados de la ecuación para ponerla en la
forma adecuada.
El procedimiento empieza con una estimación o conjetura inicial de x, que es mejorada por
iteración hasta alcanzar la convergencia. Para que converja, la derivada (dg/dx) debe ser menor
que 1 en magnitud (al menos para los valores x que se encuentran durante las iteraciones). La
convergencia será establecida mediante el requisito de que el cambio en x de una iteración a la
siguiente no sea mayor en magnitud que alguna pequeña cantidad ε.
El algoritmo a seguir para la iteración de punto fijo indica primero la ubicación de una raíz de f(x)
analizando la gráfica. Despejar x = g(x), luego obtener de x = g(x) su derivada g’(x). Calcular g’(x) de
cada caso para el valor inicial x0 para encontrar el despeje óptimo que debe cumplir -1 ≤ g’(x) ≤ 1.
8
Se inician las iteraciones con el valor inicial x0 en la función g(x) óptima hasta igual valores.
Ventajas:
- Es el método que converge con rapidez
- Cuando converge es de alta precisión.
- No necesita de un intervalo para funcionar sino únicamente un punto perteneciente al
intervalo donde esté la raíz.
Desventajas:
- No garantiza la convergencia.
- La función g(x) correcta puede ser muy compleja de encontrar.
- Hay infinidad de g(x) y no existe regla para escoger la correcta.
Ejemplo
Se tiene la siguiente ecuación x 3– x2 – x - 1= 0, la cual mediante despejes diversos podemos dejar
como siguen y evaluar cada derivada para x 0 = 0 ó para el intervalo que presenta cambio de signo
en el teorema del resto.
g1(x) = x3– x2 – 1 g’1(x)= 3x2 - 2x = 0
g2(x) = (x3 – x – 1)(1/2) g’2(x)= (1/2) (x3 – x – 1)(-1/2)(3x2 - 1)= imaginario
g3(x) = (x2 + x + 1)(1/3) g’3(x)= (1/3) (x2 + x + 1)(-2/3)(2x + 1)= (1/3)
g4(x) = (x3 – 1)/(x+1) g’4(x)= (2x3 + 3x2 +1)/(x+1)2= 1
tanto g1, g3 y g4 están entre el rango [-1, 1] y son considerados despejes óptimos. Elegiremos a
g3(x) por estar más cerca del cero que indica una convergencia más rápida.
i xi g(xi) |xi+1 – xi|
0 0 1
1 1 1.442 1
2 1.442 1.653 0.442
3 1.752 1.798 0.31
4 1.798 1.820 0.046
5 1.830 1.835 0.032
6 1.835 1.837 0.005
7 1.837 1.838 0.002
8 1.838 1.838 0.001
9 1.838 1.838 0
La solución es aproximadamente 1.838.
9
6.- Método de la Secante
El método de la secante es un método para encontrar los ceros de una función de
forma iterativa. Es una modificación de método convencional de Newton Raphson, donde en vez
de calcular la derivada de la función en el punto de estudio, teniendo en mente la definición de
derivada, se aproxima la pendiente a la recta que une la función evaluada en el punto de estudio y
en el punto de la iteración anterior.
El método de la secante parte de dos puntos (y no sólo uno como el método de Newton) y estima
la tangente (es decir, la pendiente de la recta) por una aproximación de acuerdo con la expresión:
S
u
s
t
i
t
u
y
e
n
d
o
e
s
t
a
e
x
p
r
e
s
i
ó
n
10
e
n
l
a
e
c
u
a
c
i
ó
n
d
e
l
m
é
t
o
d
o
d
e
N
e
w
t
o
n
,
o
b
t
e
n
11
e
m
o
s
l
a
e
x
p
r
e
s
i
ó
n
d
e
l
m
é
t
o
d
o
d
e
l
a
s
e
c
a
n
t
e
12
q
u
e
n
o
s
p
r
o
p
o
r
c
i
o
n
a
e
l
s
i
g
u
i
e
n
t
e
p
u
n
t
o
d
e
13
i
t
e
r
a
c
i
ó
n
:
C
u
y
a
r
e
p
r
e
s
e
n
t
a
c
i
ó
n
g
e
o
m
é
14
t
r
i
c
a
e
s
l
a
s
i
g
u
i
e
n
t
e
:
D
e
b
i
d
o
a
15
q
u
e
n
o
s
e
p
u
e
d
e
a
s
e
g
u
r
a
r
l
a
c
o
n
v
e
r
g
e
n
c
i
a
16
t
a
l
c
o
m
o
e
n
e
l
m
é
t
o
d
o
d
e
N
e
w
t
o
n
-
R
a
p
h
s
o
n
,
h
17
a
y
q
u
e
i
n
t
e
n
t
a
r
c
o
n
d
o
s
p
u
n
t
o
s
d
e
p
a
r
t
i
d
a
18
a
n
a
l
i
z
a
r
i
t
e
r
a
c
i
ó
n
i
t
e
r
a
c
i
ó
n
e
l
c
o
m
p
o
r
19
t
a
m
i
e
n
t
o
d
e
l
o
s
r
e
s
u
l
t
a
d
o
s
.
T
a
m
b
i
é
n
p
u
e
d
e
20
t
e
n
e
r
s
e
e
n
c
u
e
n
t
a
c
i
e
r
t
o
s
c
r
i
t
e
r
i
o
s
q
u
e
u
t
21
i
l
i
z
a
e
l
m
é
t
o
d
o
p
a
r
a
o
b
t
e
n
e
r
l
a
s
i
g
u
i
e
n
t
e
22
a
p
r
o
x
i
m
a
c
i
ó
n
,
p
o
r
e
j
e
m
p
l
o
,
l
a
p
e
n
d
i
e
n
t
e
d
e
23
l
a
r
e
c
t
a
s
e
c
a
n
t
e
n
o
d
e
b
e
s
e
r
n
u
l
a
p
u
e
s
n
u
24
n
c
a
i
n
t
e
r
s
e
c
t
a
r
á
c
o
n
e
l
e
j
e
n
i
t
e
n
d
i
e
n
d
o
25
i
n
f
i
n
i
t
o
(
s
i
e
m
p
r
e
c
u
a
n
d
o
l
a
r
a
í
z
s
e
e
26
n
c
u
e
n
t
r
e
a
l
e
j
a
d
a
l
o
s
p
u
n
t
o
s
i
n
i
c
i
a
l
e
s
)
y
27
q
u
e
l
a
c
o
n
v
e
r
g
e
n
c
i
a
p
u
e
d
e
l
l
e
g
a
r
s
e
r
l
28
e
n
t
a
.
D
e
t
o
d
a
s
m
a
n
e
r
a
s
,
d
e
p
e
n
d
e
r
á
d
e
l
a
f
o
29
r
m
a
q
u
e
t
i
e
n
e
l
a
f
u
n
c
i
ó
n
.
E
n
g
e
n
e
r
a
l
,
e
l
m
30
é
t
o
d
o
d
e
l
a
s
e
c
a
n
t
e
p
r
e
s
e
n
t
a
l
a
s
m
i
s
m
a
s
v
e
31
n
t
a
j
a
s
l
i
m
i
t
a
c
i
o
n
e
s
q
u
e
e
l
m
é
t
o
d
o
d
e
N
e
w
32
t
o
n
-
R
a
p
h
s
o
n
e
x
p
l
i
c
a
d
o
a
n
t
e
r
i
o
r
m
e
n
t
e
.
Ejemplo
U
t
i
33
l
i
c
e
e
l
m
é
t
o
d
o
d
e
l
a
s
e
c
a
n
t
e
p
a
r
a
e
n
c
o
n
t
r
a
34
u
n
a
r
a
í
z
r
e
a
l
d
e
l
a
e
c
u
a
c
i
ó
n
p
o
l
i
n
o
m
i
a
l
:
35
La solución es aproximadamente 0.652
