UNIVERSIDAD MILITAR NUEVA GRANADA

FACULTAD DE ESTUDIOS A DISTANCIA

CARRERA INGENIERÍA CIVIL

MÉTODOS NUMÉRICOS

Elaborado por:

Mario Franco Morales

Bogotá, 13 de octubre de 2021.

2) Usando el método de solución de ecuaciones no lineales indicado para cada problema,
aproximar la respuesta usando las iteraciones indicadas. En cada ejercicio, debe incluir
las respuestas de cada una de las iteraciones.

a) La saturación de oxígeno disuelto en agua fresca puede ser calculada por la

ecuación
1,5757∗10 5 6,6423∗107 1,2438∗1010 8,6219∗10 11
ln ( os ) =−139,34411+ − + −
T +273,15 ( T + 273,15 )2 ( T +273,15 )3 ( T +273,15 )4

donde o s es la saturación de oxígeno disuelto en agua fresca a 1 atm ( mgL ), T es la

mg
temperatura en ° C . Teniendo en cuenta que a 0 ° C la saturación es 14,6 y a 35 ° C la
L

mg
saturación es 6,9 . Estimar la temperatura requerida para obtener una concentración
L

mg
de 12 usando el método de la bisección, con un error absoluto máximo de 0,05 ° C .
L

Solución.

Para facilitar los cálculos, si a la ecuación original denomina n=T +273,15 y se efectúan
operaciones, se tiene

−139,34411 n 4+ 1,5757∗10 5 n3 −6,6423∗10 7 n 2+1,2438∗1010 n−8,6219∗1011

ln ( 12 )=
n4

→ 141,829017 n4−1,5757∗105 n3 +6,6423∗107 n2−1,2438∗1010 n+ 8,6219∗1011 =0

Se define la función

f ( n )=141,829017 n4 −1,5757∗105 n3+ 6,6423∗107 n2 −1,2438∗1010 n+8,6219∗1011=0

Se deberá seleccionar un intervalo de confianza asignando valores a n

n f(n)
273,00 -1,130E+09
275,00 -8,441E+08
280,50 -2,661E+07
290,00 1,499E+09
308,00 4,884E+09
 Se deduce que la raíz del polinomio se encuentra en el intervalo ( 280,5 ; 2 90 ).

La primera aproximación es

( 273+308 )
nr = =290,50 → f ( n0 ) =1,5838∗10 9> 0
2

Luego, se tienen las siguientes aproximaciones

( 273+290,50 )
n1 = =281,75→ f ( n 1) =1,6576∗108 >0
2

( 273+281,75 )
n2 = =277,38→ f ( n 2) =−4,9688∗108 < 0
2

( 277,38+281,75 )
n3 = =279,56→ f ( n3 )=−1,6929∗108 <0
2

( 279,56+281,75 )
n 4= =280,66 → f ( n 4 )=−2,7021∗106 <0
2

Las cuales se consolidan en la siguiente tabla

Iteración na nb nr f(na) f(nr) f(na)*f(nr) ea

1 273,00 308,00 290,50 -1,1297E+09 1,5838E+09 -1,7893E+18 9,83
2 273,00 290,50 281,75 -1,1297E+09 1,6576E+08 -1,8726E+17 1,08
3 273,00 281,75 277,38 -1,1297E+09 -4,9688E+08 5,6132E+17 3,30
4 277,38 281,75 279,56 -4,9688E+08 -1,6929E+08 8,4117E+16 1,11
5 279,56 281,75 280,66 -1,6929E+08 -2,7021E+06 4,5744E+14 0,02

El error absoluto viene dado por la siguiente fórmula

e a=|nreal −ncalculado|

Donde el valor real se obtuvo mediante el uso de la calculadora científica, resultando

nreal =280,6738923

Deberá cumplirse que |n real−ncalculado|< ε→|280,6738923−ncalculado|<0,05

El valor que cumple con la condición es n calculado=280,66 K=7,51 ° C

Cuyo error absoluto resulta

0,02<0,05

b) La ecuación de estado de Redlich-Kwong es una ecuación algebraica empírica

que relaciona temperatura, presión y volumen de los gases ideales

RT a
p= −
v−b v ( v+ b ) √ T

kJ
donde R es la constante universal de los gases 0,518 ( kg K )
, T es la temperatura

absoluta en Kelvin, p es la presión absoluta (kPa ), y v es el volumen de un kilogramo de

m3
gas ( )
kg
.

Los parámetros a y b son 12,56 y 1,86∗10−3, respectivamente. Estimar el volumen v para

una presión p=65000 kPa a una temperatura de T =−40° C , usando el método de

Newton

tomando un valor inicial de v 0=3 m 3 con un error de aproximación de 0,1.

Solución.

Se efectúan operaciones con la ecuación original

RT a RT √ T v ( v+ b )−a ( v−b ) RT √ T v 2+ RT √ T b v−av + ab

p= − = =
v−b v ( v+ b ) √ T ( v−b ) v ( v +b ) √ T √ T v ( v 2−b2 )

→ p √ T v ( v 2−b 2) =RT √ T v 2 + ( RT √ T b−a ) v + ab

→ p √ T v ( v 2−b 2) −RT √ T v 2−( RT √ T b−a ) v−ab=0

→ p √ T v 3− p √ T b2 v−RT √ T v 2−( RT √T b−a ) v−ab=0

Agrupando términos semejantes resulta

p √ T v 3−RT √ T v 2−( RT √ T b+ p √ T b2 −a ) v −ab=0

Dividiendo todos los términos entre p √ T queda

 RT 2 RT b 2 a ab
v3 −
p
v−
p( +b −
p √T
v−
p √T )
=0( i)

Se denominan los coeficientes

RT RT b 2 a ab
A= ; B= +b − ; C=
p p p √T p √T

Previo a la sustitución de valores, se debe tener presente que la temperatura debe

kJ J
convertirse a K , la presión de kPa a Pa y la constante R deberá llevarse de a
kg K kg K

T =−40° C=(−40+ 273,15 ) K =233,15 K

103 Pa
p=65000 kPa= 65000 kPa
( )
1 kPa (
=65∗106 Pa )
R=0 ,518
kJ
=0 , 518
( 1000
k J∗
J
1 kJ )
=518
J
kg K kg K kg K

Luego, se sustituyen los valores para obtener dichos coeficientes

J
A=
( 518
kg K )
∗(233,15 K)
=0,001858 03
m 3

6
65∗10 Pa kg

J
B=
( 518
kg K )
( 233,15 K )∗( 1,86∗10 −3
)
2
+ ( 1,86∗10−3 ) −
12,56
6 6
65∗10 Pa 65∗10 Pa √233,15 K

→ B=6,903∗10−6

12,56∗1,86∗10−3 −11
C= 6
=2,35∗10
65∗10 √ 233,15

Resultando la ecuación (i) como sigue

v3 −A v 2−B v −C=0

→ v 3−0,00185803 v 2−6,903∗10−6 v−2,35∗10−11 =0

Se define la función

f ( v )=v 3−0,00185803 v 2−6,903∗10−6 v−2,35∗10−11 =0

Se procede a hallar la primera derivada de f como sigue

3 2 ' 2
f ' (v )=( v − A v −B v−C ) =3 v −2 Av−B

→ f '( v )=3 v 2−2 ( 0,00185803 ) v−6,903∗10−6

→ f '( v )=3 v 2−0,00371605 v−6,903∗10−6

Se aplica el Método de Newton a partir de la ecuación

f ( vn )
v n+1=v n−
f ' ( vn )

Cuyo error deberá cumplir la condición

|v n+1−v n|<ε
Siendo la tolerancia equivalente a ε =0,1

Se tiene |v n+1−v n|<0,1

Tomando v 0=3, la primera iteración para n=0 resulta

f (v0 ) 26,983
v1 =v 0− '
=3− =2
f ( v0 ) 26,989

Cuyo error resulta: |2−3|=1> 0,1 No satisface

Segunda iteración, n=1

f (v1 ) 7 , 9 95
v 2=v 1− '
=2− =1,334
f ( v 1) 11,995

Cuyo error resulta: |1,334−2|=0,667>0,1 No satisface

Tercera iteración, n=2

 f (v2 ) 2 , 369
v3 =v 2− '
=1,334− =0,889
f ( v2) 5,331

Cuyo error resulta: |0,889−1,334|=0 , 444 >0,1 No satisface

Cuarta iteración, n=3

f ( v3 ) 0 ,702
v 4=v 3− '
=0,889− =0,593
f (v3 ) 2,369

Cuyo error resulta: |0,593−0,889|=0 , 296>0,1 No satisface

Quinta iteración, n=4

f ( v4 ) 0 , 208
v5 =v 4 − '
=0,593− =0,396
f (v4 ) 1,053

Cuyo error resulta: |0,396−0,593|=0 , 197>0,1 No satisface

Sexta iteración, n=5

f ( v5 ) 0 ,062
v 6=v 5 − '
=0,396− =0,264
f ( v5 ) 0,468

Cuyo error resulta: |0,264−0,396|=0,1 32>0,1 No satisface

Séptima iteración, n=6

f ( v6 ) 0,0 18
v7 =v 6 − '
=0,264− =0,176
f ( v6 ) 0 , 208

Cuyo error resulta: |0,176−0,264|=0 , 088< 0,1Cumpl e

A continuación, se muestra un cuadro resumen de resultados para mejor visualización

vi f(v) f'(v) vf e
3,000 26,983 26,989 2,000  
2,000 7,995 11,995 1,334 1,000
1,334 2,369 5,331 0,889 0,667
0,889 0,702 2,369 0,593 0,444
0,593 0,208 1,053 0,396 0,296
0,396 0,062 0,468 0,264 0,197
 0,264 0,018 0,208 0,176 0,132
0,176       0,088

m3
Finalmente, el volumen estimado será de 0,176 .
kg

c) La fórmula de Manning, fue propuesta por el ingeniero Robert Manning en 1889

para describir la velocidad del agua en canales abiertos y tuberías. En el caso de un canal
abierto rectangular se puede escribir como:

√ m ( Bh )5 /3
v= 2 /3
n ( B+2 h )

donde n es el coeficiente de rugosidad de la pared, m es la pendiente del muro, h la

profundidad y v el flujo. Estime la profundidad a la cual el agua fluye a una velocidad de

m3
5 , si la pared tiene un coeficiente de rugosidad de 0,03 y la pendiente m=0,0025 , y la
s
constante B=20 con un error relativo máximo del 0,5 % usando el método del punto fijo
tomando como lado inicial una profundidad de h0 =1 m.
Solución:

Para aplicar el método, es necesario efectuar un arreglo a la ecuación original de modo

que

h=g(h)

Despejando la h del numerador en la ecuación original, se tiene

( vn )3 /5 ( B+ 2h )2 /5
h=
B m 3/ 10

Sabiendo que la mayoría de los ríos son de ancho mayor que de profundidad, el valor de
B+2 h no varía mucho y se aproxima a B, por contrario el valor de B h es directamente
proporcional a h , entonces la ecuación anterior en la que se despeja h del numerador
convergerá con mayor velocidad que si se despejara del denominador de la ecuación
original.
Luego, se parte del valor h0 =1 para hallar h1:

( 5∗0,03 )3 /5 ( 20+2∗1 )2 /5
h1 = 3 /10
=0,3328
20∗( 0,0025 )

Se procede a calcular el siguiente valor de h

( 5∗0,03 )3 /5 ( 20+2∗0,3328 )2/ 5

h2 = 3 /10
=0,3246
20∗( 0,0025 )

El error será entonces

hi +1−hi hi +1−hi
e= | hi+1 |∗ 100 %< ε → e= |
h i+1 |
∗100 % <0,5 %

Para el caso en cuestión se tiene

e= |0,3246−0,3328
0,3246 |∗100 %=2,53 % >0,5 % No cumple
Se itera nuevamente

( 5∗0,03 )3 /5 ( 20+2∗0,3246 )2/ 5

h3 = 3/ 10
=0,324 5
20∗( 0,0025 )

Cuyo error resulta

e= |0,3245−0,3246
0,3245 |∗100 %=0 ,0 3 % <0,5 % C umple
Los resultados se consolidan en la siguiente tabla anexa:

i hi g(hi) = hi+1 e
1 1,0000 0,3328  
2 0,3328 0,3246 2,53%
3 0,3246 0,3245 0,03%

Finalmente, la profundidad estimada será de 0,3245 m.

 d) Se desea vender un vehículo por un precio de 69.000 .000 anunciando que, sólo
se deberá pagar 13.500 .000 anualmente por 6 años. Teniendo en cuenta que la fórmula
para calcular el pago anual del vehículo es:

i ( 1+ i )n
A=P n
( 1+i ) −1

donde A es el pago anual, P es el precio del vehículo, n el número de años e i es la tasa

de interés. Estimar el valor de la tasa de interés, i , usando el método de la Regula-Falsi
hasta alcanzar un error absoluto de 0,00001 si, la respuesta debe estar entre 0,01 y 0,3.

Solución.

Se efectúa el cambio de variable u=1+ i, de manera que i=u−1, donde la ecuación

original resulta

( u−1 )∗un
A=P n
→ A ( un−1 ) =P ( u−1 )∗un → A un− A−P un+1 + P un=0
u −1

→−P un +1+ ( A+ P ) u n− A=0

Donde

n=6 ; P=69000000 ; A=13500000

Sustituyendo los valores en la ecuación anterior, resulta

−69000000 u7 +82500000 u6−13500000=0

Dividiendo entre 100 0 000 cada término, se tiene

−69 u7 +82 , 5u 6−13 , 5=0

Se define la función

f ( u )=−69 u7 +82,5 u6 −13 ,5=0

Se deberá seleccionar un intervalo de confianza asignando valores a n

u f(u)
0,500000 -12,7500
0,750000 -8,0272
1,000000 0,0000
 1,250000 -27,8051
1,500000 -252,7031

Se deduce que la raíz del polinomio se encuentra en el intervalo ( 0 , 5 ; 1,5 0 ) y su valor

exacto es 1,00. Sin embargo, se procede a efectuar las operaciones del método de Falsa
Posición.

Asumiendo un intervalo más cerrado como (0,95 ; 1,01) Se halla la primera aproximación,
como sigue

f ( b ) ( b−a ) ( 0,098073 )( 1 , 01−0 ,9 5 )

u=b− =1 ,01− =1,004830
f ( b ) −f ( a ) ( 0,098073 )−(−1,040192 )

Luego, los resultados del cálculo se muestran en la tabla anexa:

Iteració
a b r f(a) f(b) f(a)*f(b) e
n
1 0,950000 1,010000 1,004830 -1,040192 0,098073 -0,102015  
2 0,950000 1,004830 1,002175 -1,040192 0,052943 -0,055071 0,00483
3 0,950000 1,002175 1,000946 -1,040192 0,025090 -0,026098 0,00217
4 0,950000 1,000946 1,000405 -1,040192 0,011162 -0,011611 0,00095
5 0,950000 1,000405 1,000172 -1,040192 0,004827 -0,005021 0,00041
6 0,950000 1,000172 1,000073 -1,040192 0,002061 -0,002144 0,00017
7 0,950000 1,000073 1,000031 -1,040192 0,000876 -0,000911 0,00007
8 0,950000 1,000031 1,000013 -1,040192 0,000371 -0,000386 0,00003
9 0,950000 1,000013 1,000006 -1,040192 0,000157 -0,000164 0,000013
10 0,950000 1,000006 -1,040192 0,000067 -0,000069 0,000006

El error absoluto fue obtenido como se ilustra

e=|u real−u i|<ε → e=|1−u i|< 0 , 00001

Para el caso en cuestión se tiene

e=|1−1,000006|=0,000006<0,00001 C umple

La raíz será u=1,000006, por lo que se devuelve el cambio i=u−1

→ i=1,000006−1=0 , 000006

El cual, expresado en porcentaje, resulta

i=0 , 0006 %