3TA Metodos Numericos - Solucionario
3TA Metodos Numericos - Solucionario
3TA Metodos Numericos - Solucionario
RAICES NO LINEALES
1. Determine las raíces reales de 𝑓(𝑥) = −0.4𝑥 2 + 2.2𝑥 + 4.7
a) Gráficamente
b) Empleando la fórmula cuadrática
c) Usando el método de bisección con tres iteraciones para determinar la raíz más grande. Emplee
como valores iniciales x1 = 5 y xu = 10. Calcule el error estimado εo y el error verdadero ε1 para cada
iteración.
f(x) = −0.4𝑥² + 2.2𝑥 + 4.7
Solución: 10
−𝑏 ± √𝑏 2 − 4𝑎𝑐
𝑥𝑖 = , 𝑎 = −0.4, 𝑏 = 2.2, 𝑐 = 4.7
2𝑎
c) Método Bisección
i Xo Xm X1 (Xm-Xo) f(Xo) f(Xm) f(Xo).f(Xm)
0 5 --- 10 --- --- --- ---
1 5 7.5 10 2.50000 5.70000 -1.30000 (-)
2 5 6.25 7.5 1.25000 5.70000 2.82500 (+)
3 6.25000 6.875 7.5 0.62500 2.82500 0.91880 (+)
-1000 f(x)
b) Método Bisección
i Xo Xm X1 (Xm-Xo) f(Xo) f(Xm) f(Xo).f(Xm) E%
0 0 --- 1 --- --- --- ---
1 0 0.5 1 0.50000 -2.00000 1.00000 (-) 100.00
2 0 0.25 0.5 0.25000 -2.00000 -0.46880 (+) 100.00
3 0.25 0.375 0.5 0.12500 -0.46880 0.23830 (-) 33.33
4 0.25 0.3125 0.375 0.06250 -0.46880 -0.11770 (+) 20.00
5 0.3125 0.34375 0.375 0.03125 -0.11770 0.05910 (-) 9.09
3. Determine las raíces reales de 𝑓(𝑥) = −26 + 82.3𝑥 − 88𝑥 2 + 45.4𝑥 3 − 9𝑥 4 + 0.65𝑥 5
a) Gráficamente
b) Usando el método de bisección para localizar la raíz más grande con εs = 10%. Utilice como
valores iniciales x1 = 0.5 y xo = 1.0
c) Realice el mismo cálculo en b), pero con el método de la falsa posición y εs = 0.1%
-16000
f(x)
-18000
b) Método Bisección
i Xo Xm X1 (Xm-Xo) f(Xo) f(Xm) f(Xo).f(Xm) E%
0 0.5 --- 1 --- --- --- ---
1 0.5 0.75 1 0.25000 -1.71720 2.68470 (-) 33.33
2 0.5 0.625 0.75 0.12500 -1.71720 0.83520 (-) 20.00
3 0.5 0.5625 0.625 0.06250 -1.71720 -0.33420 (+) 11.11
4 0.56250 0.59375 0.625 0.03125 -0.33420 0.27470 (-) 5.26
5 0.56250 0.57813 0.59375 0.01563 -0.33420 -0.02340 (+) 2.70
6 0.57813 0.58594 0.59375 0.00781 -0.02340 0.12720 (-) 1.33
7 0.57813 0.58203 0.58594 0.00391 -0.02340 0.05230 (-) 0.67
8 0.57813 0.58008 0.58203 0.00195 -0.02340 0.01460 (-) 0.34
9 0.57813 0.57910 0.58008 0.00098 -0.02340 -0.00440 (+) 0.17
10 0.57910 0.57959 0.58008 0.00049 -0.00440 0.00510 (-) 0.08
-1
b) Método de Bisección:
i Xo Xm X1 (Xm-Xo) f(Xo) f(Xm) f(Xo).f(Xm) E%
0 0.5 --- 2 --- --- --- ---
1 0.5 1.25 2 0.75000 -2.08630 -0.25370 (+) 60.00
2 1.25 1.625 2 0.37500 -0.25370 0.27100 (-) 23.08
3 1.25 1.43750 1.625 0.18750 -0.25370 0.02580 (-) 13.04
Por lo tanto la raíz viene a ser el valor de 1.43750
c) Método de la falsa posición:
i Xo X1 X2 (X2-Xo) f(X2) E%
0 0.5 2 --- --- ---
1 0.5 2 1.62871 1.12871 0.27558 69.30
2 1.62871 2 1.37959 0.24912 -0.05643 18.06
3 1.37959 2 1.42672 0.04713 0.01076 3.30
Solución:
a) Analíticamente:
0.9
0.9 = 0.4𝑥 → 𝑥=
0.4
0.4
b) Gráficamente, se realiza la tabulación
en base a la función dada: 0.2
x f(x) x f(x) 0
-6 -4 -2 0 2 4 6
-5 -0.58000 0 Error -0.2
-1 -1.30000 4 -0.17500 -1
-1.2
-1.4
8. Calcular la raíz cuadrada positiva de 15 usando el método de la falsa posición con εs = 0.5% Emplee
como valores iniciales xl = 3 y xa = 4
Solución:
𝑓(𝑥) = 𝑥 − √15
a) Método de la falsa posición:
i Xo X1 X2 (X2-Xo) f(X2) E%
0 3 4 --- --- ---
1 3 4 3.87298 0.87298 1.96799 22.54
2 3.87298 4 3.87298 0.00000 1.96799 0.00
9. Encuentre la raíz positiva más pequeña de la función (x está en radianes) 𝑥 2 |𝑆𝑒𝑛 𝑥̅ | = 5, usando el
método de la falsa posición. Para localizar el intervalo en donde se encuentra la raíz. Grafique
primero esta función para valores de x entre 0 y 5. Realice el cálculo hasta que εa sea menor que εs
= 1%. Compruebe su respuesta final sustituyéndola en la función original.
f(x) = 𝑥² |𝑆𝑒𝑛 𝑥 ̅ |-5
Solución:
20
10. Calcule la raíz real positiva de 𝑓(𝑥) = 𝑥 4 − 8𝑥 3 − 36𝑥 2 + 462𝑥 − 1010, utilizando el método de
la falsa posición. Use una gráfica para escoger el valor inicial y realice el cálculo con εs = 1.0%.
f(x) = 𝑥4 − 8𝑥3 − 36𝑥2 + 462𝑥 − 1010
Solución:
500
a) Gráficamente, se realiza la tabulación
en base a la función dada: 0
-3 -2 -1 0 1 2 3 4 5 6 7 8
x f(x) x f(x)
-500
-2 -1998 3 -83
-1 -1499 4 6
-1000
0 -1010 5 25
1 -591 6 34 -1500
2 -278 7 117
-2000
f(x)
-2500
b) Método de la falsa posición:
i Xo X1 X2 (X2-Xo) f(X2) E%
0 3 4 --- --- ---
1 3 4 3.93258 0.93258 5.99656 23.71
2 3.93258 4 3.87619 0.05639 5.07105 1.45
3 3.87619 4 3.88137 0.00518 5.15570 0.13
4 3.88137 4 3.88091 0.00046 5.14815 0.01
Solución:
a) Analíticamente:
3.3
𝑥 3.3 = 79 𝑥= √79 𝒙 = 𝟑. 𝟕𝟓𝟖𝟕
𝑔𝑚
12. La velocidad de caída de un paracaidista está dada por: 𝑣 = 𝑐 (1 − 𝑒 −(𝑐/𝑚)𝑡 )
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 v = 35 m/s en t = 7 s. Con el método de la falsa posición
determine “m” a un nivel de εs = 0.1%.
Solución:
9.8𝑚 14
−( )(7)
35 = (1 − 𝑒 𝑚 )
14
98
−( )
𝑓(𝑚) = 0.7𝑚 (1 − 𝑒 𝑚 ) − 35
x f(x) x f(x)
61 -0.86463 66 0.73400
62 -0.53355 67 1.03734
63 -0.20827 68 1.33549
64 0.11132 69 1.62859
65 0.42537 70 1.91674
Método de la falsa posición:
i Xo X1 X2 (X2-Xo) f(X2) E%
0 63 64 --- --- ---
1 63 64 63.65169 0.65169 0.00064 1.02
2 63.65169 64 63.64969 0.00200 0.00000 0.00
3 63.64969 64 63.64969 0.00001 0.00000 0.00
14. Con el método de iteración simple de punto fijo localice la raíz de 𝑓(𝑥) = 𝑆𝑒𝑛(√𝑥) − 𝑥. Use un
valor inicial de xo = 0.5 y haga iteraciones hasta que εs ≤ 0.01%.
Solución:
𝑔(𝑥) = 𝑆𝑒𝑛(√𝑥)
15. Utilice a) la iteración de punto fijo y b) el método de Newton-Raphson para determinar la raíz de
𝑓(𝑥) = −0.9𝑥 2 + 1.7𝑥 + 2.5 usando xo = 5. Efectué el cálculo hasta que εo sea menor que εs =
0.01%
También realice una prueba de error en su respuesta final.
Solución:
2.5 2.5
𝑔(𝑥) = 𝑔′ (𝑥) = 𝑔′ (𝑥) = 7.0268 𝑁𝑜 𝑐𝑢𝑚𝑝𝑙𝑒
0.9𝑥 − 1.7 (0.9𝑥 − 1.7)2
Resolviendo se obtiene el siguiente cuadro:
i Xo X1 (X1-Xo) Ea%
1 5 3.49603 1.50397 43.02
2 3.49603 3.06291 0.43312 14.14
3 3.06291 2.92631 0.13660 4.67
4 2.92631 2.88188 0.04442 1.54
5 2.88188 2.86729 0.01460 0.51
6 2.86729 2.86248 0.00481 0.17
7 2.86248 2.86089 0.00159 0.06
8 2.86089 2.86036 0.00052 0.02
9 2.86036 2.86019 0.00017 0.01
10 2.86019 2.86013 0.00006 0.00
b) Método de Newton-Raphson:
100
a) Gráficamente, se realiza la tabulación
en base a la función dada: 80
x f(x) x f(x)
60
-2 -34.00 4 -10.00
-1 -12.50 5 -9.50 40
0 -2.00 6 -2.00
20
1 0.50 7 15.50
2 -2.00 8 46.00 0
-2 -1 0 1 2 3 4 5 6 7 8 9
3 -6.50 9 92.50
-20
f(x)
-40
b) Método de Newton-Raphson:
17. Emplee el método de Newton-Raphson para determinar la raíz real de 𝑓(𝑥) = −2.0 + 6𝑥 − 4𝑥 2 +
0.5𝑥³ usando valores iniciales de a) 4.2 y b) 4.43. Discuta y use métodos gráficos y analíticos para
explicar las peculiaridades de los resultados.
Solución:
𝑓′(𝑥) = 6 − 8𝑥 + 1.5𝑥²
𝑓′(𝑥) = 6 − 8𝑥 + 1.5𝑥²
-2 -34.00 4 -10.00 80
-1 -12.50 5 -9.50
60
0 -2.00 6 -2.00
1 0.50 7 15.50 40
2 -2.00 8 46.00
20
3 -6.50 9 92.50
0
-2 -1 0 1 2 3 4 5 6 7 8 9
-20
En esta función existe una peculiaridad en
cuanto a las raíces reales, como se puede f(x)
-40
observar el intervalo [0,1], como el
intervalo [6,7], como se muestra en el,
gráfico y el cuadro, que hay cambio de signos, considerando los valores inicial de 4.2 y 4.43, las
raíces son 0.47457, la primera con 10 iteraciones y la segunda con 25 iteraciones.
18. Determine la menor raíz real de 𝑓(𝑥) = −11 − 22𝑥 + 17𝑥 2 − 25𝑥³: a) gráficamente y b) usando
el método de la secante para un valor de εs con tres cifras significativas.
300
a) Gráficamente, se realiza la tabulación f(x)
en base a la función dada: 250
1 -18.50 6 -71.00 0
-3 -2 -1 0 1 2 3 4 5 6
-50
Tomando los valores iniciales [-1,0]
-100
b) Método de la Secante:
i Xo X1 X2 Es%
1 -1 0 -0.26506 -
2 0 -0.26506 -0.41227 35.71
3 -0.26506 -0.41227 -0.37928 8.70
4 -0.41227 -0.37928 -0.38130 0.53
5 -0.37928 -0.38130 -0.38133 0.008
19. Localice la primera raíz positiva de: 𝑓(𝑥) = 𝑠𝑒𝑛 𝑥 + 𝐶𝑜𝑠(1 + 𝑥 2 ) − 1, donde x está en radianes.
Use cuatro iteraciones con el método de la secante con valores iniciales de
a) xi-1 = 1.0 y xi = 3.0,
b) xi-1 = 1.5 y xi = 2.5. Para localizar la raíz
c) Use el método grafico para verificar los resultados.
Solución:
0.00
c) Gráficamente, se realiza la tabulación -3 -1 1 3 5
en base a la función dada: -0.50
-1 -2.26 4 -2.03
-2.00
0 -0.46 5 -1.31
f(x)
1 -0.57 6 -0.51 -2.50
20. Calcule la raíz real de 𝑥 3.3 = 79, con el método de la secante modificado que cumpla con εs = 0.1%.
Intente diferentes valores de y analice los resultados.
Solución:
110
a) Gráficamente, se realiza la tabulación
en base a la función dada: 90
x f(x) x f(x) 70
-2 -60.1 3 -0.1 50
-1 -24.1 4 5.9 30
0 -6.1 5 23.9 10
-50 f(x)
-70
𝑓(𝑥𝑜 )
𝑥1 = 𝑥𝑜 −
𝑓′(𝑥𝑜 )
f(x)
-50.000
23. La función 𝑥 3 + 2𝑥 2 − 5𝑥 + 3 tiene una raíz doble en x = 1. Use a) el método estándar de Newton-
Raphson b) el método de Newton-Raphson modificado de la ecuación para resolver en la raíz x = 1.
Compare y analice la velocidad de convergencia usando xo = 0.2.
Solución:
𝑓(𝑥) = 𝑥 3 + 2𝑥 2 − 5𝑥 + 3
𝑓′(𝑥) = 3𝑥 2 + 4𝑥 − 5
𝑓(𝑥) = 𝑥 3 + 2𝑥 2 − 5𝑥 + 3
𝑓′(𝑥) = 3𝑥 2 + 4𝑥 − 5
𝑓′′(𝑥) = 6𝑥 + 4
26. El método de “divide y promedia”, es un antiguo método para aproximar la raíz cuadrada de
𝑥+𝑎/𝑥
cualquier número positivo a. se puede formular como 𝑥 = 2 . Demuestre que esta fórmula es
equivalente al algoritmo de Newton-Raphson para resolver 𝑥 = √𝑎.
Solución:
𝑥 = √𝑎
𝑓(𝑥) = 𝑥 2 − 𝑎
Derivando se obtiene:
𝑓 ′ (𝑥) = 2𝑥
𝑥𝑖 2 − 𝑎
𝑥𝑖+1 = 𝑥𝑖 −
2𝑥𝑖
𝒙𝒊 + 𝒂/𝒙𝒊
𝒙𝒊+𝟏 = 𝑳𝒒𝒒𝒅
𝟐
27. a) Aplique el método de Newton-Raphson a la función 𝑓(𝑥) = 𝑇𝑎𝑛ℎ(𝑥 2 − 9) para evaluar su raíz
real ya conocida en x = 3. Con valor inicial use xo = 3.1 y realice un mínimo de 4 iteraciones b) ¿El
método muestra convergencia hacia su raíz real? Dibuja una gráfica con los resultados para cada
iteración efectuada.
Solución:
𝑓(𝑥) = 𝑇𝑎𝑛ℎ(𝑥 2 − 9)
1.0
0.5
0.0
-11 -9 -7 -5 -3 -1 1 3 5 7 9 11
-0.5
-1.0
f(x)
-1.5
28. El polinomio 𝑓(𝑥) = 0.0074𝑥 4 − 0.28𝑥 3 + 3.355𝑥 2 − 12.183𝑥 + 5, tiene una raíz real entre 15 y
20.
Aplique a esta función el método de Newton Raphson usando como valor inicial xo = 16.15
Explique sus resultados.
Solución:
𝑓(𝑥𝑜 )
𝑥1 = 𝑥𝑜 −
𝑓′(𝑥𝑜 )
30. Los reactores de flujo tipo tapón (es decir, aquellos que en el fluido va de un extremo al otro con
una mezcla mínima a lo largo del eje longitudinal) se usan para convertir reactantes en producto.
Se ha determinado que la eficiencia de la conversión algunas veces se mejora recirculando una
porción de la corriente del producto, de tal forma que regrese a la entrada para un paso adicional
a través del reactor (ver figura inferior). La razón de recirculando se define como:
Suponga que se está procesando una sustancia química A para generar un producto B. Para el caso
en que A forma a B de acuerdo con una reacción auto catalítico (es decir, en la cual uno de los
productos actúa como catalizador o estimulante en la reacción). Es posible demostrar que una
razón óptima de recirculación debe satisfacer:
1 + R(1 − X Af ) R+1
Ln =
(1 − X Af ) R[1 + R(1 − X Af )]
Donde XAf es la fracción del reactante A que se convierte en el producto B. La razón óptima de
recirculación corresponde a un reactor de tamaño mínimo necesario para alcanzar el nivel deseado
de conversión. Utilice un método numérico para determinar la razón de recirculación necesaria, de
manera que se minimice el tamaño del reactor para una conversión fraccional de XAf = 0.9
Alimentación Producto
Reactor de flujo
Tipo Tapón
Reciclaje
Representación esquemática de un reactor de flujo tipo tapón con recirculación
Solución:
1 + R(0.1) R+1
Ln =
(0.1) R[1 + R(0.1)]
R+1
f(R) = − Ln(10 + R)
R + 0.1R2
Solución:
32. La siguiente ecuación permite calcular la concentración de un químico en un reactor donde se tiene
una mezcla completa: 𝐶 = 𝐶𝑒𝑛𝑡 (1 − 𝑒 −0.04𝑡 ) + 𝐶𝑜 𝑒 −0.04𝑡
Si la concentración inicial es Co = 4 y la concentración de entrada es Cent = 10. Calcule el tiempo
requerido para que C sea el 93% de Cent.
Solución:
Ce
K=
Ca 2 Cb
Cco + x
K=
(Cao − 2x)2 (Cbo − x)
Solución:
4+x 4+x
0.015 = f(x) = − 0.015
(42 − 2x)2 (30 − x) (42 − 2x)2 (30 − x)
39. La ecuación de Ergun, que se da abajo sirve para describir el flujo de un líquido a través de un lecho
empacado. P es caída de presión, es la densidad del fluido, Go es la velocidad másica (el cociente
del flujo de masa dividido entre el área de la sección transversal), Dp es el diámetro de las partículas
dentro del lecho, es la viscosidad del fluido, L es la longitud del lecho y es la fracción vacía del
lecho.
∆𝑃. 𝜌. 𝐷𝑝 . 𝜀 3 (1 − 𝜀)
2 = 150 + 1.75
𝐺𝑜 . 𝐿(1 − 𝜀) 𝐷𝑝 . 𝐺𝑜
( 𝜇 )
Dados los siguientes valores para los parámetros encuentre fracción vacía del lecho
𝐷𝑝. 𝐺𝑜 ∆𝑃. 𝜌. 𝐷𝑝
= 1000 = 10
𝜇 𝐷𝑜 2 𝐿
Solución:
10. 𝜀 3 (1 − 𝜀) 10. 𝜀 3
= 150 + 1.75 = 0.15(1 − 𝜀) + 1.75
(1 − 𝜀) (1000) (1 − 𝜀)
47. En ingeniería ambiental, la siguiente ecuación sirve para calcular el nivel de oxígeno en un rio aguas
abajo desde una descarga de aguas residuales:
Donde x es la distancia aguas abajo en kilómetros. Determine la distancia aguas abajo donde el nivel
de oxigeno se encuentra a una lectura de 5. (Sugerencia: este valor está dentro de los 2 km de la
descarga) Determine una respuesta con un 1% de error.
Solución:
48. La concentración de bacterias contaminantes “c” en un lago decrece de acuerdo con la relación
Solución:
49. En ingeniería marítima, la ecuación de una ola estacionaria reflejada en un puerto está dada por
2𝜋𝑥 2𝜋𝑡𝑤
𝜆 = 16, 𝑡 = 12, 𝑣 = 48, ℎ = ℎ𝑜 [𝑆𝑒𝑛 ( 𝜆
) 𝐶𝑜𝑠 ( 𝜆 ) + 𝑒 −𝑥 ] Encuentre el valor positivo más
bajo de x si: h = 0.4ho
Solución:
52. El valor acumulado de una cuenta de ahorros puede calcularse con la ecuación de anualidad vencida
𝑃[(1 + 𝑖)𝑛 − 1]
𝐴=
𝑖
En esta ecuación A es el monto de la cuenta, P es la cantidad que se deposita periódicamente e i es
la tasa de interés por periodo para los n periodos de depósito. A un ingeniero le gustaría tener una
cuenta de ahorros con un monto de $ 750,000 dólares al momento de retirarse dentro de 20 años,
y puede depositar $ 1,500 dólares mensuales para lograr dicho objetivo. ¿Cuál es la mínima tasa de
interés a que puede invertirse ese dinero, suponiendo que es un interés compuesto mensual?.
Solución: