UNIVERSIDAD CATÓLICA DE TEMUCO

FACULTAD DE INGENIERÍA
DEPTO. DE CIENCIAS MATEMÁTICAS Y FÍSICAS

Asignatura : Cálculo Numérico, MAT-1123.

Profesor : Emilio Cariaga L.
Periodo : 1er. Semestre 2018.

APROXIMACIÓN DE CEROS DE FUNCIONES

Definición: Sea f una función real de variable real, y sea x̃ perteneciente

al dominio de f . Se dice que x̃ es un cero de f ssi f (x̃) = 0.

En este capı́tulo abordaremos el problema que consiste en calcular una

aproximación del cero de una función, pues no siempre es posible resolver
analı́ticamente la ecuación f (x) = 0, esto es, despejar la incógnita en términos
de los datos del problema.
En todo lo que sigue usaremos indistintamente las frases aproximar el cero
de una función, o resolver la ecuación.

Ejemplo: para la función cuadrática f (x) = ax2 + bx + c, con a 6= 0,

se puede plantear el problema de determinar sus ceros, esto es, resolver la
ecuación f (x) = 0. Como se sabe los ceros de la función cuadrática están
dados por p
b ± b2 4ac
x1,2 = .
2a
Los valores anteriores, esto es, los ceros x1 y x2 se pueden denominar
soluciones analı́ticas de la ecuación f (x) = 0. El problema que motiva este
capı́tulo es que la posibilidad de despejar la incógnita de una ecuación es
bastante excepcional. Es muy frecuente que debamos conformarnos sólo con
una aproximación. Sin embargo, veremos que los métodos numéricos nos
proveen de soluciones aproximadas de excelente calidad.

En la siguiente sección presentamos el método de Newton-Raphson, el

cual es ampliamente utilizado para resolver el problema numérico planteado.
Por supuesto, existe una amplia gama de métodos para resolver el problema
de aproximar un cero de una función uni y multivariada. Sin embargo, en
este curso nos focalizamos en este único método en razón de sus créditos, y
por ser muy utilizado en variadas librerias computacionales.

1
1 Método de Newton-Raphson.
Isaac Newton (1641-1727), Joseph Raphson (1648-1715).

Sea f una función al menos dos veces diferenciable sobre un intervalo

I = [a, b], sobre el cual se sabe que existe x̃ 2 I tal que f (x̃) = 0.
Sea xn 2 I, con n = 0, 1, 2, 3, ..., un número real muy cercano a x̃. Sea
x 2 I un número real vecino con xn . La aproximación lineal L(x) de f (x)
construida en torno a xn está dada por

L(x) = f (xn ) + f 0 (xn )(x xn ).

Asumiendo f 0 (xn ) 6= 0, para cada n = 0, 1, 2, ..., se tiene que

L(x) = 0
ssi

f (xn )
x = xn .
f 0 (xn )
Esta última igualdad motiva la siguiente definición del Método de New-
ton:
Método de Newton-Raphson: sea f 2 C 2 [a, b] para la cual se sabe
que axiste x̃ 2 [a, b] tal que f (x̃) = 0. Suponga que f 0 (x) 6= 0, para cada

2
x 2 [a, b]. La sucesión de números reales {xn ; n = 0, 1, 2, ...} ⇢ [a, b]:

f (xn )
xn+1 = xn ,
f 0 (xn )

cuando se usa para aproximar el valor de x̃, se denomina aproximación o

método de Newton-Raphson.

Los problemas matemáticos asociados a la definición de una sucesión de

números reales {xn ; n = 0, 1, 2, ...}, que será utilizada para aproximar x̃, son:

1. Problema de Convergencia: este problema consiste en establecer las

condiciones matemáticas bajo las cuales la sucesión {xn ; n = 0, 1, 2, ...}
es convergente a x̃, esto es, bajo cuáles hipótesis

lim xn = x̃.
n!1

2. Problema de Estimar el Error: este problema consiste en estimar

el error de aproximación en = x̃ xn ; n = 0, 1, 2, ....

El siguiente teorema responde las dos preguntas anteriores:

Teorema: sea f 2 C 2 [a, b], siendo x̃ 2 [a, b] un cero simple de f . Existe

una vecindad V de x̃, V ⇢ [a, b], y una constante C > 0, tales que si x0 2 V ,
entonces,
(1) La iteración de Newton-Raphson {xn ; n = 0, 1, 2, ...} converge a x̃, y
(2) El error de aproximación em = x̃ xm ; m = 0, 1, 2, ..., satisface, para todo
n = 0, 1, 2, ...,
|en+1 |  Ce2n .
El siguiente teorema responde ambos problemas mencionados, con inde-
pendencia del valor inicial x0 utilizado:

Teorema: si f 2 C 2 (R2 ), es creciente, convexa y tiene un cero, entonces

el cero es único y la iteración de Newton convergerá a él a partir de cualquier
punto inicial.

1.1 Ejemplos:
p
1. Construya una aproximación de 2 utilizando el método de Newton-
Raphson.

3
 p
Solución: sea x = 2, a partir de lo cual x2 = 2, ó x2 2 = 0. Esto
2
motiva definir
p la función f (x) = x 2, con lo cual el problema de
aproximar 2 se reduce al problema de calcular o aproximar un cero
de la función f (x) = x2 2.
La fórmula iterativa de Newton-Raphson está dada por

f (xn )
x1+n = xn ; n = 0, 1, 2, ...,
f 0 (xn )
o sea, para n = 0, 1, 2, ...

x2n 2 x2 + 2
x1+n = xn = n .
2xn 2xn
Si se elige como valor inicial x0 = 1 se obtienen los siguientes valores
numéricos. Puede notar que los valores obtenidos convergen rápidamente
al valor exacto.

n xn
0 1.000000000000000
1 1.500000000000000
2 1.416666666666667
3 1.414215686274510
4 1.414213562374690
5 1.414213562373095
6 1.414213562373095

Se concluye que

lim xn = 1.41421356237309...,
n!1

ó
p
2 = 1.4142135623730...

2. Se sabe que el polinomio p(x) = 4x3 2x2 + 3 posee una raiz x̃ en el

intervalo [ 2, 1]. Se pide aproximarla utilizando el método de Newton-
Raphson. Utilice x0 = 1 como valor inicial.
Solución: en este caso para n = 0, 1, 2, ...

4
 f (xn )
x1+n = xn
f 0 (xn )
4x3n 2x2n + 3
= xn
12x2n 4xn
8x3n 2x2n 3
= .
12x2n 4xn

Si se elige como valor inicial x0 = 1 se obtienen los siguientes valores

numéricos. Puede notar que los valores obtenidos convergen rápidamente
al valor exacto.

n xn
0 1.000000000000000
1 0.812500000000000
2 0.770804195804196
3 0.768832384255760
4 0.768828085869608
5 0.768828085849211
6 0.768828085849211

Se concluye que

lim xn = 0.76882808584921...,
n!1

ó

x̃ = 0.76882808584921...

1.2 Comentarios Adicionales:

1. En términos prácticos la recomendación es iniciar el análisis del prob-
lema f (x) = 0 graficando la función f en si misma lo cual en la ac-
tualidad no presenta ninguna dificultad dada la gran cantidad de pro-
gramas computacionales comerciales y gratuitos disponibles. Incluso
en algunos casos es posible obtener uno o dos decimales significativos
haciendo zoom en el gráfico.

5
 2. Una pregunta práctica que surge al momento de utilizar una sucesión
como herramienta matemática de aproximación es cuál elemento de
ésta utilizar como aproximación, o dicho de otra manera cuándo dejar
de iterar. Algunos de los criterios prácticos utilizados habitualmente en
diversos programas computacionales consisten en dada una tolerancia
✏ > 0 dejar de iterar cuando (n = 0, 1, 2, 3, ...):

(a) |xn+1 xn | < ✏

|xn+1 xn |
(b) |xn |
· 100% < ✏

(c) f (xn+1 ) < ✏.

Usualmente (según Scarborough (1966), citado por Chapra-Canale (2011,
6ta ed.)), si se selecciona ✏ como ✏% = 0.5⇤102 d , en donde d = 1, 2, 3, ...
es el número de dı́gitos o cifras significativos que al menos deberı́a
poseer la aproximación xn+1 . En general, se recomienda utilizar el cri-
terio (b), basado en el error relativo porcentual, pues es más robusto
ante comportamientos extremos de la sucesión.
3. Método de la Secante: una de las dificultades del método de Newton-
Raphson podrı́a consistir en el cálculo de la derivada f 0 . Un método
alternativo que sólo utiliza evaluaciones de la función f es el método
de la secante definido como

f (xn 1 )(xn 1 xn 2 )
xn = xn 1 ,
f (xn 1 ) f (xn 2 )
para n = 2, 3, ... Se invita al estudiante a practicar esta técnica para
algunos de los ejercicios ya planteados. En general, el método de la
secante converge un poco más lento que el método de Newton-Raphson.

2 Ejercicios en Contexto Docente

Calcule un cero para la función f (x) o resuelva la ecuación dada en el intervalo
indicado.
1. ex + 2 x
+ 2 cos x 6 = 0, en [1, 2].
2. ln(x 1) + cos(x 1) = 0, en [1.3, 2].
3. 2x cos 2x (x 2)2 = 0, en [2, 3], y [3, 4].

6
 4. (x 2)2 ln x = 0, en [1, 2], y [e, 4].

5. f (x) = x2 4x + 4 ln(x), en [1, 2] y [2, 4].

6. f (x) = 2x + 3 cos x ex , en [0, 1].

1
7. f (x) = x tan(x), en [0, ⇡/2].
1
8. f (x) = x 2x , en [0, 1].
x
9. f (x) = 2 + ex + 2 cos(x) 6, en [1, 3].

10. f (x) = (x3 + 4x2 + 3x + 5)/(2x3 9x2 + 18x 2), en [0, 4].

11. f (x) = x tan(x), en [1, 2].

12. f (x) = x2 4x sin(x) + (2 sin(x))2 .

13. En el intervalo [5.5; 6.5]:

f (x) = x8 36x7 + 546x6 4536x5 + 22449x4 67284x3

+ 118124x2 109584x + 40320.

14. Calcule una aproximación del valor de

435.000
1.564.000 = 1.000.000 · e + · (e 1).

15. La suma de dos números es 20. Si cada uno se agrega a su raiz cuadrada,
el producto de las dos sumas es 155,55. Determine los dos números.

7
3 Ejercicios en Contexto Profesional.
Todos los problemas enunciados a continuación deben ser resueltos utilizando
todas las técnicas vistas en cátedra y laboratorio, tales como, (i) método
gráfico, (ii) comando fzero de Matlab, (iii) método de Newton-Raphson (eje-
cutado en la lı́nea de comandos, ejecutado a partir de un programa .m, eje-
cutado con calculadora cientı́fica, etc...), (iv) otros algoritmos,...
Además, en todos los casos debe garantizar un ciero nivel de exactitud de la
aproximación calculada.

1. Un abrevadero de longitud L tiene una sección transversal en forma de

semicı́rculo con radio r. Cuando se llena de agua hasta una distancia
h de la parte superior, el volumen V de agua es
V 1 h
= ⇡r2 r2 · arcsin( ) h(r2 h2 )1/2 .
L 2 r
Suponga L = 10[pie], r = 1[pie], V = 12.4[pie3 ]. Determine la profun-
didad h del agua en el abrevadero hasta 0.01[pie].

2. Una partı́cula parte del reposo sobre un plano inclinado uniforme, cuyo
ángulo ✓ cambia con una rapidez constante de d✓ dt
= ! < 0. Al final de
t segundos, la posición del objeto está dada por

g et! e t!
x(t) = · ( sin(t!)).
2! 2 2
Suponga que la partı́cula se desplazó 1, 7[pie] en 1[s]. Encuentre la
rapidez ! con que ✓ cambia. Suponga g = 32, 17[pie/s2 ].

3. En estudios hechos sobre la recolección de energı́a solar enfocando un

campo de espejos planos sobre un colector central, deduce una ecuación
para el factor de concentración geométrica C (L.L.Vant-Hull,Solar En-
ergy,18,33 (1976)):

⇡(h/ cos A)2 F

C= ,
0.5⇡D2 (1 + sin A 0.5 cos A)

en donde A es el ángulo del borde del campo, F es la cobertura frac-

cional del campo con espejos, D es el diámetro del colector, y h es la
altura del colector. Encuentre A si h = 300, C = 1200, F = 0.8 y
D = 14.

8
4. Un objeto que cae verticalmente en el aire está sujeto a una resistencia
viscosa y también a la fuerza de gravedad. Suponga que dejamos caer
un objeto de masa m desde una altura h0 y que la altura del objeto
después de t[s] es
mg m2 g
h(t) = h0 t + 2 (1 e kt/m ),
k k
2
en donde, g = 32.17[pie/s ], h0 = 300[pie], m = 0.25[lb], k = 0.1[lb
s/pie]. Calcule el tiempo que tarda este peso en caer al suelo.
5. Un jugador A dejará en cero (por una puntuación de 21 a 0) al jugador
B en un partido de raquetbol con una probabilidad de
1+p p
P = ·( )21 ,
2 1 p + p2
en donde p denota la probabilidad de que A gane un intercambio de
tiros (independientemente del servicio). Determine el valor de p para
el cual P = 1/2. Interprete en términos del juego los valores de p y P .
6. En el diseño de vehı́culos para todo tipo de terreno, es necesario tener
en cuenta las fallas cuando se trata de librar dos tipos de obstáculos.
Una es la falla por rozamiento, y ocurre cuando el vehı́culo intenta
cruzar un obstáculo que hace que su fondo toque el suelo. La otra
recibe el nombre de falla por colisión de la defensa delantera y ocurre
cuando el vehı́culo desciende por una zanja y la defensa delantera toca
el suelo. El ángulo máximo ↵ que puede alcanzar un vehı́culo cuando
es el ángulo máximo en que no ocurre la falla por rozamiento satisface
la ecuación

A sin(↵) cos(↵) + B sin2 (↵) C cos(↵) E sin(↵) = 0,

en donde,
A = l sin( 1 ),
B = l cos( 1 ),
C = (h + 0.5D) sin( 1 ) 0.5D tan( 1 ),
E = (h + 0.5D) cos( 1 ) 0.5D.
(a) Se sabe que si l = 89[pulg], h = 49[pulg], D = 55[pulg], y 1 =
11.5 , el ángulo ↵ será aproximadamente igual a 33 . Se pide verificar
este resultado. Note que l es la distancia entre los centros de ambas
ruedas, h es la altura superior en que se apoya una de las ruedas, y
D/2 es el radio de las ruedas.
(b) Aproxime el valor de ↵ con D = 30[pulg], asumiendo los mismos
valores numéricos para los otros parámetros.

9
 7. Una relación para el factor de compresibilidad de los gases reales está
dada por (R.DeSantis,I.E.C.Fundamentals,15,184(1976))

1 + y + y2 y3
z= ,
(1 y)3

con y = b/4⌫, siendo b la corrección de van der Waals, y ⌫ el volumen

molar. Si z = 0.892, ¿cuál es el valor de y?.

8. El factor de fricción para un flujo de suspensión de partı́culas fibrosas

se ha relacionado con los números de Reynolds a través de

1 1 p 5.6
p = ln(RE · f ) + (14 ),
f k k
con f factor de fricción, RE número de Reynolds, y k constante deter-
minada por la concentración de la suspensión. Se sabe que para una
suspensión con 0.08% de concentración k = 0.28. Se pide calcular el
valor de f si RE = 3750 (Lee,Du↵y,A.I.Ch.E.Journal,Julio,1976 ).

9. En la ecuación de Redlich-Kwong se ha medido que P = 87.3; T =

486.9; ⌫ = 2.005; R = 1.98; y A(T = 486.9) = 0.0837. Calcular b:

R·T A(T )
P = .
⌫ b ⌫(⌫ b)
10. Las temperaturas en el interior de un material con fuentes de calor
incorporadas se determinan a partir de la relación:
p
e t/2
· cosh 1 (et/2 ) = Lcr /2.

Aproxime t si se sabe que Lcr = 0.88.

11. La velocidad ⌫ de caı́da de un paracaidista de masa m está dada por

g·m (c/m)·t
⌫= · (1 e ),
c
donde g = 9.8[m/s2 ], y el coeficiente de rozamiento está dado por
c = 14[kg/s]. Se pide determinar la masa m de tal modo que a los 7[s]
su velocidad sea igual a 35[m/s].

10
12. La concentración de saturación del oxı́geno disuelto en agua dulce puede
ser calculada a través de la siguiente relación

1.575701 · 105 6.642308 · 107

ln Osf = 139.34411 +
Ta Ta2
1.243800 · 1010 8.621949 · 1011
+ ,
Ta3 Ta4

en donde Osf es la concentración de saturación de oxı́geno disuelto en

agua dulce en 1atm(mg/L), y Ta es la temperatura absoluta en [K].
Se pide aproximar el valor de Ta para Osf = 10.
Para aguas naturales tı́picas con temperatura templada la ecuación an-
terior es válida para rangos de temperatura entre 0 C (en este extremo
Osf = 14.621[mg/L]) y 35 C (en este extremo Osf = 6.949[mg/L]).

13. Un balance de masa para un lago bien mezclado puede escribirse como:

dc p
V · =W Q·c k·V · c,
dt
con V = 106 [m3 ], Q = 105 [m3 /año], W = 106 [g/año], k = 0.2[ ]. Se
pide calcular el valor de la concentración c en estado estable ( dc
dt
= 0).
Se sugiere utilizar c = 4[g/m3 ] como valor inicial de las iteraciones.

11
4 Método de Newton Multivariado
Considere el problema de calcular P̃ = P̃ (x̃, ỹ), tal que, F (P̃ ) = 0, en donde,
F (x, y) = (f1 (x, y), f2 (x, y)), o sea, resolver el sistema no lineal de dos ecua-
ciones:

f1 (x, y) = 0
f2 (x, y) = 0.

Con el objeto de formular el método de Newton para este problema con-

sidere el desarrollo en serie de Taylor bivariada aplicado a las funciones f1 y
f2 :
@f1 @f1
0 ⇡ f1 (x + x, y + y) ⇡ f1 (x, y) + x + y
@x @y
@f2 @f2
0 ⇡ f2 (x + x, y + y) ⇡ f2 (x, y) + x + y ,
@x @y

siendo (x + x, y + y) un punto muy cercano a P̃ . Matricialmente esta

última aproximación se puede escribir como:
  
f1x (x, y) f1y (x, y) x f1 (x, y)
= ,
f2x (x, y) f2y (x, y) y f2 (x, y)
o sea,
J(P ) · P = F (P ),
0
en donde P = P (x, y), P = [ x, y] , y

f1x (x, y) f1y (x, y)
J(P ) = .
f2x (x, y) f2y (x, y)
A partir de lo anterior se define el método de Newton (bivariado) como la
sucesión de aproximación {Pn = Pn (xn , yn ); n = 0, 1, 2, 3, ...}, definida para
n = 0, 1, 2, 3, ..., como

Pn+1 = Pn + Pn ,
en donde, Pn es la solución del sistema de ecuaciones lineales

J(Pn ) · Pn = F (Pn ).

Note que lo anterior se extiende naturalmente a sistemas no lineales de 3

o más incógnitas.

12
 Ejemplo: aproxime la solución exacta P̃ del sistema no lineal de ecua-
ciones

f1 (x, y) = x2 + y 2 1 = 0
f2 (x, y) = y x2 = 0,
qp p
5 1 5 1
en donde, P̃ = (± 2
, 2
) ⇡ (±0.7861513777, 0.6180339887), uti-
lizando el método de Newton bivariado.

En este caso la sucesión de aproximación queda definida a través de

  
xn+1 xn xn
= + ,
yn+1 yn yn

con n = 0, 1, 2, 3..., en donde [ xn , yn ]0 es la solución del sistema de

ecuaciones lineales,
  
f1x (xn , yn ) f1y (xn , yn ) xn f1 (xn , yn )
· = ,
f2x (xn , yn ) f2y (xn , yn ) yn f2 (xn , yn )
o sea,
  
2xn 2yn xn x2n + yn2 1
· = .
2xn 1 yn y x2n
Note que el sistema anterior se puede resolver fácilmente utilizando la
regla de Gabriel Cramer.
Si se utiliza Po (1, 1) como vector inicial se obtienen los siguientes primeros
elementos de la sucesión de aproximación

n xn yn
0 1 1
1 0.83333333........... 0.66666666...........
2 0.78809523........... 0.61904761...........
3 0.78615406........... 0.61803444...........
4 0.78615137........... 0.61803398...........
6 0.786151377757423 0.618033988749895
Los cálculos anteriores son para aproximar la solución exacta del primer
cuadrante.

13
4.1 Ejercicios en Contexto Docente
Para cada uno de los sistemas de ecuaciones no lineales dados a continuación
se pide iterar hasta que ||Pn+1 Pn ||1 < 10 6 . En cada caso se sugiere un
valor inicial.

1.

3x2 y 2 = 0
3xy 2 x3 1 = 0
(x0 , y0 ) = (1, 1)

2.

ln(x2 + y 2 ) sin(xy) = ln 2 + ln ⇡
ex y + cos(xy) = 0
(x0 , y0 ) = (2, 2)

3.

x3 + x2 y xz + 6 = 0
ex + ey z = 0
y 2 2xz = 4
(x0 , y0 , z0 ) = ( 1, 2, 1)

4.

6x 2 cos(yz) 1 = 0
p
2
9y + x + sin(z) + 1.06 + 0.9 = 0
60z + 3e xy + 10⇡ 3 = 0
(x0 , y0 , z0 ) = (0, 0, 0)

14
4.2 Estudio de Casos: Los gases de invernadero y la
lluvia.
Fuente: Métodos numéricos para ingenieros, S.C.Chapra y R.P.Canale, 6ta.
Ed., McGrawHill, 2011.

Los niveles atmosféricos de diversos gases de invernadero han ido aumen-

tando durante los últimos 50 años. Además del calentamiento global, los
gases de invernadero también pueden influir en la quı́mica atmosférica . Una
pregunta que se puede hacer es cómo la tendencia de aumento en el dióxido
de carbono está afectando el pH del agua de lluvia. Fuera de las zonas ur-
banas e industriales, está bien documentado que el dióxido de carbono es el
principal determinante del pH de la lluvia. El pH es la medida de la activi-
dad de los iones de hidrógeno, y por lo tanto, de la acidez. Para soluciones
acuosas diluidas, se puede calcular como pH = log10 [H + ], en donde [H + ]
es la concentración molar de iones de hidrógeno.

El objetivo central de este estudio es: calcular el pH (en áreas

limpias siempre cae entre 2 y 12 ) del agua de lluvia a partir de [H + ], para
lo cual se dispone de datos de la presión parcial de CO2 en la atmosfera,
esto es, pCO2 , desde el año 1958 hasta el año 2003. Por ejemplo, en 1958:
pCO2 = 315[pmm], mientras que en 2003: pCO2 = 375[ppm]

El siguiente sistema no lineal de ecuaciones determina la quı́mica del agua:

[H + ][HCO3 ]
K1 = 106
KH pCO2
[H ][CO32 ]
+
K2 =
[HCO3 ]
Kw = [H + ][OH ]
KH pCO2
cT = + [HCO3 ] + [CO32 ]
106
0 = [HCO3 ] + 2[CO32 ] + [OH ] [H + ],

en donde KH = 10 1.46 es la constante de Henry, y K1 = 10 6.3

, K2 =
10 10.3 y Kw = 10 14 son los coeficientes de equilibrio.

15
 Las 5 incógnitas son:

1. cT : carbono inorgánico total.

2. [HCO3 ]: bicarbonato.
3. [CO32 ]: carbonato.
4. [H + ]: ion hidrógeno.
5. [OH ]: ion hidroxilo.

Note que pCO2 , en [ppm], denota la presión parcial de CO2 en la atmósfera.

Para datos recabados en Mauna Loa, Hawai, desde 1958 hasta 2003, se pudo
ajustar el polinomio:

pCO2 = 0.011852(t 1980.5)2 + 1.356975(t 1980.5) + 339,

en donde, t = 1953, 1954, ..., 2002, 2003.

A partir de la información dada se pide:

1. Demostrar que algebraicamente el sistema de 5 ecuaciones se puede

reducir a la siguiente ecuación en la incógnita [H + ]:

K1 K2 K 1 Kw
0= · KH · pCO2 + 2 6 · KH · pCO2 + + [H + ].
106 +
· [H ] +
10 · [H ] 2 [H ]

2. Resuelva la ecuación anterior para los años 1958 y 2003, utilizando el

método de Newton-Raphson, el método de la bisección y el comando
fzero de Matlab.
3. Utilice las aproximaciones obtenidas en el punto anterior para calcular
el pH del agua de lluvia en los años indicados.
4. Construya una tabla de valores en donde la primera columna sean los
años: 1960, 1970, 1980, 1990, y 2000. La segunda columna el valor de
pCO2 . La tercera columna el valor de [H + ]. La cuarta columna el valor
de pH. El método para calcular [H + ] puede elegirlo libremente.
5. Resuelva el sistema no lineal original de 5 ecuaciones para el año 2003
utilizando el método de Newton multivariado. Reporte el resultado de
al menos 10 iteraciones. Compare con los resultados obtenidos previa-
mente.

16
4.3 Estudio de Casos: Flujo Turbulento.
Fuente: Análisis Numérico, C.Gerald, Ed. Alfaomega, 1991.

Para el flujo turbulento de fluidos en una red interconectada la tasa de

flujo V de un nodo a otro es proporcional a la raı́z cuadrada de la diferencia
en la presión de los nodos. El problema consiste en calcular la presión en
cada nodo: pi , i = 1, 2, 3, 4, las cuales deben satisfacer el siguiente sistema de
ecuaciones no lineales:

p p p
0.3 500 p1 = 0.2 p1 p2 + 0.1 p1 p3
p p p
0.2 p1 p2 = 0.1 p2 p4 + 0.2 p2 p3
p p p
0.1 p1 p3 + 0.2 p2 p3 = 0.1 p3 p4
p p p
0.1 p2 p4 + 0.1 p3 p4 = 0.2 p4 0.

Note que los valores de b representan factores de conductancia en la

p
relación vij = bij pi pj .

17
4.4 Estudio de Casos: Biomatemática
Fuente: Análisis Numérico, R.L.Burden y J.D.Faires, 7ma. ed., Thomson
Learning, 2002.

El experimento biológico consiste en la determinación de la temperatura

máxima del agua XM en la que varias especies de hidra pueden sobrevivir sin
que su esperanza de vida disminuya. Una forma de resolver este problema
consiste en aplicar un ajuste ponderado de mı́nimos cuadrados de la forma
a
y = f (x) =
(x b)c

a un conjunto de datos experimentales. Los datos x de los datos se refieren

a la temperatura del agua. La constante b es la ası́ntota de la gráfica de f ,
y por tanto, es una aproximación a XM .

1. Demostrar que la elección de a, b, y c, para minimizar

n 
X 2
a
w i yi
i=1
(xi b)c

se reduce a resolver el sistema no lineal (todas las sumatorias son para

i = 1, 2, ..., n):

X 1 X
w i yi
a =
(xi b)2c (xi b)c
X w i yi X 1 X w i yi X 1
c+1
· = ·
(xi b) (xi b)2c (xi b)c (xi b)2c+1
X wi yi ln(xi b) X 1 X w i yi X ln(xi b)
c
· = ·
(xi b) (xi b)2c (xi b)c (xi b)2c

2. Con los siguientes datos, resuelva el sistema no lineal para las especies.
Utilice los pesos wi = ln yi :

i 1 2 3 4
yi 2.40 3.80 4.75 21.60
xi 31.8 31.5 31.2 30.2

18
4.5 Estudio de Casos: Tránsito Partı́culas Granulares
Fuente: Análisis Numérico, R.L.Burden y J.D.Faires, 9na. ed., Cengage
Learning, 2011, p.13.

Sea ✓ = (✓1 , ..., ✓N ). Considere el siguiente sistema de N ecuaciones no

lineales siendo ✓ el vector de variables independientes. Para n = 1, 2, ..., N 1,
sin ✓n+1 sin ✓n
fn (✓) = (1 µwn+1 ) (1 µwn ) = 0,
vn+1 vn
y
N
X
fN (✓) = y tan ✓i X = 0,
i=1

en donde, para n = 1, 2, ..., N,

n
X 1
vn2 = v02 + 2gn y 2µ y
j=1
cos ✓j
y
N
X 1
wn = yvn .
v3
i=1 i
cos ✓i
Además, g = 32.17pies/s2 , µ = 0, X = 2, y = 0.2, y v0 = 0.
Se pide:

1. Escribir el sistema de ecuaciones no lineal resultante para los casos N =

2, 3, 4. Debe reemplazar g = 32.17pies/s2 , µ = 0, X = 2, y = 0.2, y
v0 = 0.

2. Escriba la matriz jacobiana para N = 3.

3. Se pide ejecutar 10 iteraciones del método de Newton, esto es, cal-

cular hasta ✓10 , para el caso N = 5 con aproximación inicial ✓0 =
(1, 1, 1, 1, 1) + ✏, siendo

d1 d2 d3 d4 d5
✏ = (10 , 10 , 10 , 10 , 10 ),

en donde di , i = 1, 2, 3, 4, 5, son los 5 primeros dı́gitos de su RUN. Por

favor, si di = 0 ó di = 1 use di = 2.

19
5 ANEXO: ECUACIONES CON SOLUCIÓN
ANALÍTICA
Obtenga la solución exacta o analı́tica de las siguientes ecuaciones y sistemas.
Se recomienda verificar sus resultados utilizando una CAS (Computer Algebra
System) tal como MAXIMA (gratuito, http://maxima.sourceforge.net/ ) o
MAPLE (comercial).
1/3 1/x 11+3/x
1. 1/3+1/x
= 11 8/x
.
1/x+1 1/x 1
1/x 1 1/x+1 4
2. 1 x = 2+x
.
1+x

3. xx x x
= 3(1 + x x ).
4 18x2 +86
4. 6x = 7776.
0.5
5. 2x 5x+1 = 10 8
.

6. x2/n + 6 = 5x1/n .

7. 1 + 8x6/5 + 9x2/5 .
4
8. 9 + x = 10x 2 .
64
9. 4x+1 + 4x
= 257.
log(35 x3 )
10. log(5 x)
= 3.
p p
11. a x b x= pa b .
b x
q p q p
p p
12. x x 2x 2 x + x + x + 2 = 0.

13.

x/y = 3/2
2
x y 2 = 45.

14.

(x2 + y 2 )(x + y) = 1080

(x2 + y 2 )(x y) = 540.

20
15.

ax ay = a15
x y = 3.

16.

x = (y + a)/3 a/2
y = (x + b)/3 b/2.

17.
p p p
12
x
aya = a7
p
x
py
a3 : a4 = 1.

18.

x + y/a + z/b = 5a + 4b
x + y/b + z/a = 4a + 5b
y + z = 7ab.

19.

xy/(5x + 4y) = 6
yz/(3y + 2z) = 6
xz/(3x + 2z) = 8.

21
BIBLIOGRAFÍA
Los problemas planteados han sido tomados de los siguientes textos:

1. Análisis Numérico, R.L.Burden, J.D.Faires, México:International Thom-

son Editores, 2002.

2. Métodos Numéricos para Ingenieros, S.Chapra, R.Canale, McGraw

Hill, 3ra. ed., 1999.

3. Algebra, Curso de Matemáticas Elementales, Francisco Pröschle, Edición

Occidente, Chile, 2006.

22
TAREA 1

1. Debe resolver el Problema 6, Página 9, de este apunte:

2. En el diseño de vehı́culos para todo tipo de terreno, es necesario tener

en cuenta las fallas cuando se trata de librar dos tipos de obstáculos.
Una es la falla por rozamiento, y ocurre cuando el vehı́culo intenta
cruzar un obstáculo que hace que su fondo toque el suelo. La otra
recibe el nombre de falla por colisión de la defensa delantera y ocurre
cuando el vehı́culo desciende por una zanja y la defensa delantera toca
el suelo. El ángulo máximo ↵ que puede alcanzar un vehı́culo cuando
es el ángulo máximo en que no ocurre la falla por rozamiento satisface
la ecuación

A sin(↵) cos(↵) + B sin2 (↵) C cos(↵) E sin(↵) = 0,

en donde,
A = l sin( 1 ),
B = l cos( 1 ),
C = (h + 0.5D) sin( 1 ) 0.5D tan( 1 ),
E = (h + 0.5D) cos( 1 ) 0.5D.
(a) Se sabe que si l = 89[pulg], h = 49[pulg], D = 55[pulg], y 1 =
11.5 , el ángulo ↵ será aproximadamente igual a 33 . Se pide verificar
este resultado. Note que l es la distancia entre los centros de ambas
ruedas, h es la altura superior en que se apoya una de las ruedas, y
D/2 es el radio de las ruedas.
(b) Aproxime el valor de ↵ con D = 30[pulg], asumiendo los mismos
valores numéricos para los otros parámetros.

23
3. Debe utilizar el método de Newton-Raphson.

4. Debe utilizar Matlab u Octave.

5. La Tarea 1 es individual.

6. Plazo máximo de entrega: 23:59 hrs del Domingo 1 de Abril.

7. Debe enviar vı́a e-mail, o subir a Educa, el archivo en formato PDF,

completamente tipeado en un procesador de texto idóneo (Word, Latex,
etc...).

8. Su informe debe contener a lo menos:

- Portada.
- Enunciado.
- Gráfico de la función auxiliar utilizada.
- Estimación inicial: ↵0 . Debe argumentar cómo la obtuvo.
- Una tabla con los valores numéricos ↵0 , ↵1 , ↵2 ,...que fueron calculados
con Newton-Raphson.
- Una evaluación del nivel de exactitud (o error) de su respuesta final.
- Una conclusión final con sus respuestas, junto con una breve discusión.
- Bibliografı́a.
- Anexo: programas, códigos, y/o funciones en Matlab u Octave uti-
lizados.
- Otros que usted estime pertinentes o sean indicados por su profesor.

..............................................
M.J.Andrade, S.Berres, E.Cariaga, D.Sepúlveda

24