CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

 
 
 
CAPÍTULO 3 
SOLUCIÓN  DE  ECUACIONES  NO 
LINEALES 
 
 
 
3.1 Introducción 
Desde muchos años hemos aprendido a usar la fórmula cuadrática: 

− b ± b 2 − 4ac
x=      (3.1) 
2a
Para resolver 
f ( x) = ax 2 + bx + c = 0     (3.2) 
A los valores calculados con la relación (3.1), se le llama las “raíces” de la ecuación 
(3.2), que representa los valores de x que verifican la ecuación (3.2) igual a cero. 
Debido a esto, algunas veces a las raíces se les conoce como ceros de la ecuación. 
Aunque  la  formula  cuadrática  es  útil  para  resolver  la  ecuación  (3.2),  existen 
muchas  funciones  en  la  ingeniería  donde  las  raíces  no  se  pueden  determinar  tan 
fácilmente. Por ejemplo, incluso una función tan simple como  f ( x) = e − x − x  no se 
puede resolver en forma analítica. En estos casos, los métodos numéricos descritos 
a continuación, proporcionaran medios eficientes para obtener la respuesta. 
 

67
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

3.2 Métodos cerrados 
Estos  métodos  aprovechan  el  hecho  de  que  una  función  cambia  de  signo  en  la 
vecindad de una raíz. A esta técnica se les llama métodos cerrados, o de intervalos, 
por que se necesita de dos valores iniciales para la raíz. Como su nombre lo indica, 
dichos  valores  iniciales  deben  “encerrar”,  o  estar  a  ambos  lados  de  la  raíz.  Los 
métodos  particulares  descritos  a  continuación,  emplean  diferentes  estratégicas 
para  reducir  sistemáticamente  el  tamaño  del  intervalo  y  así  converger  a  la 
respuesta correcta. 
Como  preámbulo  de  estas  técnicas  se  utilizarán  los  métodos  gráficos  para 
representar  tanto  las  funciones  como  sus  raíces.  Además  los  métodos  gráficos 
permiten  determinar  valores  iniciales,  también  son  útiles  para  determinar  las 
propiedades  de  las  funciones  y  el  comportamiento  de  los  diversos  métodos 
numéricos. 
3.2.1 Método gráfico 
Un método simple para obtener una aproximación a la raíz de la  ecuación f(x)=0 
consiste  en  graficar  la  función  y  observar  dónde  cruza  al  eje  x.  Este  punto,  que 
representa el valor de x para el cual f(x)=0, ofrece una aproximación inicial de la 
raíz. 

Ejemplo 3.1: La función  f ( x ) = sen (10 x ) + cos (3 x ) , tiene varias raíces en 

el  rango  [0,  5],  utilice  grafica  por  computadora  para  comprender  mejor  el 
comportamiento de esta función. 
>> fplot('sin(10*x)+cos(3*x)',[0,5]),grid 
2

1.5

1
Y =sin(10*x)+cos(3*x)

0.5

-0.5

-1

-1.5

-2
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
X

68
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Se puede apreciar que tiene 13 raíces en este intervalo, a excepción del intervalo 
[4, 4.5] que puede ser una o dos raíces adicionales. 
>> fplot('sin(10*x)+cos(3*x)',[4,4.5]),grid 
1.6

1.4

1.2

1
Y =sin(10*x)+cos(3*x)

0.8

0.6

0.4

0.2

-0.2
4 4.05 4.1 4.15 4.2 4.25 4.3 4.35 4.4 4.45 4.5
X

Ampliando el grafico podemos observar que existen dos raíces entre 4.2 y 4.3 
 
>> fplot('sin(10*x)+cos(3*x)',[4.2,4.3]),grid 

0.12

0.1

0.08
Y =sin(10*x)+cos(3*x)

0.06

0.04

0.02

-0.02
4.2 4.21 4.22 4.23 4.24 4.25 4.26 4.27 4.28 4.29 4.3
X
 
Finalmente las dos raíces diferentes se encuentran entre x=4.22 y x=4.27. 
 

69
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Ejemplo 3.2: Localizar las raíces de: f ( x ) = x + 2x − 4  
2

>> fplot('[x.^2+2*x‐4]',[‐4,4]),grid 

20

15

10
y=x.2+2*x-4

-5
-4 -3 -2 -1 0 1 2 3 4
x
 
x = 4 − 2x  
2

Si separamos la función principal en dos sub funciones y graficamos, obtendremos 
también las raíces en la intersección de las nuevas funciones. 

   f1 ( x) = x 2   f 2 ( x) = 4 − 2 x  
>> fplot('[x.^2,4‐2*x]',[‐4,4]),grid 

16

14

12

10

8
x.2,4-2*x

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x
 

70
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Ejemplo 3.3: Analizar gráficamente la función:  f ( x) = e − x − x  
>> fplot('exp(‐x)‐x',[‐2,2]),grid 
y=exp(-x)-x
10

6
Eje y

-2
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
Eje x
 
Como en el ejemplo anterior, separamos la función principal en dos sub funciones 
y graficamos, obtenemos las raíces en la intersección de las nuevas funciones. 
−x
e ;   =x
f1 ( x) = e − x   2   f ( x) = x
>> fplot('exp(‐x)',[‐2,2],'r') 
>> hold on 
>> fplot('x',[‐2,2],'b'),grid 
y=exp(-x); y=x
8

4
Eje y

-1

-2
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
Eje x
 

71
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

3.2.2 Método de Bisección 
Al aplicar las técnicas gráficas, si se observa que una función cambia de signo en la 
vecindad de una raíz. En general, si f(x) es real y continua en el intervalo que va 
desde xl hasta xu y f(xl) y f(xu) tienen signos opuestos, es decir: 

f ( xl ) f ( xu ) < 0  
Entonces hay al menos una raíz real entre xl y xu 
 
Los métodos de búsqueda incremental aprovechan esta característica localizando 
un  intervalo  en  el  que  la  función  cambia  de  signo.  Entonces,  la  localización  del 
cambio  de  signo  (y,  en  consecuencia,  de  la  raíz)  se  logra  con  más  exactitud  al 
dividir  el  intervalo  en  varios  subintervalos.  Se  investiga  cada  uno  de  estos 
intervalos  para  encontrar  el  cambio  de  signo.  El  proceso  se  repite  se  repite  y  la 
aproximación a la raíz mejora cada vez más en la medida que los subintervalos se 
dividen en intervalos cada vez más pequeños. 
 
El método de bisección, conocido también como de corte binario, de partición de 
intervalos o de Bolzano, es un tipo de búsqueda incremental en el que el intervalo 
se divide siempre a la mitad. Si la función cambia de signo sobre un intervalo, se 
evalúa el valor de la función en el punto medio. La posición de la raíz se determina 
situándola en el punto medio del subintervalos, dentro del cual ocurre un cambio 
de signo. El proceso se repite hasta obtener una mejor aproximación. 
Una representación gráfica del método de bisección. 
 
f(x)

f(xu)

xl xr x
0
xu
f(xr)
f(xl)
 
 

72
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Procedimiento del método de bisección: 
Paso 1: Elija valores inicial inferior, xl, y superior, xu, que encierren la raíz, de tal 
forma que la función cambie de signo en el intervalo. Esto se verifica comprobando 
que  f ( xl ) ⋅ f ( xu ) < 0 . 

Paso 2: Una aproximación de la raíz xr se determina mediante: 

x l + xu
xr =  
2
Paso  3:  Realice  las  siguientes  evaluaciones  para  determinar  en  qué  subintervalo 
está la raíz: 
a) Si  f ( xl ) ⋅ f ( x r ) < 0 ,  entonces  la  raíz  se  encuentra  dentro  del  subintervalo 

inferior o izquierdo. Por lo tanto, haga  xu = x r  y vuelva al paso 2. 

b) Si  f ( xl ) ⋅ f ( x r ) > 0 ,  entonces  la  raíz  se  encuentra  dentro  del  subintervalo 

superior o derecho. Por lo tanto, haga  x l = x r  y vuelva al paso 2 

c) Si  f ( xl ) ⋅ f ( x r ) = 0 , la raíz es igual a xr; termina el cálculo. 

Se puede calcular el error relativo porcentual  ε a de la siguiente manera: 

− xr
nuevo anterior
xr
εa = nuevo
× 100  
xr

Donde  xrnuevo  es  la  raíz  en  la  iteración  actual  y  xranterior  es  el  valor  de  la  raíz  en  la 
iteración anterior. 
Ejemplo 3.4: 
Utilice  el  método  de  Bisección  para  calcular  la  raíz  de  la  función 
f ( x) = e − x − x ,  empleando  como  valores  de  intervalo  xl=0,  xu=1.  Con  una 
aproximación de 4 cifras significativas. 
 
Solución: 
Empleando el procedimiento citado del método de bisección se tiene: 

Iter.  xl  xu  f(xl)  f(xu)  f(xl)*f(xu)  xr  f(xr)  f(xl)*f(xr) εa  

1  0  1  1.0000  ‐0.6321  ‐0.6321  0.5000  0.1065  0.10653   
2  0.5000  1  0.1065  ‐0.6321  ‐0.0673  0.7500  ‐0.2776  ‐0.02958 33.33 
3  0.5000  0.7500  0.1065  ‐0.2776  ‐0.0296  0.6250  ‐0.0897  ‐0.00956 20.00 
4  0.5000  0.6250  0.1065  ‐0.0897  ‐0.0096  0.5625  0.0073  0.00078  11.11 
5  0.5625  0.6250  0.0073  ‐0.0897  ‐0.0007  0.5938  ‐0.0415  ‐0.00030 5.26 
6  0.5625  0.5938  0.0073  ‐0.0415  ‐0.0003  0.5781  ‐0.0172  ‐0.00013 2.70 
7  0.5625  0.5781  0.0073  ‐0.0172  ‐0.0001  0.5703  ‐0.0050  ‐0.00004 1.37 

73
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

8  0.5625  0.5703  0.0073  ‐0.0050 

‐0.000036150  0.5664  0.0012  0.00001  0.69 
9  0.5664  0.5703  0.0012  ‐0.0050 
‐0.000005734  0.5684  ‐0.0019  0.00000  0.34 
10  0.5664  0.5684  0.0012  ‐0.0019 
‐0.000002201  0.5674  ‐0.0004  0.00000  0.17 
11  0.5664  0.5674  0.0012  ‐0.0004 
‐0.000000434  0.5669  0.0004  0.00000  0.09 
12  0.5669  0.5674  0.0004  ‐0.0004 
‐0.000000146  0.5671387  0.000007  0.00000  0.043 
13  0.56714  0.5674  0.0000  ‐0.0004 
‐0.000000003  0.56726  ‐0.000184  0.00000  0.022 
14  0.56714  0.5673  0.0000  ‐0.0002 
‐0.000000001  0.567200  ‐0.000088  0.00000  0.011 
15  0.56714  0.5672  0.0000  ‐0.0001  ‐0.000000001  0.567169  ‐0.000041  0.00000  0.005 

La relación de la tolerancia porcentual prefijada, para 4 cifras significativas será: 
es = (0.5 × 10 2− 4 )% = 0.005 %  
Por lo que los cálculos se debe efectuar hasta que  ε a ≤ es  

 
Ejemplo 3.5: 
Utilizando el Método de Bisección, desarrolle y pruebe un programa en MATLAB, 

permita obtener la raíz de la función  f ( x) = e − x − x , así mismo el programa 
debe  permitir  ingresar  diferentes  valores  del  intervalo  xl,  xu,  y  número  de  cifras 
significativas,  para  así  permitir  aproximar  la  raíz  con  diferente  precisión. 
Finalmente probar el programa con los intervalo xl=0, xu=1 y con aproximaciones 
de 2, 3, 4, 5…15 cifras significativas. 
Solución: 

Se  deberá  crear  la  función  f ( x) = e − x − x   en  Matlab,  y  se  debe  guardar  el 
archivo con el mismo nombre de la función. 
Para acceder a la ventana de codificación en Matlab a través de la barra de menú 
es.  File/New/M­file,  al  cual  aparecerá  una  hoja  en  la  que  se  debe  codificar  la 
función como sigue: 
function y=funcion(x) 
y=exp(‐x)‐x; 

Seguidamente  se  debe  guardar  dicho  archivo  con  el  nombre:  funcion.m,  para 
posteriormente  poder  hacer  referencia  y  evaluarlo  desde  otro  punto  u  otro 
archivo. 

74
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Seguidamente se debe codificar el algoritmo del método de bisección en un nuevo 
archivo m, y luego guardar con un nombre referente al método de bisección, para 
nuestro ejemplo será: biseccion.m 
Algoritmo del Método de Bisección 
Leer: xl, xu, n 
Es=(0.5*10^(2‐n) 
Ea=100 
xr=0 
i=0 
Mientras Ea>Es hacer 
xa=xr 
xr=(xl+xu)/2) 
f(xl)*f(xr)<0 

xu=xr  xl=xr 
Ea=abs((xr‐xa)/xr)*100  Ea=abs((xr‐xa)/xr)*100 

i=i+1 

Imprimir: xr, Ea, i 
Código del Método Bisección (biseccion.m) 
clear;clc; 
disp('CALCULO DE LA RAIZ DE UNA ECUACION NO LINEAL') 
disp('POR EL METODO CERRADO DE BISECCION') 
n=input('Ingrese precision de cifras significativas n = '); 
xl=input('Ingrese limite inferior xl = '); 
xu=input('Ingrese limite superior xu = '); 
Es=(0.5*10^(2‐n)); 
Ea=100; 
xr=0; 
i=0; 
while Ea>Es 
    xa=xr; 
    xr=(xl+xu)/2; %FORMULA DE BISECCION 
    if funcion(xl)*funcion(xr)<0 
        xu=xr; 
        Ea=abs((xr‐xa)/xr)*100; 
    elseif funcion(xl)*funcion(xr)>0 
        xl=xr;   
        Ea=abs((xr‐xa)/xr)*100; 
    end 
    i=i+1; 
end 
fprintf('Raiz Solucion : %12.15f\n',xr) 
fprintf('Error aproximado : %12.10f\n',Ea) 
fprintf('Numero iteraciones : %12.0f\n',i) 

75
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Prueba del Programa

>> biseccion 
CALCULO DE LA RAIZ DE UNA ECUACION NO LINEAL 
POR EL METODO CERRADO DE BISECCION 
Ingrese precision de cifras significativas n = 4 
Ingrese limite inferior xl = 0 
Ingrese limite superior xu = 1 
Raiz Solucion:      0.567153930664063 
Error aproximado:      0.0026904141 
Numero iteraciones:           16 

3.2.3 Método de Falsa Posición 
Consiste  en  unir  f ( xl ) y f ( xu ) con  una  línea  recta.  La  intersección  de  esta  línea 

con  el  eje  de  las  x  representa  una  mejor  aproximación  de  la  raíz.  Lo  que  da  una 
falsa  posición  de  la  raíz,  del  cual  proviene  el  nombre  del  método  de  la  falsa 
posición,  o  del  latín,  regula  falsi.  También  se  le  conoce  como  método  de 
interpolación lineal. 

f(x)

f(xu)

xl xr
0 xu x

f(xl)
 
Representación gráfica del método de la falsa posición, con los triángulos semejantes sombreados 
se obtiene la fórmula para el método. 
 
Usando triángulos semejantes, la interpolación de la línea recta con el eje de las x 
se estima mediante la siguiente relación: 
f ( xl ) f ( xu )
Semejanza de triángulos     =  
x r − xl x r − xu
De la cual se despeja xr 

76
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

f ( xu )( xl − xu )
x r = xu −  
f ( x l ) − f ( xu )
Esta es la fórmula de la falsa posición. El valor de xr calculado con esta ecuación, 
reemplazará después a cualquiera de los valores iniciales, xl o xu, y dando un valor 
de la función con el mismo signo de f(xr). Con el mismo criterio que el método de 
bisección. De esta manera, los valores xl y xu siempre encierran la verdadera raíz. 
El  proceso  se  repite  hasta  que  la  aproximación  de  la  raíz  sea  adecuada.  El 
algoritmo es idéntico al de la bisección, excepto en la ecuación del método. 
 
Ejemplo 3.6: 
Utilice  el  método  de  Falsa  Posición  para  calcular  la  raíz  de  la  función 

f ( x) = e − x − x ,  empleando  como  valores  de  intervalo  xl=0,  xu=1.  Con  una 

aproximación de 4 cifras significativas. 
 
Solución: 

Iter.  xl  xu  f(xl)  f(xu)  f(xl)*f(xu)  xr  f(xr)  f(xl)*f(xr)  εa  

1  0  1  1.0000  ‐0.6321  ‐0.6321  0.61270  ‐0.0708  ‐0.07081    
2  0  0.61270  1.0000  ‐0.0708  ‐0.0708  0.57218  ‐0.0079  ‐0.00789  7.0814 
3  0  0.57218  1.0000  ‐0.0079  ‐0.0079  0.56770  ‐0.0009  ‐0.00088  0.7888 
4  0  0.56770  1.0000  ‐0.0009  ‐0.000877  0.56721  ‐0.0001  ‐0.00010  0.0877 
5  0  0.56721  1.0000  ‐0.0001  ‐0.000098  0.56715  ‐0.000011  ‐0.00001  0.0098 
6  0  0.56715  1.0000  0.0000  ‐0.000011  0.56714  ‐0.000001  ‐0.000001  0.0011 

Ejemplo 3.7: 
Desarrolle y pruebe un programa en MATLAB, que utilizando el Método de Falsa 

Posición,  permita  calcular  la  raíz de  la  función  f ( x) = e − x − x ,  así  mismo  el 
programa  debe  permitir  ingresar diferentes valores  del  intervalo  xl,  xu,  y  número 
de cifras significativas, para permitir así aproximar la raíz con diferente precisión. 
Finalmente probar el programa con los intervalo xl=0, xu=1 y con aproximaciones 
de 2, 3, 4, 5…15 cifras significativas. 
 
 
 
 

77
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Solución: 
Algoritmo del Método de Falsa Posición 
Leer: xl, xu, n 

Es=(0.5*10^(2‐n) 
Ea=100 
xr=0 
i=0 
Mientras Ea>Es hacer 
xa=xr 
xr=xu‐(f(xu)*(xl‐xu))/(f(xl)‐f(x)) 
f(xl)*f(xr)<0 

xu=xr  xl=xr 
Ea=abs((xr‐xa)/xr)*100  Ea=abs((xr‐xa)/xr)*100 

i=i+1 
Imprimir: xr, Ea, i 

Código del Método Falsa Posición (regulafalsi.m) 
clear;clc; 
disp('CALCULA LA RAIZ DE UNA ECUACION NO LINEAL') 
disp(' POR EL METODO CERRADO DE LA FALSA POSICION') 
n=input('Ingrese precision de cifras significativas n = '); 
xl=input('Ingrese limite inferior xl = '); 
xu=input('Ingrese limite superior xu = '); 
Es=(0.5*10^(2‐n)); 
Ea=100; 
xr=0; 
i=0; 
while Ea>Es 
    xa=xr; 
    xr=xu‐(funcion(xu)*(xl‐xu))/(funcion(xl)‐funcion(xu));%USO 
DE LA FORMULA DE METODO RAGULA FALSI 
    if funcion(xl)*funcion(xr)<0 
        xu=xr; 
        Ea=abs((xr‐xa)/xr)*100; 
    elseif funcion(xl)*funcion(xr)>0 
        xl=xr;   
        Ea=abs((xr‐xa)/xr)*100; 
    end 
    i=i+1; 
end 
fprintf('La raiz es : %12.15f\n',xr) 
fprintf('Error aproximado : %12.10f\n',Ea) 
fprintf('Numero de iteraciones : %12.0f\n',i) 

78
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Prueba del Programa 
>> regulafalsi 
CALCULA LA RAIZ DE UNA ECUACION NO LINEAL 
POR EL METODO CERRADO DE LA FALSA POSICION 
Ingrese precision de cifras significativas n = 4 
Ingrese limite inferior xl = 0 
Ingrese limite superior xu = 1 
La raiz es : 0.567144060375102 
Error aproximado : 0.0010850621 
Numero de iteraciones :      6 

3.3 Métodos abiertos 
Los  métodos  abiertos  en  contraste  a  los  métodos  cerrados,  se  basan  en  fórmulas 
que requieren únicamente de un solo valor de inicio x o que inicien con un par de 
ellos, pero que no necesariamente encierran la raíz. Estos algunas veces divergen o 
se alejan de la raíz verdadera a medida que se avanza con el cálculo. Sin embargo, 
cuando  los  métodos  abiertos  convergen,  en  general  lo  hacen  mucho  más  rápido 
que los métodos cerrados. 
 
3.3.1 Iteración simple de punto fijo 
Emplea  una  fórmula  para  predecir  la  raíz.  Esta  fórmula  se  desarrolla  como  una 
iteración simple de punto fijo (también llamada iteración de un punto o sustitución 
sucesiva o método de punto fijo), consiste en reordenar la ecuación f(x) = 0 de tal 
modo que x esté del lado izquierdo de la ecuación: 
x = g (x)  
Esta  transformación  se  realiza  mediante  operaciones  algebraicas  o  simplemente 
sumando x a cada lado de la ecuación original. 
La utilidad de la última ecuación, es que proporciona una fórmula para predecir un 
nuevo valor de x en función del valor anterior de x. De esta manera, dado un valor 
inicial  para  la  raíz  xi,  la  ecuación  anterior  se  utiliza  para  obtener  una  nueva 
aproximación xi+1, expresada por la formula iterativa 
xi +1 = g ( xi )  

El error aproximado de esta ecuación se calcula usando el error normalizado: 
xi +1 − xi
εa = × 100%  
xi +1

79
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Ejemplo 3.8: 
Use una iteración simple de punto fijo para localizar la raíz de 

f ( x) = e − x − x  
La función puede expresarse como: 

x = e− x  
O de modo con subíndices como: 

xi+1 = e − xi  
Empezando  desde  un  valor  xi=0,  se  calcula  xi+1  reiteradas  veces  hasta  la 
aproximación deseada. 
 

i  xi  xi+1  εa  
0  0  1  100.0 
1  1  0.3679  171.828 
2  0.3679  0.6922  46.854 
3  0.6922  0.5005  38.309 
4  0.5005  0.6062  17.447 
5  0.6062  0.5454  11.157 
M  M  M  M 
10  0.5649  0.5684  0.624 
M  M  M  M 
19  0.5671  0.56714771  0.004 

Ejemplo 3.9: 
Use una iteración simple de punto fijo para localizar la raíz de 
f ( x) = e − x − x  
 
Algoritmo del Método de Punto Fijo 
Leer: xi, n 
Es=(0.5*10^(2‐n) 
Ea=100 
i=0 

Mientras Ea>Es hacer 
xa=xi 
xi=exp(‐xi) 
Ea=abs((xi‐xa)/xi)*100 
i=i+1 
Imprimir: xi, Ea, i 
 

80
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Código del Método de Punto Fijo (puntofijo.m) 
clear; clc;    %Permite eliminar variables en uso 
disp('     CALCULA LA RAIZ DE ECUACION NO LINEAL') %Muestra el contenido 
disp('          POR EL METODO ABIERTO ')   %Muestra el contenido 
disp('      ITERACION SIMPLE PUNTO FIJO')   %Muestra el contenido 
n=input('Ingrese  precisión  de  cifras  significativas  n  =  ');  %Solicita  un 
valor n 
xi=input('Ingrese valor inicial xo = ');   %Solicita un valor xi 
Es=(0.5*10^(2‐n));  %Formula  que  permite  calcular  la  tolerancia  porcentual 
prefijada 
Ea=100;   %Asigna el valor de 100 a la variable Ea 
i=0;     %Asigna el valor de 0 a la variable i 
while Ea>Es  %Verifica si el valor de la variable Ea es mayor que la de Es 
    xa=xi;    %El dato de la variable xi es asignado a la variable xa 
    xi=exp(‐xi);    %Se aplica la formula de punto fijo 
    Ea=abs((xi‐xa)/xi)*100;  %Se  calcula  el  error  relativo  porcentual 
aproximado 
    i=i+1;    %Se adiciona en una unidad el número de contador 
end 
fprintf('La  raiz  es  :  %12.15f\n',xi)  %Imprime  xr  (raíz  buscada)  con  un 
formato 
fprintf('Numero de iteraciones : %12.0f\n',i)   %Imprime i (iteraciones) con un 
formato 
 
Corrida y prueba del programa “Iteración simple de punto fijo” 
 
>> puntofijo 
CALCULA LA RAIZ DE ECUACION NO LINEAL 
          POR EL METODO ABIERTO  
      ITERACION SIMPLE PUNTO FIJO 
Ingrese precisión de cifras significativas n = 8 
Ingrese valor inicial xo = 0 
La raíz es : 0.567143289516084 
Numero de iteraciones:           36 

3.3.2 Método de Newton Raphson 
El  método  de  Newton‐Raphson  es  uno  de  los  algoritmos  más  ampliamente 
utilizados  para  localizar  raíces.  La  fórmula  del  método  de  Newton  Raphson,  se 
deduce a partir de la siguiente interpretación geométrica. 

f(x)
Pendiente =f’(xi)

f(xi)

0
xi+1 xi x
 

81
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Si  el  valor  inicial  para  la  raíz  es  xi  ,  entonces  se  puede  trazar  una  tangente  en  el 
punto (xi ,f(xi)) de la curva: Por lo general, la abscisa xi+1 del punto de intersección 
de la tangente con el eje x, representa una mejor aproximación de la raíz. 
Del grafico se tiene: 
f ( xi ) − 0
f ' ( xi ) =  
xi − xi +1
Reordenando y despejando xi+1, se tiene: 
f ( xi )
xi +1 = xi −  
f ' ( xi )
Un método alternativo para deducir la formula de Newton‐Raphson, es la serie de 
Taylor de primer orden, el cual se escribe como: 

f ( xi +1 ) ≅ f ( xi ) + f ' ( xi ) ⋅ ( xi +1 − xi )  
En la intersección con el eje x, f(xi+1) debe ser igual a cero, es decir: 

0 = f ( xi ) + f ' ( xi ) ⋅ ( xi +1 − xi )  
De donde se puede despejar xi+1, así 
f ( xi )
xi +1 = xi −  
f ' ( xi )
Que es la relación similar al obtenido de la interpretación geométrica. 
 
Ejemplo 3.10: 
Utilice el método de Newton Raphson para localizar la raíz de   f ( x) = e − x − x = 0  y 
empleando  como  valor  inicial  o  semilla  xi=0.  Con  una  aproximación  de  4  cifras 
significativas. 
 
Solución: 
Se sabe que la derivada de la función f(x) es:  f ' ( x) = −e − x − 1 

i xi f(xi) f'(xi) xi+1 εa

0 0 1 -2.0000 0.500 100.00
1 0.5000000 0.10653065971 -1.6065 0.56631100319722 11.709290977
2 0.5663110 0.00130450981 -1.5676 0.56714316503486 0.146728708
3 0.5671432 0.00000019648 -1.5671 0.56714329040978 0.000022106

82
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

La raíz de la ecuación es 0.567143 con un error relativo porcentual aproximado de 
0.000022106%, que nos indica que es exacto hasta la sexta cifra significativa. 
 
Ejemplo 3.11: 
Desarrolle, depure y pruebe un subprograma en MATLAB que permita determinar 
la  raíz  de  f ( x) = e − x − x = 0 ,  utilizando  el  método  de  Newton  Raphson  y 
empleando como valor inicial o semilla xi=0. 
Algoritmo del Método de Newton – Raphson 
Leer: xi, n 
Es=(0.5*10^(2‐n) 
Ea=100 
i=0 

Mientras Ea>Es hacer 
xa=xi 
fxi=exp(‐xi)‐xi 
dxi=‐exp(‐xi)‐1 
xi=xi‐fxi/dxi 
Ea=abs((xi‐xa)/xi)*100 
i=i+1 

Imprimir: xi, Ea, i 

 
Código del Método de Newton Raphson (newton.m) 
clear;clc; 
disp('     CALCULA LA RAIZ DE ECUACION NO LINEAL') 
disp('          POR EL METODO ABIERTO ') 
disp('             NEWTON RAPHSON') 
n=input('Ingrese precision de cifras significativas n = '); 
xi=input('Ingrese valor inicial xi = '); 
Es=(0.5*10^(2‐n));  %tolerancia porcentual prefijada 
Ea=100;          %Error relativo porcentual aproximado inicial 
i=0; 
while Ea>Es 
    xa=xi; 
    fxi=exp(‐xi)‐xi;        %Funcion f(x) evaluada en xi 
    dxi=‐exp(‐xi)‐1;        %Derivada f(x) evaluada en xi 
    xi=xi‐fxi/dxi;          %Algoritmo de Newton Raphson 
    Ea=abs((xi‐xa)/xi)*100;  %Calcula  el  error  relativo 
porcentual aproximado 
    i=i+1;         %Cuenta el numero de iteraciones efectuadas 
end 
fprintf('La raiz es : %12.15f\n',xi)       %Imprime la raiz x 
fprintf('Numero de iteraciones : %12.0f\n',i)   %Imprime el num 
de iteraciones 
Corrida y prueba del programa “Newton Raphson”

83
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

>> newton 
     CALCULA LA RAIZ DE ECUACION NO LINEAL 
          POR EL METODO ABIERTO  
             NEWTON RAPHSON 
Ingrese precisión de cifras significativas n = 5 
Ingrese valor inicial xi = 0 
La raíz es: 0.567143290409781 
Numero de iteraciones:            4 

3.3.3 Método de la secante 
Esta  técnica  es  similar  a  la  del  método  de  Newton  Raphson.  Sin  embargo  un 
problema  potencial  en  la  implementación  del  método  de  Newton  Raphson,  es  la 
evaluación  de  la  derivada.  Aunque  no  es  un  inconveniente  para  los  polinomios  y 
muchas  otras  funciones,  existen  algunas  funciones  cuyas  derivadas  resultan  muy 
engorrosas  y  difíciles  de  calcular.  Sin  embargo,  para  salvar  esta  dificultad  el 
método de la Secante usa una diferencia finita divida hacia atrás en lugar de una 
derivada de la función para estimar la pendiente. 
La  representación  gráfica  del  método  de  la  secante  es  como  se  muestra  en  la 
siguiente figura. 

f(x) f ( xi −1 ) − f ( xi )
f ' ( xi ) ≅
xi −1 − xi

f(xi)

f(xi-1)

0
xi+1 x   xi-1 xi
La aproximación de la derivada f´(xi) está dado por la recta secante que corta a la 
función en (xi ,f(xi)) y (xi‐1 ,f(xi‐1)), siendo xi y xi‐1 puntos muy cercamos, obteniéndose 
la siguiente relación: 
f ( xi −1 ) − f ( xi )
f ' ( xi ) ≅  
xi −1 − xi
Esta  aproximación  de  la  derivada  por  la  diferencia  finita  divida  hacia  atrás,  se 
reemplaza  en  la  ecuación  o  algoritmo  de  Newton  Raphson,  para  obtener  la 
siguiente ecuación iterativa del método de la secante: 

84
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

f ( xi )( xi −1 − xi )
xi +1 = xi −  
f ( xi −1 ) − f ( xi )
Ejemplo 3.12: 
Utilice  el  método  de  la  Secante  para  localizar  la  raíz  de    f ( x) = e − x − x = 0   y 
empleando  como  valor  inicial  o  semilla  xi‐1=0    xi=1.  Con  una  aproximación  de  4 
cifras significativas. 

Iteración. xi-1 xi f(xi-1) f(xi) xi+1 ε a (%)

1 0 1 1 -0.63212056 0.61269984 63.21206
2 1 0.61270 -0.632121 -0.07081395 0.56383839 8.66586
3 0.61270 0.56384 -0.070814 0.00518235 0.56717036 0.58747
4 0.56384 0.56717 0.005182 -0.00004242 0.56714331 0.00477

Ejemplo 3.13: 
Desarrolle, depure y pruebe un subprograma en MATLAB que permita determinar 
la  raíz  de  f ( x) = e − x − x = 0 ,  utilizando  el  Método  de  la  Secante  y  empleando 
como valor inicial o semilla xi‐1=0 y xi=1. 
Algoritmo del Método de la Secante 

Leer: xi, xo, n 
Es=(0.5*10^(2‐n) 
Ea=100 
i=0 

Mientras Ea>Es hacer 
xa=xi 
fxo=exp(‐xo)‐xo 
fxi=exp(‐xi)‐xi 
xi=xi‐fxi*(xo‐xi)/(fxo‐fxi) 
Ea=abs((xi‐xa)/xi)*100 
xo=xa 
i=i+1 
Imprimir: xi, Ea, i 

Código del Método de la Secante (secante.m)

clear;clc; 
disp('     CALCULA LA RAIZ DE ECUACION NO LINEAL') 
disp('          POR EL METODO ABIERTO ') 
disp('             DE LA SECANTE') 
n=input('Ingrese precision de cifras significativas n = '); 
xo=input('Ingrese valor inicial xi‐1 = '); 
xi=input('Ingrese valor inicial xi = '); 
Es=(0.5*10^(2‐n));     %tolerancia porcentual prefijada 

85
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Ea=100;            %Error relativo porcentual aproximado inicial 
i=0; 
while Ea>Es 
    xa=xi; 
    fxo=exp(‐xo)‐xo;         %Funcion f(x) evaluada en xi‐1 
    fxi=exp(‐xi)‐xi;         %Funcion f(x) evaluada en xi 
    xi=xi‐fxi*(xo‐xi)/(fxo‐fxi);  %Algoritmo  del  Metodo  de  la 
Secante 
    Ea=abs((xi‐xa)/xi)*100;    %Calcula  el  error  relativo 
porcentual aproximado 
    xo=xa; 
    i=i+1;                    %Cuenta el numero de iteraciones 
efectuadas 
end 
fprintf('La raiz es : %12.15f\n',xi)          %Imprime la raiz x 
fprintf('Numero de iteraciones : %12.0f\n',i)  %Imprime el num 
de iteraciones 

Corrida y prueba del programa “Secante” 
>> secante 
     CALCULA LA RAIZ DE ECUACION NO LINEAL 
          POR EL METODO ABIERTO  
             DE LA SECANTE 
Ingrese precisión de cifras significativas n = 5 
Ingrese valor inicial xi‐1 = 0 
Ingrese valor inicial xi = 1 
La raíz es: 0.567143290409705 
Numero de iteraciones:            5 

3.4 Sistema de ecuaciones no lineales 
Consiste  en  obtener  las  raíces  de  un  conjunto  de  ecuaciones  simultáneas  no 
lineales 
f1 ( x1 , x 2 ,K , x n ) = 0
f 2 ( x1 , x 2 ,K , x n ) = 0
 
M M
f n ( x1 , x 2 ,K , x n ) = 0
La solución de este sistema de ecuaciones simultáneas, consiste de un conjunto de 
valores de xi que simultáneamente verifican todas las ecuaciones. 
Ejemplo:  Para  el  caso  de  dos  ecuaciones  dos  variables  haciendo  x1=x  y  x2=y  se 
tiene. 
x 2 + xy = 10  

y + 3 xy 2 = 57  

86
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Las que son dos ecuaciones simultaneas no lineales con dos incógnitas, las que se 
pueden expresar como: 
u ( x, y ) = x 2 + xy − 10 = 0  

v( x, y ) = y + 3xy 2 − 57 = 0  
La  mayoría  de  los  métodos  para  determinar  la  solución  de  este  sistema  de 
ecuaciones  no  lineales,  son  extensiones  de  los  métodos  abiertos  para  resolver 
ecuaciones simples como: iteración de punto fijo, Newton Raphson y otros. 
 
3.4.1 Iteración simple de punto fijo 
El  método  de  iteración  simple  de  punto  fijo  para  resolver  ecuaciones  simples, 
puede modificarse para resolver dos ecuaciones simultáneas no lineales. 
f ( x) = 0  
x = g (x)  

xi +1 = g ( xi )  

De modo similar, el algoritmo de iteración simple de punto fijo para dos ecuaciones 
simultaneas no lineales es: 

xi +1 = g ( xi , y i )
y i +1 = h( xi , y i )  
 
Ejemplo 3.14: 
Desarrolle, depure y pruebe un subprograma en MATLAB que permita determinar 
las raíces de las ecuaciones: 
u ( x, y ) = x 2 + xy − 10 = 0  

v( x, y ) = y + 3xy 2 − 57 = 0  
Utilizando el Método de Iteración de Punto Fijo y empleando como valor inicial o 
semilla xi=1.5 y yi=3.5 
Solución 
Despejando x de la ecuación:  u ( x, y ) = x 2 + xy − 10 = 0  

Y despejando y de la ecuación:  v( x, y ) = y + 3xy 2 − 57 = 0  
Se tiene: 
x = 10 − x ⋅ y   O con subíndices  xi +1 = 10 − xi y i  

87
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

57 − y 57 − yi
y=    O con subíndices  yi +1 =  
3x 3 xi

Utilizando  alternadamente  una  y  otra  de  las  ecuaciones  se  obtiene  la  siguiente 
tabla de resultados. 

Iteración (i)  xi  yi   xi+1  yi+1  ε a x (%)  ε a y (%) 

1  1.5  3.5  2.179449  2.860506 31.175  22.356 
2  2.179449  2.860506 1.940534  3.049551 12.312  6.199 
3  1.940534  3.049551 2.020456  2.983405 3.956  2.217 
4  2.020456  2.983405 1.993028  3.005704 1.376  0.742 
5  1.993028  3.005704 2.002385  2.998054 0.467  0.255 
6  2.002385  2.998054 1.999185  3.000666 0.160  0.087 
7  1.999185  3.000666 2.000279  2.999773 0.055  0.030 
8  2.000279  2.999773 1.999905  3.000078 0.019  0.010 
9  1.999905  3.000078 2.000033  2.999973 0.006  0.003 
10  2.000033  2.999973 1.999989  3.000009 0.002  0.001 
11  1.999989  3.000009 2.000004  2.999997 0.0007  0.0004 
12  2.000004  2.999997 1.999999  3.000001 0.0003  0.0001 
13  1.999999  3.000001 2.000000  3.000000 0.00009  0.00005 

Código en MATLAB
clear;clc; 
disp('CALCULA LA RAIZ DE UN SISTERMA DE DOS ECUACIONES NO LINEALES') 
disp('          POR EL METODO ABIERTO DE') 
disp('      ITERACION SIMPLE PUNTO FIJO') 
n=input('Ingrese precision de cifras significativas n = '); 
xo=input('Ingrese valor inicial xo = '); 
yo=input('Ingrese valor inicial yo = '); 
Es=(0.5*10^(2‐n)); 
Eax=100; 
Eay=100; 
i=0; 
%EL SISTEMA DE ECUACIONES ES 
%U(x,y)=x^2+x*y‐10=0 
%V(x,y)=y+3x*y^2‐57=0 
while Eax>Es & Eay>Es 
    xi=sqrt(10‐xo*yo);          %xi+1=g(xi,yi) 
    yi=sqrt((57‐yo)/(3*xi));    %yi+1=h(xi,yi) 
    Eax=abs((xi‐xo)/xi)*100;    %Calcula el error relativo aproximado 
para x 
    Eay=abs((yi‐yo)/yi)*100;    %Calcula el error relativo aproximado 
para y 
    xo=xi; 
    yo=yi; 
    i=i+1; 
end 
fprintf('La raiz "x" es : %12.15f\n',xi) 
fprintf('La raiz "y" es : %12.15f\n',yi) 
fprintf('Numero de iteraciones : %12.0f\n',i) 

88
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Corrida y prueba del programa “punto fijo”

>> puntofijo2 
CALCULA LA RAIZ DE UN SISTERMA DE DOS ECUACIONES NO LINEALES 
          POR EL METODO ABIERTO DE 
      ITERACION SIMPLE PUNTO FIJO 
Ingrese precisión de cifras significativas n = 5 
Ingrese valor inicial xo = 1.5 
Ingrese valor inicial yo = 1.5 
La raíz "x" es: 2.000002516109123 
La raíz "y" es: 2.999997946009110 
Numero de iteraciones:           13 

3.4.2 Newton Raphson 
El método de Newton Raphson para una sola ecuación, se baso en la Serie de 
Taylor de primer orden. 

f ( xi+1 ) ≅ f ( xi ) + f ' ( xi )( xi+1 − xi )  

Donde xi es el valor inicial de la raíz y xi+1 es el valor en el cual la recta tangente 
intersecta  al  eje  x,  la  que  por  definición  en  la  intersección  hace  f(xi+1)  es  igual  a 
cero. 
La ecuación anterior se reordena para obtener: 

0 = f ( xi ) + f ' ( xi )( xi+1 − xi )  
f ( xi )
xi +1 = xi −  
f ' ( xi )
Que  es  la  ecuación  de  Newton  Raphson  para  una  sola  ecuación  ya  visto 
anteriormente. 
La forma del algoritmo de Newton Raphson para múltiples ecuaciones se obtiene 
en  forma  idéntica.  Sin  embargo,  se  debe  usar  una  serie  de  Taylor  de  múltiples 
variables para tomar en cuenta el hecho de que más de una variable independiente 
contribuye a la determinación de la raíz. 
En  el  caso  de  dos  variables,  una  serie  de  Taylor  de  primer  orden  se  escribe  para 
cada ecuación no lineal como: 
∂u i ∂u
u i +1 = u i + ( xi +1 − xi ) + i ( y i +1 − y i )  
∂x ∂y
∂vi ∂v
vi +1 = u i + ( xi +1 − xi ) + i ( y i +1 − y i )  
∂x ∂y

89
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

De la misma manera como en el caso de una sola ecuación, la raíz aproximada 
correspondiente a los valores de x y y, donde  ui +1  y  vi +1  son iguales a cero, en tal 

sentido se reordena la ecuación anterior como: 
∂u i ∂u ∂u ∂u
xi +1 + i y i +1 = −u i + xi i + y i i  
∂x ∂y ∂x ∂y
∂vi ∂v ∂v ∂v
xi +1 + i y i +1 = −vi + xi i + y i i  
∂x ∂y ∂x ∂y
Debido a que se conoce todos los valores con subíndice i (correspondiente al 
último valor estimado), las únicas incógnitas son xi+1 y yi+1. Entonces, las dos 
ecuaciones anteriores son dos ecuaciones lineales con dos incógnitas. En 
consecuencia se puede usar manipulaciones algebraicas (por ejemplo la regla de 
Cramer) para resolverlo. 
Finalmente el algoritmo de Newton Raphson para dos ecuaciones no lineales 
simultáneas es: 
∂vi ∂u ∂u i ∂v
ui ⋅ − vi ⋅ i vi ⋅ − ui ⋅ i
∂y ∂y ∂x ∂x  
xi +1 = xi − ;       y i +1 = yi −
∂u i ∂vi ∂u i ∂vi ∂u i ∂vi ∂u i ∂vi
⋅ − ⋅ ⋅ − ⋅
∂x ∂y ∂y ∂x ∂x ∂y ∂y ∂x
Donde  al  denominador  de  cada  una  de  estas  ecuaciones,  se  conoce  formalmente 
como el determinante Jacobiano del sistema. 
Ejemplo 3.15: 
Desarrolle, depure y pruebe un subprograma en MATLAB que permita determinar 
las raíces de las ecuaciones: 
u ( x, y ) = x 2 + xy − 10 = 0  

v( x, y ) = y + 3xy 2 − 57 = 0  
Empleando  el  Método  de  Newton  Raphson  y  empleando  como  valor  inicial  o 
semilla xi=1.5 y yi=3.5

εa εa
i  xi  yi  ui  vi  ∂ui/∂x ∂ui/∂y ∂vi/∂x ∂vi/∂y Jacabiano xi+1 yi+1 x(%) y(%)
1  1.0  4.0  ‐5.0  ‐5.0  6.0  1.0  48.0  25.0  102.0  2.17647  1.94118  54.05  106.06 
2  2.17647  1.94118  ‐1.03806  ‐30.45  6.2941  2.1765 11.3045 26.3495 141.243  1.90083  3.21524  14.50  39.63 
3  1.90083  3.21524  ‐0.2752  5.1663  7.0169  1.9008 31.0133 37.6698 205.374  1.99913  2.99717  4.92  7.28 
4  1.99913  2.99717  ‐0.01177  ‐0.128  6.9954  1.9991 26.9490 36.9503 204.608  2.00000  3.00000  0.044  0.095 
5  2.00000  3.00000  3.2E‐06  9E‐05  7.0000  2.0000 27.0000 37.0000 205.000  2.00000  3.00000  0.00002 0.00009

90
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Código en MATLAB
clear;clc; 
disp('CALCULA LA RAIZ DE UN SISTERMA DE DOS ECUACIONES NO LINEALES') 
disp('          POR EL METODO ABIERTO DE') 
disp('                  NEWTON RAPHSON') 
n=input('Ingrese precision de cifras significativas n = '); 
xi=input('Ingrese valor inicial xi = '); 
yi=input('Ingrese valor inicial yi = '); 
Es=(0.5*10^(2‐n)); %tolerancia porcentual prefijada 
Eax=100;           %Error relativo porcentual aproximado para x 
Eay=100;           %Error relativo porcentual aproximado para y 
i=0; 
%EL SISTEMA DE ECUACIONES ES 
%U(x,y)=x^2+x*y‐10=0 
%V(x,y)=y+3x*y^2‐57=0 
while Eax>Es & Eay>Es %Verifica si los errores Eax y Eay cumplen con la Es 
    xa=xi; 
    ya=yi; 
    Ui=xi^2+xi*yi‐10;   %Funcion U(x,y) evaluada en xi, yi 
    Vi=yi+3*xi*yi^2‐57; %Funcion V(x,y) evaluada en xi, yi 
    dUidx=2*xi+yi;      %Derivada U(x,y) respecto de x, evaluada en xi, yi 
    dUidy=xi;           %Derivada U(x,y) respecto de y, evaluada en xi, yi 
    dVidx=3*yi^2;       %Derivada V(x,y) respecto de x, evaluada en xi, yi 
    dVidy=1+6*xi*yi;    %Derivada V(x,y) respecto de y, evaluada en xi, yi 
    Jacobiano=dUidx*dVidy‐dUidy*dVidx;   %Determinante Jacobiano del Sistema 
    xi=xi‐(Ui*dVidy‐Vi*dUidy)/Jacobiano; %Algoritmo de Newton Raphson para 
xi 
    yi=yi‐(Vi*dUidx‐Ui*dVidx)/Jacobiano; %Algoritmo de Newton Raphson para 
yi 
    Eax=abs((xi‐xa)/xi)*100;             %Calcula el error relativo 
aproximado para x 
    Eay=abs((yi‐ya)/yi)*100;             %Calcula el error relativo 
aproximado para y 
    i=i+1;                               %Cuenta el numero de iteraciones 
efectuadas 
end 
fprintf('La raiz "x" es : %12.15f\n',xi)        %Imprime la raiz x 
fprintf('La raiz "y" es : %12.15f\n',yi)        %Imprime la raiz y 
fprintf('Numero de iteraciones : %12.0f\n',i)   %Imprime el numero de 
iteraciones

Corrida y prueba del programa “Newton Raphson” 
>> newton2 
CALCULA LA RAIZ DE UN SISTERMA DE DOS ECUACIONES NO LINEALES 
          POR EL METODO ABIERTO DE 
                  NEWTON RAPHSON 
Ingrese precisión de cifras significativas n = 5 
Ingrese valor inicial xi = 1.5 
Ingrese valor inicial yi = 3.5 
La raíz "x" es: 1.999999999999980 
La raíz "y" es: 3.000000000000075 
Numero de iteraciones:            4 
 
 

91
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

3.5 Localización de raíces con software 
Los paquetes o software de cómputo tienen gran capacidad para localizar las raíces 
de ecuaciones, en este texto se presentan los más comunes. 
3.5.1 EXCEL 
La hoja de cálculo Excel, se utiliza para localizar la raíz mediante prueba y error. 
−x
Por  ejemplo  si  se  quiere  encontrar  la  raíz  de  f ( x) = e − x = 0 ,  primero  se 
introduce un valor de x en una celda. Después se destina otra celda para f(x) donde 
se obtendrá el valor de la función para la x de la primera celda. Se puede variar el 
valor de la celda en x hasta que la celda en f(x) se aproxime a cero. Este proceso se 
mejora  usando  la  capacidad  de  graficación  de  Excel  para  obtener  un  buen  valor 
inicial. 
Aunque  Excel  facilita  el  método  de  prueba  y  error,  también  posee  dos 
herramientas estándar que sirve para la localización de raíces: Buscar objetivo y 
Solver.  Ambos  son  útiles  para  ajustar  sistemáticamente  los  valores  iniciales, 
Buscar objetivo se utiliza expresamente para llevar la ecuación a un valor (en este 
caso, cero) mediante la variación de un solo parámetro. 
Ejemplo 3.16: 
Haciendo uso de la herramienta Buscar objeto (Goal Seek) de Excel, determinar la 
−x
raíz de la Ecuación  f ( x) = e − x = 0  
Solución: 
Pasos a seguir: 
1. Con  Excel  crear  una  celda  que  tenga  el  valor  de  la  función  en  cuestión  y 
hacer,  después,  el  valor  dependiente  de  otra  celda.  Para  nuestro  ejemplo 
A11 y B11. 

92
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

2. Una  vez  hecho  esto,  del  menú  herramientas  seleccionar  “Buscar  objetivo”. 
Ahora  aparece  una  ventana  de  diálogo  pidiendo  se  especifique  una  celda 
para  un  valor  al  modificar  otra  celda.  Para  nuestro  ejemplo,  el  valor 
propuesto se escribe en la celda A11 y la función resultante en la celda B11. 
La ventana de diálogo “Buscar objetivo” será: 
 

 
Al hacer click en Aceptar, se muestra una ventana con el valor actual de la función. 
 

 
 
Al  aceptar,  la  celda  A11  de  la  hoja  de  de  cálculo,  se  ha  modificado  con  el  nuevo 
valor de x (raíz), que hace muy próximo a cero la función evaluada. 
 

 
Siendo la solución de la raíz x=0.5671… que hace la función f(x) ≅ 0. 
 

93
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Ejemplo 3.17: 
Use la herramienta Solver de Excel, para determinar la raíz de la Ecuación 
f ( x) = e − x − x = 0  
Solución. 
Pasos a seguir: 
1. Tabular  la  función  f(x)  para  diferentes  valores  de  x,  de  tal  forma  que  la 
función f(x) pase de un valor positivo (+) a otro negativo (‐). 
2. Graficar  la  función,  para  visualizar  e  identificar  el  valor  más  próximo  a  la 
raíz solución, como el mostrado en la siguiente figura. 

 
 
3. Cargar  el  complemento  Solver  de  Excel,  para  lo  cual  en  el  menú 
Herramientas, haga click en Complementos… 

94
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

4. En  el  cuadro  Complementos  disponibles,  active  la  casilla  de  verificación 
que aparece junto a Solver y haga clic en Aceptar. 

 
5. Ejecutar Solver. Para lo cual en el menú Herramientas, haga clic en 
Solver… 

 
6. Las celdas solicitadas en la ventana Solver se llenan como: 

95
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Al hacer click en el botón Resolver, se mostrará una ventana de dialogo como: 

 
 
Al hacer click en el botón Aceptar, se mostrará la ventana de Excel con la celda 
A11 modificado, mostrando la raíz solución de la ecuación como: 

 
Solución de x=0.567143351958934 
 
Ejemplo 3.18: 
Resolver un Sistema de dos Ecuaciones No Lineal, con herramienta Solver de Excel. 
u ( x, y ) = x 2 + xy − 10 = 0  

v( x, y ) = y + 3xy 2 − 57 = 0  
Procedimiento a seguir: 
1. Ingresar  los  valores  iniciales  de  “x”  e  “y”  en  las  celdas  B1  y  B2  de  Excel 
respectivamente. 
2. Evaluar las funciones u(x,y) y v(x,y) respectivamente en las celdas B3 y B4 en 
función de “x” y “y” ubicadas en las celdas B1 y B2, como se muestra en la figura 
3. En la celda B6, hacer una sumatoria de cuadrados de las funciones u(x,y) y 
v(x,y), a fin de que la celda B6 sea la que agrupe los resultados de las dos 
funciones evaluadas. 

96
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

 
 
4. Ejecutar Solver, como se muestra en la figura. 

 
 
5. Finalmente se obtiene las raíces del sistema de ecuaciones, como: 
 

 
Respuesta 
x=2.000029761   ≅ 2 
y=2.999983754   ≅ 3 
 
 
 

97
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

3.5.2 MATLAB 
MATLAB es capaz de localizar raíces en ecuaciones algebraicas y trascendentales, 
Siendo excelente para la manipulación y localización de raíces en los polinomios. 
La  función  fzero  está  diseñada  para  localizar  la  raíz  de  una  función.  Una 
representación simplificada de su sintaxis es: 
fzero(f, Xo, opciones) 
Donde  f  es  la  función  que  se  va  a  analizar,  Xo  es  el  valor  inicial  y  opciones  son 
parámetros  de  optimización  (estos  pueden  cambiarse  al  usar  la  función 
optimset).  Si  no  se  anotan  las  opciones  se  emplean  los  valores  por  omisión. 
Observe que se pueden emplear uno o dos valores iniciales, asumiendo que la raíz 
esta dentro del intervalo. 
 
Ejemplo 19: 
Utilice la función fzero del MATLAB para encontrar la raíz de la siguiente función. 
1 1
f ( x) = + − 6 
( x − 0.3) + 0.01 ( x − 0.9) 2 + 0.04
2

Dentro del intervalo xl=‐1 y xu=2. 
 
Solución: 
Lo  primero  que  podemos  hacer  en  MATLAB  es  graficar  la  función,  a  fin  de  ver  el 
comportamiento  de  la  misma,  para  lo  cual  en  la  ventana  de  comandos  del 
MATLAB se debe digitar lo siguiente: 
 
>> x=‐1:0.01:2; 
>> y=1./((x‐.3).^2+.01)+1./((x‐.9).^2+0.04)‐6; 
>> plot(x,y); grid; 
 
Las  tres  anteriores  líneas  de  código  permiten  hacer  lo  siguiente:  la  primera  línea 
permite crear la variable x desde ‐1 hasta 2 cada 0.01 intervalos, la segunda línea 
permite  evaluar  la  variable  x  en  la  función  f(x)  llamado  en  este  caso  como  y,  y 
finalmente la tercera línea permite graficar las variables x y  y haciendo uso de la 
función plot(x,y) de MATLAB y el ultimo código grid permite dale un grillado al 
gráfico, luego de estas líneas el gráfico debe quedar como: 

98
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Raíces

 
Ahora  podemos  ver  que  la  función  f(x)  corta  al  eje  x  en  dos  puntos  que  hacen  la 
función  cero,  por  lo  cual  existen  dos  raíces  en  el  intervalo  ‐1  y  2,  uno  de  ellos 
próximo a cero (0) y el otro próximo a uno (1). 
 
Para encontrar las raíces se usa función fzero y se registra su resultado en las 
variables x1 y x2 respectivamente. 
>> x1=fzero('1./((x‐.3).^2+.01)+1./((x‐.9).^2+0.04)‐6',0) 
x1 = 
   ‐0.1316 
>> x2=fzero('1./((x‐.3).^2+.01)+1./((x‐.9).^2+0.04)‐6',1) 
x2 = 
    1.2995 
Ahora podemos evaluar esta raíces, es decir x1 y x2 en la función f(x), tal como se 
muestra a continuación. 
 
>> y1=1./((x1‐.3).^2+.01)+1./((x1‐.9).^2+0.04)‐6; 
>> y2=1./((x2‐.3).^2+.01)+1./((x2‐.9).^2+0.04)‐6; 

Finalmente podemos graficar estos resultados, conjuntamente con la función

con el siguiente código.

99
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

>> plot(x,y,x1,y1,'*',x2,y2,'+');grid 
 
Luego del cual, se muestra el siguiente grafico. 

 
Siendo finalmente los resultados 
x1= ‐0.1316 
x2= 1.2995 
 

100
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

3.6 Problemas propuestos 
Métodos cerrados 
1. Escriba  un  programa  script  que  calcule  las  raíces  reales  de  una  función 
cuadrática  0. Llame al fichero raicescuad. Cuando el fichero se 
ejecute, éste debe pedir al usuario que introduzca los valores de las consonantes 
a, b y c. Para las raíces de la ecuación, el programa calculará el discriminante D: 

 4  
• Si D > 0, el programa visualizará un mensaje del tipo:”La ecuación tiene dos 
raíces”, y los valores de las raíces se visualizarán en la línea siguiente. 
• Si D = 0, el programa visualizará un mensaje del tipo: “La ecuación tiene una 
raíz”, y el valor de la raíz se visualizará en la línea siguiente. 
• Si D < 0, el programa visualizará un mensaje del tipo: “La ecuación no tiene 
raíces reales”. 
Ejecute el fichero script en la Ventana de Comandos tres veces para calcular las 
soluciones de las siguientes ecuaciones: 
a) 2 8 3 0 
b) 15 10 5 0 
c) 18 12 2 0 
2. Determine las raíces reales de 
26 82.3 88 45.4 9 0.65 : 
a) Gráficamente 
b) Usando  el  método  de  bisección  para  localizar  la  raíz  más  grande  con 
10%. Utilice como valores iniciales  0.5    1.0. 
c) Realice el mismo cálculo que en b), pero con el método de la falsa posición y 
0.1%. 
3. Determine la raíz real de ln  0.7: 
a) Gráficamente 
b) Empleando  tres  iteraciones  en  el  método  de  bisección  con  los  valores 
iniciales  0.5    2 
c) Usando tres iteraciones de la falsa posición, con los mismos valores iniciales 
de b). 
4. Determine la raíz real de  0.9 0 4 / : 

101
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

a) Analíticamente 
b) Gráficamente 
c) Empleando tres iteraciones en el método de la falsa posición, con los valores 
iniciales de 1 a 3. Calcule el error aproximado   y el error verdadero   en 
cada iteración. 
5. La velocidad v de caída de un paracaidista está dada por  

/
1  

Donde  g  =  9.8.  Para  el  paracaidista  con  un  coeficiente  de  arrate  c  =  14  Kg/s, 
calcule la masa m de éste de tal forma que la velocidad sea de v = 35 m/s en t = 7 
s. Con el método de la falsa posición determine m a un nivel de  0.1 %. 
6. Una  viga  se  carga  como  se  muestra  en  la  figura  P6.  Utilice  el  método  de  la 
bisección para encontrar la posición en la viga en la que no hay momento. 

100 lb/ft 100 lb

3’ 3’ 2’ 4’
 
Figura P6 
Métodos abiertos 
7. Con el método de iteración simple de punto fijo localice la raíz de 

√  
Use un valor inicial de  0.5 y haga iteraciones hasta que  0.001 %. 
8. Determine las raíces de  2.0 6 4 0.5 : 

a)  Gráficamente  y  b)  usando  el  método  de  Newton‐Raphson  que  cumpla  con 
  0.01% 
 
9. Localice la primera raíz positiva de 

102
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

  cos 1 1 
Donde x está en radianes. Use cuatro iteraciones con el método de la secante con 
valores  iniciales  de  a)  1.0  y  3.0,  y  b)  1.5  y  2.5,  para 
localizar la raíz. c) Use el método gráfico para verificar los resultados. 
10. Calcule  las  raíces  de  las  siguientes  ecuaciones  simultáneas  no  lineales 
usando  a)  el  método  de  iteración  de  punto  fijo  y  b)  el  método  de  Newton‐
Raphson: 

0.5 
5  
Emplee los valores iniciales de x = y = 1.0 y analice los resultados. 
 
11. Determine las raíces de las ecuaciones simultáneas no lineales 

4 4 4 
16 
Utilice  una  aproximación  grafica  para  obtener  los  valores  iniciales.  Determine 
una  mejor  aproximación  con  el  método  de  Newton‐Raphson  usando  las  dos 
ecuaciones descritas en la sección 3.4.2. 
12. a)  Aplique  el  método  de  Newton‐Raphson  a  la  función  tanh  9  
para evaluar su raíz real ya conocida en x = 3. Con valor inicial use  3.1 y 
realice un mínimo de cuatro iteraciones. 

b)  ¿El método muestra convergencia hacia su raíz real? Dibuja una  gráfica con 

los resultados para cada iteración efectuada. 
13. El polinomio 
0.0074 0.284 3.355 12 183 5 
Tiene una raíz real entre 15 y 20. Aplique a esta función el método de Newton‐
Rephson usando como valor inicial  16.15. Explique sus resultados. 
14. Utilice  el  método  de  la  secante  con  la  función  que  describe  un  círculo 
1 2 16  para  encontrar  una  raíz  real  positiva.  Tome  como 
valor  inicial  3 y  0 5.  Aproxime  la  solución  del  primer  y  cuarto 
cuadrante.  Al  encontrar  la  solución  para  f(x)  en  el  cuarto  cuadrante  asegúrese 
de tomar el valor negativo de la raíz cuadrada. ¿Por qué la solución diverge? 

103
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

15. v
La velocidad   de caída de un paracaidista esta dado por: 

g ⋅m
v=
c
(
1 − e − ( c / m ) ⋅t ) 
Donde  g=9.8.  Para  el  paracaidista  con  un  coeficiente  de  arrastre  c=  14  kg/s, 

calcule la masa m de éste de tal forma que la velocidad sea de  =35 m/s en t=7  v
s. Con el método de la secante, determine  m a un nivel de Es=0.1%. Es posible 
que  la  masa  m  se  encuentre  entre  50  y  100  Kg,  si  es  así  utilizar  estos  valores 
como iniciales para el método indicado. 
16. Para hallar la raíz cúbica de un número  a , se puede aplicar el algoritmo de 
Newton Raphson: 

Como  x = 3 a   ;  x = a 1 / 3  

Entonces  f ( x) = x 3 − a = 0 ;  f ' ( x) = 3x 2  
Reemplazando en el algoritmo de Newton 

1⎛ a ⎞
xi +1 = ⎜⎜ 2 ⋅ xi + 2 ⎟⎟  
3⎝ xi ⎠

Utilizando  el  diagrama  de  flujo  de  la  figura,  para  determinar  la  raíz  cúbica  del 
número 10 a un nivel de Es = 0.1% 
Inicio

xi=1

xi=(2*xi+a/xi2)/3

V Ea>Es

F
Imprime: xi

Fin
 
Aplicación a la ingeniería 
17. El desplazamiento de una estructura está definido por la siguiente ecuación 
para una vibración amortiguada: 

8 cos  

104
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

Donde  k  =  0.5  y  3.  a)  Utilice  el  método  gráfico  para  realizar  una  primera 
estimación del tiempo requerido para que el desplazamiento disminuya a 4. b) 
Use el método de Newton‐Raphson para calcular la raíz con  0.01 %. c) Con 
el método de la secante determine la raíz aj  0.01 %. 
18. La  fórmula  que  define  la  fuerza  por  unidad  de  área,  P/A,  que  causa  un 
máximo esfuerzo   en una columna que tiene una relación de esbeltez  /  es: 

 
1 /   0.5 / /
Si  200 000 kPa, ⁄ 0.2 y  250  kPa,  calcule  P/A  para  L/r  =100. 
Recuerde que sec 1/ cos . 
19. Un cable catenario es aquel que cuelga entre dos puntos que no están en la 
misma línea vertical. Como se muestra en la figura P19a, no está sujeto a otras 
cargas más que su propio peso. Así que su peso w (N/m) actúa como una carga 
uniforme  por  unidad  de  longitud  a  lo  largo  del  cable.  Un  diagrama  de  cuerpo 
libre de una sección  AB se representa en la figura P19b, donde  TA y  TB son las 
fuerzas  de  tensión  en  sus  extremos.  Considerando  los  equilibrios  de  fuerzas 
horizontal  y  vertical,  se  puede  obtener  el  siguiente  modelo  de  ecuación 
diferencial del cable. 

1  

B TB

A
w y0 W = ws
TA

x
a) b)
 
Figura P19. a) fuerza que actúa en una sección AB de un cable flexible colgante. La 
carga  es  uniforme  a  lo  largo  del  cable  (pero  no  es  uniforme  a  lo  largo  de  la 
distancias horizontal x). b) Un diagrama de cuerpo libre de la sección AB. 
 

105
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

20. En  la  figura  P20a    se  muestra  una  viga  uniforme  sujeta  a  una  carga 
distribuida  creciente  linealmente.  La  ecuación  para  calcular  la  curva  elástica 
resultante es (véase la figura P20b) 

2       (P20) 

Utilice el método de la bisección para determinar el punto de máxima deflexión 
(es  decir,  el  valor  de  x  donde  dy/dx  =  0).  Después,  sustituya  este  valor  en  la 
ecuación  (P8.18)  para  determinar  el  valor  de  máxima  deflexión.  Use  los 
siguiente  valores  de  los  parámetros  en  sus  cálculos:  L=  450  cm,  E  =  50  000 
Kn/cm2, I = 30 000 cm4 y  1.75 KN/cm. 

w0

a)

(x = L, y = 0)
(x = 0, y = 0)
x

b)  
Figura P20 
21. Una viga que descansa libremente en dos soportes, recibe una carga como 
se  muestra  en  la  figura  P21.  Usando  funciones  de  singularidad,  se  puede 
expresar el esfuerzo cortante a lo largo de la viga mediante la ecuación. 

20 0 5 15 8 57 
La función de singularidad se define como: 
  donde 
 
0                donde 
Encuentre el (los) punto(s) el esfuerzo cortante sea igual a cero. 
 

106
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

20 kips/ft

150 kip-ft
15 kips

5’ 2’ 1’ 2’
 
Figura P21 
 
22. En  la  hidráulica  de  canales es  frecuente,  calcular  el  tirante  crítico del  flujo 
con  fines  a  análisis  y  diseño.  Si  se  sabe  que  el  régimen  critico  en  canales 
trapezoidales está gobernado por las siguientes ecuaciones: 

A3 Q 2 2 ⋅ z ⋅ A3
Función  f ( y ) = −   Derivada  f ' ( y ) = 3 ⋅ A 2 −  
T g T2
Área: A = (b + z ⋅ y ) y ;    Espejo;  T = b + 2 ⋅ z ⋅ y  
Utilizando  el  Algoritmo  de  Newton  Raphson.  Determinar  el  tirante  crítico  “y” 
que se encuentra implícito en las ecuaciones del régimen crítico, para el cálculo 
iniciar de y=1.0 m y utilizar los siguientes datos: Q= Caudal de 150.4 m3/s, b= 
base de 7.4 m y z= talud de 3, donde Área (A) y Espejo (T) están en función de: 
y, b y z. 
23. La formula de Manning para un canal trapezoidal se escribe como: 

S 1 / 2 ((b + z ⋅ y ) y )
5/3
Q=  
(
n b + 2 ⋅ y z2 +1 2/3 )
Donde Q = flujo (m3/s), S = pendiente (m/m), n = coeficiente de rugosidad, b = 
ancho (m), y=tirante normal (m) y z = talud. Aplicando la formula de Manning y 
el  Algoritmo  de  la  Secante.  Determinar  el  tirante  normal  al  centímetro, 
iniciando los cálculos con los valores de yi‐1= 1.0 y yi=1.01, si se conoce que el 
acho  b=1.5  m,  el  talud  z=1.5,  el  coeficiente  de  rugosidad  n=0.014,  pendiente 
S=0.002 m/m y el canal conduce un caudal de 2 m3/s. 
24. El factor de fricción f de Darcy es utilizado para obtener la perdida de carga 
en  tuberías  con  flujo  de  fluido.  Determinar  el  factor  “f”,  utilizar  la  ecuación  de 
Colebroock­White y el Algoritmo de la Secante, con valores iniciales de f de 0.01 y 
0.02  con  una  aproximación  de  er ≤ 0.001% .  Para  el  cálculo  de  f  utilizar  los 

107
 CAPÍTULO 3: SOLUCIÓN DE ECUACIONES NO LINEALES

siguientes datos: Q= caudal de agua de 0.042 m3/s, Ks= rugosidad de tubería de 
PVC  es  1.5x10 -6 m ,  d=diámetro  de  0.1524  m.  y  ν = viscosidad  cinemática  del 
agua es  1.14 × 10 -6 m 2 /s . 

La ecuación de Colebroock‐White a utilizar es: 

1 ⎛ Ks 2.51 ⎞
= −2 log10 ⎜ + ⎟ 
⎜ ⎟
f ⎝ 3.71 ⋅ d Re⋅ f ⎠
Adicionalmente  Re  =  Número  de  Reynolds  (adimensional  y  constante) 
V ⋅d
Re = .  La  ecuación  de  continuidad  es: Q = A ⋅ V   donde  A  =  área  y  V= 
v
velocidad del flujo en la tubería. 
 
Formulas  adicionales  de  la  función  y  derivada  de  la  formula  de  Colebroock‐
White: 
1
Si hacemos  x =  
f

⎛ Ks 2.51x ⎞
g ( x) = −2 log 10 ⎜ + ⎟− x 
⎝ 3.71d Re ⎠

⎛ 2.51 ⎞
⎜ ⎟
2 ⎜ Re ⎟ −1 
g ' ( x) = −
ln 10 ⎜ Ks 2.51x ⎟
⎜ + ⎟
⎝ 3.71d Re ⎠
Algoritmo de Newton – Raphson 
g ( xi ) − xi
xi +1 = xi −  
g ' ( xi ) − 1

108