5.

9) Encuentre la raíz positiva más pequeña de la función (x está en radianes)

x 2|cos √ x =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 ea sea menor que error= 1%. Compruebe su respuesta final
sustituyéndola n la función original.

X F(X)
0 -5
1 -4.00015
2 -1.00121
3 3.99589
4 10.99024
5 19.98096
Metodo de la falsa posicion
Xi 1=a0 −¿

a 0=2 F ( a0 ) =−1.00121b 0=3 F ( b 0 )=3,99589

Aplicamos la ecuacion del metodo de falsa posicion de la siguiente manera:

1 interacion:
I −1.00121× ( 3−2 )
1=2−
[ 3.99589− ( −1.00121 ) ]
I 1=2.20036
Se evaluada 3.6819 en la ecuacion dada por el ejercicio.

x 2|cos √ x|−5=0

2.2004 2|cos √ 2.2004|−5=−0.15986> error dado por el ejercicio se vueleve hacer el mismo procedimiento hasta
a0 b0 F (a0 ) F (b0 ) Xi 1 F ( Xi 1 )=¿
1iteración 2 3 -1.00121 3.99589 2.2004 -0.15986

2iteración 2.2004 3 -0.15986 3.99589 2.2312 -0.02344

3iteración 2.2312 3 -0.02344 3.99589 2.2357 -0.00335

Se deja de iterrar ya que 0.00381 es menor que el error pedido por el ejercicio, esto
quiere decir que la solucion de raiz aproximadamente esta en 2.23567.

5.13 La velocidad u de un paracaidista que cae está dada por

gm
v= ¿
c
donde g = 9.81 m/s2. Para un paracaidista con coeficiente de resistencia de c = 15 kg/s,
calcule la masa m de modo que la velocidad sea y = 36 m/s en t = 10 s. Utilice el método
de la falsa posición para determinar m a un nivel de es = 0.1%.
solucion del problema 5.13
scilab
 Metodo de falsa posicion.

X F(X)
0 -36
10 -29.460
20 -22.927
30 -16.512
40 -10.455
50 -4.928
60 0.01898
70 4.409

Xi 1=a0 −¿
1 primera iteracion.

a 0=50 F ( a 0) =−4.928b 0=60 F ( b 0 )=0.01898

 −4.928× ( 60−50 )
Xi 1=50−
[ 0.01898−(−4.928 ) ]
=60.093

F ( Xi 1 )=0.0623
a0 b0 F (a0 ) F (b0 ) Xi 1 F ( Xi 1 )
1iteración 50 60 -4.928 0.01898 60.093 0.0623

2iteración 50 60.093 -4.928 0.0623 59.967 0.0036

3iteración 50 59.967 -4.928 0.0036 59.959 0.000130

5.15 Como se ilustra en la figura P5.15, la velocidad del agua, y (m/s), en la descarga de
un tanque cilíndrico a través de un tubo largo se puede calcular como

2 gH
v=√ 2 gH × tanh( √ ×t )
2L
donde g = 9.81 m/s2, H = carga hidrostática inicial (m), L = longitud de tubo (m) y t =
tiempo transcurrido (s). Determine la carga hidrostática necesaria para obtener u = 5 m/s
en 2.5 s para un tubo de 4 m de longitud a) gráficamente, b) por bisección y c) con
posición falsa. Utilice los valores iniciales de x1 = 0 y xu= 2 m, con un criterio de detención
de s = 1%. Revise sus resultados.

Solucion del problema 5.15

Metodo de biseccion.
C bn +an
n=an+
2

1 iteracion.
a o=0 y b =2o

C 2+0
0=0+ =1
2

F ( C 0 )=¿ √2 ×9.81 ×1 × tanh ( √ 2× 29.81×

L
1
×2.5 )−5

F ( C 0 )=−1.0937<0 esto quiere decir que este valor no es lacarga hidraulica .

2 iteracion.
a 1=1 y b =2o

C 2+1
1=1 + =1.5
2

F ( C 1 )= √ 2 ×9.81 ×1.5 × tanh ( √ 2× 9.81×

2L
1.5
×2.5 )−5

F ( C 1 )=0.071

El valor de carga hidraulica del tanque mostrado en el problema 5.13 es

aproximadamente de 1.5 metro por el metodo de la biseccion.
Metodo de falsa posicion.

X F(x)
0 -5
1 -1.0937
2 1.0192
3 -2.5461

Xi 1=a0 −¿
1 primera iteracion.

a 0=1 F ( a0 ) =−1.0937b 0=2 F ( b0 ) =1.0192

−1.0937 × ( 2−1 )
Xi 1=1−
[ 1.0192−(−1.0192 )] =1.517

F ( Xi 1 )=√ 2× 9.81× 1.517 × tanh ( √2 × 9.81×1

2L
.517
× 2.5 )−5

F ( Xi 1 )=0.1065
 a0 b0 F (a0 ) F (b0 ) Xi 1 F ( Xi 1 )
1iteració 1 2 -1.0937 1.0192 1.517 0.1065
n
2iteració 1.517 2 0.1065 1.0192 1.460 -0.0103
n
3iteració 1.460 2 -0.0103 1.0192 1.465 -0.0019
n

5.17 Suponga el lector que está diseñando un tanque esférico (véase la figura P5.17) para
almacenar agua para un poblado pequeño en un país en desarrollo. El volumen de líquido
que puede contener se calcula con

[ 3 R−h ]
V =π h2 ×
3
donde V = volumen [m3], h = profundidad del agua en el tanque [m] y R = radio del tanque
[m]. Si R = 3 m, ¿a qué profundidad debe llenarse el tanque de modo que contenga 30
m3? Haga tres iteraciones con el método de la falsa posición a fin de obtener la
respuesta. Determine el error relativo aproximado después de cada iteración. Utilice
valores iniciales de 0 y R.

Solucion 5.17
Metodo de falsa posicion.

X F(X)
0 -30
1 -21.622
2 -0.678
3 26.549
4 53.776

Xi 1=a0 −¿
1 primera iteracion.

a 0=2 F ( a0 ) =−0.678b 0=3 F ( b 0 )=26.549

−0.678 × (3−2 )
Xi 1=2−
[ 26.549−(−0.678 ) ]
=2.025

[ 3× 3−2.025 ]
F ( Xi 1 )=π ( 2.025 )2 × −30
3
F ( Xi 1 )=−0.0482
a0 b0 F (a0 ) F (b0 ) Xi 1 F ( Xi 1 )
1iteració 2 3 -0.678 26.549 2.025 -0.0482
n
2iteració 2.025 3 -0.0482 26.549 2.027 0.00239
n
3iteració 2.027 3 0.00239 26.549 2.027 0.00239
n

5.18 La concentración de saturación de oxígeno disuelto en aguadulce se calcula con la

ecuación (APHA, 1992)

1.575701 ×105 6.642308× 107 1.243800 ×1010 8.621949× 1011

ln o sf =−139.34411+ − + −
Ta T 2a T 3a T 4a

donde O sf = concentración de saturación de oxígeno disuelto en agua dulce a 1 atm

(mg/L) y Ta = temperatura absoluta (K). Recuerde el lector que Ta = T + 273.15, donde T
= temperatura (°C). De acuerdo con esta ecuación, la saturación disminuye con el
incremento de la temperatura. Para aguas naturales comunes en climas templados, la
ecuación se usa para determinar que la concentración de oxígeno varía de 14.621 mg/L a
0°C a 6.413 mg/L a 40°C. Dado un valor deconcentración de oxígeno, puede emplearse
esta fórmula y el métodode bisección para resolver para la termperatura en °C.
a) Si los valores iniciales son de 0 y 40°C, con el método de la bisección, ¿cuántas
iteraciones se requerirían para determinarla temperatura con un error absoluto de
0.05°C?
Se necesita hacer 3 iteraciones para poder llegar a una temperatura con error absoluto de
0.05°C.2

Temperatur Osf
a
0°C 14.621
40°C 6.413
C bn +an
n=an+
2

1 iteracion.
a o=0 y b =40
o

C 40+ 0
0=0+ =20
2

T a=20+273.15=293.15
Iteramos entre las temperatura dada el 20°C y asi hallar el Osf, para 20°C tenemos que
Osf=10.517 y ln(10.417)= 2.353.

1.575701 ×105 6.642308 ×107 1.243800 ×1010 8.621949 ×1011

0 ¿−139.34411+ − + − −ln ( 10.417 )=¿
293.15 293.152a 293.153a 293.154a
¿−0.146
ao bo I n+ 273.15 osf Ln(osf) F( I n)

1 iteracion 0 40 20 10.517 2.352 -0.146

2 iteracion 20 40 23 9.9014 2.2926 -0.143
3 iteracion 23 40 54.5 4.756 1.5595 0.0452
b) Desarrolle y pruebe un programa para el método de bisección a fin de determinar
T como función de una concentración dada de oxígeno, con un error absoluto
preespecificado como en el inciso a). Dadas elecciones iniciales de 0 y 40°C,
pruebe su programa para un error absoluto de 0.05°C para los casos siguientes:
osf = 8, 10 y 12 mg/L. Compruebe sus resultados.
Grafica.
Este resultado es para un osf de 8.10 mg/L.
Este resultado es para un osf de 12 mg/L.
5.19 De acuerdo con el principio de Arquímedes, la fuerza de flotación es igual al peso de
fluido desplazado por la porción sumergida de un objeto. Para la esfera ilustrada en la
figura P5.19, use la bisección para determinar la altura h de la porción que queda encima
del agua. Utilice los siguientes valores para su cálculo: r = 1 m, densidad de la esfera =
200 kg/m3 y densidad del agua = 1 000kg/m3. Observe que el volumen de la porción de la
esfera por encima del agua se puede calcular mediante

π h2
V= ( 3r −h)
3

Solucion 5.19
5.20 Realice el mismo cálculo que en el problema 5.19, pero para el frústum de un cono,
como se ilustra en la figura P5.20. Use los siguientes valores para su cálculo: r1 = 0.5 m,
r2 = 1 m, h = 1 m, densidad del frústum = 200 kg/m3, y densidad del agua = 1 000 kg/m3.
Observe que el volumen del frústum está dado por
 πh 2 2
V= (r + r + r × r )
3 1 2 1 2

6.19 Suponga el lector que está diseñando un tanque esférico (véase la figura P6.19) de
almacenamiento de agua para un poblado pequeño de un país en desarrollo. El volumen
del líquido que puede contener se calcula con

[ 3 R−h ]
V =π h2 ×
3
Donde V = volumen [m3], h = profundidad del agua en el tanque [m] y R = radio del
tanque [m]. Si R = 3 m, ¿a qué profundidad debe llenarse el tanque de modo que
contenga 30 m3? Haga tres iteraciones del método de Newton-Raphson para determinar
la respuesta. Encuentre el error relativo aproximado después de cada iteración. Observe
que el valor inicial de R convergerá siempre.

Solución 6.19

[ 3 R−h ]
V =π h2 ×
3
Donde:
R=3 ,h=? ,V =30 m3

9 × π ×h 2−π × h3−90=0
 Aplicamos el metodo dado por el problema que es por el metodo de Newton-
Raphson.

f ( y cn )
y c n = y c n−
f '( y cn )

f ( y c n ) =9 × π ×h2 −π × h3−90

f ' ( y c n )=18 × π × h−3 × π ×h 2

1 primera iteracion.
Supongamos en el h va de 0 a 3
0+3
y n= =1.5 siendo1.5 nuestro primero numero para iterar .
2
(−36.986)
n :0 → y c 0 =1.5−
53.014
y c 0=2.196

f ( y c 0)=π ¿

f ( y c 0 ) =4.361

2 iteracion.

(13.081)
n :1 → y c 0=2.196−
78.731
y c 0=2.029

f ( y c 1)=π ¿

f ( y c 1 )=0.0783

3 iteracion.

(0.235)
n :2 → y c 0=2.029−
75.955
y c 0=2.027

f ( y c 2)=π ¿

f ( y c 2 ) =0.00239
 Aproximadamente profundidad del agua en el tanque es de 2.207 con un error de
0.00239.
6.20 La ecuación de Manning se puede escribir para un flujo en un canal abierto
rectangular como:
5
√ S ( BH ) 3
Q= 2
3
n(B+2 H )
Donde Q = flujo [m3/s], S = pendiente [m/m], H = profundidad [m] y n = coeficiente de
rugosidad de Manning. Desarrolle un esquema de iteración para despejar H de esta
ecuación dado Q = 5, S = 0.0002, B = 20 y n = 0.03. Pruebe que su esquema converge
para todos los valores iniciales mayores que cero o iguales a cero.
Solución
5
3
√ 0.0002 ( 20 × H )
2
−5=0
3
0.003(20+2 H)
5 2
3 3
√ 0.0002 (20 × H ) −5 ×0.003 ( 20+2 H ) =0
Se toma un valor de H inicial de 0.68.
5 2
3 3
F ( H )= √0.0002 ( 20 × H ) −5 × 0.003 ( 20+ 2 H )

5 2
3 3
F ( X o )= √0.0002 ( 20 ×0.68 ) −5× 0.003 ( 20+2× 0.68 )

F ( X 0 )=−2.310.

5 2
F ( X 1 )=√ 0.0002 ( 20× 0.72 ) 3 −5 ×0.003 ( 20+2 ×0.72 ) 3

F ( X 0 )=−2.654
Método de la secante.
1 iteración.

f ( b0 ) × ( b0−a0 )
C 1=b0 −
f ¿¿
−2.654 × ( 0.68−0.72 )
C 1=0.72−
−2.310+ 2.654
C 1=0.411
5 2
F ( C 1 )= √ 0.0002 ( 20 ×0.411 ) 3 −5× 0.003 ( 20+2× 0.411 ) 3

F ( C 1 )=0.537

2 iteración.

−0.537× ( 0.72−0.411 )
C 1=0.411−
−2.652−(−0.537)
C 1=0.332
5 2
F ( C 1 )= √ 0.0002 ( 20 ×0.332 ) 3 −5 ×0.003 ( 20+2 ×0.332 ) 3

F ( C 1 )=−0.050

3 iteración.

−0.048 × ( 0.332−0.323 )
C 1=0.323−
−0.050−(−0.048)
C 1=0.322
5 2
3 3
F ( C 1 )= √ 0.0002 ( 20 ×0.322 ) −5 ×0.003 ( 20+2 ×0.322 )

F ( C 1 )=−0.048

6.23 Encuentre las raíces de las ecuaciones simultáneas que siguen:

x 2+ y 2=16
Use un enfoque gráfico para obtener los valores iniciales. Encuentre estimaciones
refinadas con el método de Newton-Raphson para dos ecuaciones, que se describe en la
sección 6.6.2.

6.25 El balance de masa de un contaminante en un lago bien mezclado se expresa así:

dc
V =W −Qc−kV √ c
dt
Dados los valores de parámetros V = 1 × 106 m3, Q = l × 105 m3/año, W = l × 106 g/año y
k = 0.25 m0.5/g0.5/año, use el método de la secante modificado para resolver para la
concentración de estado estable. Emplee un valor inicial c = 4 g/m3 y d = 0.5. Realice tres
iteraciones y determine el error relativo porcentual después de la tercera iteración.

Método de la secante por scilab.

6.26 Para el problema 6.25, la raíz puede localizarse con iteración de punto fijo como:
2
W −Qc
c= ( kV )
O bien como:

W −kV √ c
C=
Q
De las que sólo una convergerá para valores iniciales de 2 < c < 6. Seleccione la correcta
y demuestre por qué siempre lo será.

8.15 En la ingeniería ambiental (una especialidad de la ingeniería civil), la ecuación

siguiente se emplea para calcular el nivel de oxígeno c (mg/L) en un río aguas abajo de la
descarga de un drenaje:

c=10−20( e−0.2 x −e−0.75 x )

Donde x es la distancia aguas abajo en kilómetros.
a) Determine la distancia aguas abajo de la corriente, a la cual el nivel de oxígeno cae
hasta una lectura de 5 mg/L. (Sugerencia: Está dentro de 2 km de la descarga.) Encuentre
la respuesta con un error de 1%. Obsérvese que los niveles de oxígeno por debajo de 5
mg/L por lo general son dañinos para ciertas especies de pesca deportiva, como la trucha
y el salmón.
b) Calcule la distancia aguas abajo a la cual el oxígeno se encuentra al mínimo. ¿Cuál es
la concentración en dicha ubicación?
8.16 La concentración de bacterias contaminantes c en un lago disminuye de acuerdo con
la ecuación.

c=70 e−1.5 t +25 e−0.075t

Determine el tiempo que se requiere para que la concentración de bacterias se reduzca a
9 con el uso de a) el método gráfico, y b) el método de Newton-Raphson, con un valor
inicial de t = 10 y criterio de detención de 0.5%. Compruebe los resultados que obtenga.

F ( c )=70 e−1.5 t +25 e−0.075t

F ' ( c )=−105 e−1.5 t −1.875 e−0.075t

Método de newton-Raphson.

f (xi )
xi+1=xi ×
f '( xi )

70 e−1.5 (10 )+25 e−0.075 (10)

X 1+ x 1=xi=10− −9
−105 e−1.5 (10 )+.1875 e−0.075 (10)
xi=10.66
Xi+1 10.66
Xi+2 10.505
Xi+3 10.617
Xi+4 10.6404
Xi+5 10.6406
Xi+6 10.6405
Xi+7 10.6467
8.21 En ingeniería oceanográfica, la ecuación de una ola estacionaria reflejada en un
puerto está dada por l = 16, t = 12, v = 48:

2 πx 2 πtv
[ ( ) ( ) ]
h=ho sen
δ
cos
δ
+ e− x

δ =16 t=12 y v=48

Resuelva para el valor positivo más bajo de x, si h = 0.4 h0.
0.4=sen ( 216πx ) cos ( 2 π∗48
16 )
−x
+e

πx
sen
( 8)
cos(72 π )+ e −x

−0.4=0
f (x )

Método de falsa posición.

8.22 Suponga el lector que compra una pieza de equipo en $20 000 como pago inicial y
$4000 por año durante 6 años. ¿Qué tasa deinterés estaría pagando? La fórmula que
relaciona el valor presente P, los pagos anuales A, el número de años n y la tasa de
interés i, es la que sigue:

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