Universidad Aut

onoma de Nuevo Le
on

UANL-FCT

Facultad de Ciencias de la Tierra

M
etodos Num
ericos
Evidencia 01

Soluci
on de Ecuaciones no Lineales
January 27, 2016

Instrucciones: En cada uno de los problemas encuentre por los metodos de bisecci
on, regla falsa, sectante,
punto fijo y Newton-Raphson y alg
un metodo adicional
no visto en clase, las races, extremos relativos y puntos de inflexi
on en el intervalo dado, haga gr
aficos en R
para cada una de las funciones se
nalando donde se encuentran estos puntos, las gr
aficas de su derivada y su segunda derivada, grafique sobre la funci
on dada los puntos
m
aximos o mnimos y estime con cinco cifras decimales
el valor de la variable independiente donde esta tiene un
m
aximo o mnimo y los puntos de inflexi
on, que valor
tiene este. Escriba un enunciado sobre la conclucion del
problema. En caso de que la funci
on no tenga extremos
relativos ni puntos de inflexi
on en el intervalo indicado,
da argumentos s
olidos de por que no los tiene. Haga
adem
as una descripci
on de las funciones.
1. f (x) = 1+2x3x2 ex +2x3 (sin x)ex/5 , x [4, 20]
2. f (x) = sin x cos(x2 + 1) 1,
3. f (x) = x2 ln x 9x + 18,

x [, ]
x (0, 8]

14. f (x) = 2000x3 e5x cos 2x + 1,

15. f (x) = x 2 cos x,

x [2, 2]

16. f (x) = 5xex + cos 5x,

x [0.2, 2]

17. f (x) = x2 cos 10x + 1,

18. f (x) = 23000x3 e0.01x + sin 5000x,
[0.05, 0.05]
19. f (x) = 5x cos x + 2,

21. f (x) = 8x + cos(x + ) x2 ,

22. f (x) = ln x

x2
5
+ ,
3
3

23. f (x) = x2 ex + ln x 3,

x
x [, ]

x [0, 3]
x (0, 4]
x (0, 40]

24. f (x) = sin(0.1x) + e2x 0.6,

x [1, 100]
x [9, 11]

x [, ]

25. f (x) = cos 3x + 5e0.01x 3,

5. f (x) = 2x tan x 10,

x [, ]

26. f (x) = x2 cos x 6x ln x 25,

x (0, 8]

x [1, 2]

20. f (x) = e2x sin(0.1x) e10x cos(0.1x) + 0.3, x

[0, 1/2]

4. f (x) = x3 2x2 sin x,

6. f (x) = x ln x 10,

x [0.15, 2.2]

27. f (x) = (x 2)2 ln x,

x [30, 40]
x (0, 4]

7. f (x) = 2x0.2 ex tan x,

x [0, /2)

28. f (x) = cos x 3x,

8. f (x) = 2x0.6 (cos x) log10 x,

x [35, 55]

29. f (x) = cos 5x + 5 sin 15x + e0.05x 3, x [17, 18]

9. f (x) = x3 ex + 4x2 10,

x [2, 3]

10. f (x) = x3 cos x 4x2 1,

x [20, 40]

30. f (x) = log10 (50x 40) x1.2 + 20 cos x + 12, x

(4/5, 10]

11. f (x) = 2x 20 sin x,

12. f (x) = 1000x2 (tan x2 )e5x 10,
3

13. f (x) = x 2x cos 2x 12,

x [2, 2]
x [0, /2]
x [0, 4]

x [2, 2]

31. p(x) = x7 10x6 + 14x5 + 140x4 511x3 + 350x2 +

496x 480
32. p(x) = x6 10x5 60x4 + 646x3 + 491x2 7116x +
6048

Compiladores: R. Soto-Villalobos, J.A. Briones-Carrillo, B.E. De Leon-Jasso

Universidad Aut
onoma de Nuevo Le
on

UANL-FCT

Facultad de Ciencias de la Tierra

33. p(x) = x7 12x6 + 42x5 + 12x4 375x3 + 840x2

748x + 240

59. p(x) = x6 +11.25x5 +47.875x4 +97.5x3 +99.625x2 +

48.75x + 9

34. p(x) = x7 4x6 62x5 +200x4 +1009x3 2356x2

3828x + 5040

60. p(x) = x7 + 16.3x6 + 91.31x5 + 233.329x4 +

294.262x3 + 184.7227x2 + 53.15896x + 5.25504

35. p(x) = x5 x4 43x3 + 61x2 + 342x 360

61. p(x) = x5 + 2x4 + 22x3 + 160x2 + 69x + 132

36. p(x) = x5 2x4 23x3 112x2 44x + 1680

62. p(x) = x7 + 28x6 + 322x5 + 1960x4 + 6769x3 +

13132x2 + 13068 + 504

37. p(x) = x6 17x5 +80x4 90x3 +249x2 853x+630

38. p(x) = x6 +19x5 +154x4 +664x3 +1543x2 +1765x+
750
39. p(x) = x5 + 43x4 + 666x3 + 4502x2 + 12293x + 8415

63. p(x) = x7 + 9.21x6 + 32.892x5 + 57.528x4 +

50.772x3 + 20.502x2 + 2.6x + 0.024
64. p(x) = x4 + 29x3 + 288x2 + 1116x + 1296

40. p(x) = x8 + 52x7 + 1079x6 + 11638x5 + 71159x4 +

252088x3 + 503401x2 + 513822x + 201960

65. p(x) = x4 + 17x3 + 101x2 + 247x + 210

41. p(x) = x5 + 26x4 + 254x3 + 1156x2 + 2433x + 1890

67. p(x) = x6 + 59x5 + 1330x4 + 14290x3 + 74129x2 +

164851x + 103740

42. p(x) = x5 + 4x4 375x3 + 2510x2 856x 14784

68. p(x) = x8 + 43x7 + 784x6 + 7882x5 + 47509x4 +

174307x3 + 375066x2 + 422568x + 181440

43. p(x) = x5 + 3x2 63x3 + 13x2 + 726x 1080

44. p(x) = x5 + 27x2 + 269x3 + 1197x2 + 2250x + 1296
6

66. p(x) = x5 17x4 + 107x3 307x2 + 396x 180

69. p(x) = x4 + 34x3 + 421x2 + 2244x + 4320

45. p(x) = x + 34x + 440x + 2690x + 7755x +

9036x + 3564

70. p(x) = x5 33x4 + 406x3 2262x2 + 5353x 3465

46. p(x) = x4 + 23x3 + 182x2 + 568x + 756

71. p(x) = x5 + 37x4 + 511x3 + 3242x2 + 9280x + 9100

47. p(x) = x7 + 4x6 + 6x5 + 8x4 + 145x3 + 6768x2 +

2425x + 6454

72. p(x) = x6 + 53x5 + 1082x4 + 10670x3 + 51809x2 +

110237x + 68068

48. p(x) = x10 + 9x9 + 54x8 + 76x7 + 34x6 + 23x5 +

567x4 + 89x3 + 23x2 + 34x + 23

73. p(x) = x5 12x4 + 35x3 + 80x2 + 3x + 77

49. p(x) = x 6x 12x + 90x + 39x 324x

28x + 240
50. p(x) = x5 x4 27x3 + 13x2 + 134x 120
4

51. p(x) = x + 184x + 11386x + 316704x + 2995200

52. p(x) = x6 + 49x5 + 945x4 + 9035x3 + 43874x2 +
96216x + 60480
10

74. p(x) = 4x5 + 56x4 12x3 + 34x2 + 9x 14

75. p(x) = x5 13.999x4 +74.99x3 159.959x2 +1.95x+

150
76. p(x) = 5x4 + 55.998x3 260.982x2 + 463.934x
199.95
77. p(x) = 23x6 4x5 + 23x4 + +21x3 + 9x2 + 5x 8
78. p(x) = 4x6 + 7x5 + x4 4x3 + 8x2 + 12x 5

53. p(x) = x + 103x + 4722x + 126798x +

2206785x6 + 25985127x5 + 209407148x4 +
1138752932x3 + 3991620384x2 + 8123774400x +
7264857600
54. p(x) = x4 + 30x3 + 305x2 + 1200x + 1404
55. p(x) = x4 + 27x3 + 226x2 + 648x + 448
56. p(x) = x6 + 45x5 + 802x4 + 7236x3 + 34792x2 +
84384x + 80640

79. p(x) = 2.453x6 7.434x5 1.23x4 + 21.98x3 +

2.679x2 + 0.53x 2.456
80. p(x) = x4 2.0379x3 15.4245x2 + 15.6696x +
35.4936
81. p(x) = 1.234x6 + 0.277x5 5.235x4 4.21x3 +
5.348x2 5.12x 3.4565

57. p(x) = x8 + 64x7 + 1708x6 + 24640x5 + 208054x4 + Miscelanea de Problemas

0.5cm1038016x3 +2924172x2 +4098240x+2027025
58. p(x) = x6 + 42x5 + 700x4 + 5880x3 + 25948x2 +
56448x + 46080

82. Apply Newtons method to find a root near x = 2

of x4 5x3 12x2 + 76x 79 = 0

Compiladores: R. Soto-Villalobos, J.A. Briones-Carrillo, B.E. De Leon-Jasso

Universidad Aut
onoma de Nuevo Le
on

UANL-FCT

83. Apply Newtons method to find two roots near

x = 1.5 of x6 = x4 + x3 + 1
84. Apply Newtons method to find all the roots of
x3 13x2 22x + 3 = 0
85. Apply Newtons method to find a root near 1.5 of
x = 12 sin x
86. Apply Newtons method to split the fourth-degree
polynomial x4 + 5x3 + 3x2 5x 9 = 0 in two
quadratic factors.
87. Apply Newtons method to find the smallest real
root of ex = sin x

Facultad de Ciencias de la Tierra

99. Lee y Duffy (1976) relacionan el coeficiente

de friccion para el flujo de una suspensi
on de
partculas fibrosas con el n
umero de Reynolds mediante la siguinente ecuacion emprica:


p
1
5.6
1
= ln(RE f ) + 14
k
k
f
Em esta relacion f es el cooeficiente de fricci
on, RE
es el n`
umero de Reynolds y k es una constante determinada por la concentraci`on de la suspensi
on.
Para una suspension de 0.08% de concentraci
on,
k = 0.028. Cual es el valor de f si RE = 3750?
100. Se tiene
x00 + x + 2y 0 + y = f (t)

88. Apply Newtons method to find to twelve digits a

root near 2 to x4 x = 10
89. Find de Maximum value of y(x) near x = 1, given
that sin(xy) = y x
90. Find the smallest positive root of
1x+

x2
x3
x4

+
= 0
2
2
(2!)
(3!)
(4!)2

91. Where does the curve y = e

tion point?

ln x have an inflec-

92. Find the minimum for positive x of y = (tan x)/x2 .

93. Find a root near x = 2 of x = sin 2x
94. Find a small positive root of x1.8632 = 5.2171x
2.1167
95. Find any real roots of x3 = x + 4

x00 x + y = g(t)
x(0) = x0 (0) = y(0) = 0
Para resolver este par de ecuaciones diferenciales
de segundo orden simultaneas aplicando la transformada de Laplace se hace necesario factorizar la
experesion:
s3 s2 + 3s + 1
de modo que para obtener la transformada inversa
pueden usarse fracciones parciales. Cu`ales son estos factores?
101. Con base en el trabajo de Frank-Kamenetski
(1955), las temperaturas en el interior de un material con fuentes de calor incrustadas puede determinarse si se resuelve la siguiente ecuacion:
r
1
1 0.5t
0.5t
Lct
e
cosh (e ) =
2
Dado Lct = 0.008, encontrar t.

96. Find two roots near x = 1 of 2x4 + 16x3 + x2

74x + 56 = 0
97. DeSantis(1976) dedujo una relaci
on para el factor
de compresibilidad de gases reales de la forma
z=

1 + y + y2 y3
(1 y)3

con y = b/4v, donde b es la correcci

on de Van Der
Walls y v es el volumen molar. Si z = 0.892 cu`al
es el valor de y?
98. En estudios sobre recolecci
on de energa solar al
enfocar un campo de espejos planos en un colector central, un imvestigador obtuvo la siguiente
ecuaci
on para el factor de concentraci
on geometrica
C:
(h/ cos A)2 F
C=
2
0.5D (1 + sin A 0.5 cos A)

102. La solucion de problemas con valor en la frontera

por medio de series de Fourier a menudo comprende la detarminacion de races de ecuaci
ones
trasendentes para evaluar lso coeficientes. Por
ejemplo,
y 00 + y = 0,

y(0) = 0,

y(1) = y 0 (1)

comprende la solucion tan z = z. Encuantra tres

valores distintos de z = 0.
103. La razon de flujo de agua en una corriente a
menudo se meide con la instalacion de un arepresa.
Esto supone la construccion de un dique a travez
de la corriente con una muesca en forma de V
cerca del centro a punta de la muesca hacia
abajo). Si se despreia la velocidad contracorriente,
el flujo Q(pies3 /s) esta relacionado con la distancia
h(pies) de la superficie del agua a contra corriente

Compiladores: R. Soto-Villalobos, J.A. Briones-Carrillo, B.E. De Leon-Jasso

Universidad Aut
onoma de Nuevo Le
on

UANL-FCT

Facultad de Ciencias de la Tierra

al punto de la V y con el
angulo (grados) entre
los lados de la muesca mediante la ecuacion

 
p
8
tan
2gh2.5
Q = 0.59
15
2

de la seccion transversal se relacionan con la profundidad y por medio de:

donde g es la constante gravitacional. si Q = 200,

elabora uan tabla que muestre como esta relacionada h con par avalores de entre 20 y 130
grados.

Resuelva para la profundidad crtica.

104. In the year 1225 Leonardo of Pisa estudied the ecuation

p(x) = x3 + 2x2 + 10x 20
and produced x = 1.386808107. No body knows
by what method Leonardo found this value but it
is remarkable result for his time. Apply Newtons
method to Leonardos equation.
105. Show that the second-drgree term is removed from
the general cubic equation x3 + ax2 + bx + c = 0
by traslation x = y a/3. In 1545 Cardano published this formula for solving the cubic equation
x3 + bx + c = 0:

"

1
x= c+
2

#1/3
1 2 1 2
c + b

4
9
"
#1/3
r
1
1 2 1 2
c+
c + b
2
4
9

Apply it to find at least the real root x = 1 of

x3 15x 4 = 0. Apply Newtons method to Cardanos equation.
106. Encuentre la
raz positiva m
as peque
na de la
funci
on x2 | cos x| = 5
107. La velocidad v de un paracaidista que cae esta dada
por
gm
v=
(1 e(c/m)t )
c
donde g = 9.8 m/s2 . Para un paracaidista con coeficiente de arrastre de c = 15 kg/s, calcule la masa
m de modo que la velocidad sea v = 35 m/s en
t = 9s.
108. Por un canal trapezoidal fluye agua a una tasa de
Q = 20 m3 /s. La profundidad crtica y para dicho
canal satisface la ecuaci
on
y =1

Q2
B
gA3c

donde g = 9.81 m/s2 , Ac ,

area de la seccion
transversal (m2 ), y B es el ancho del canal en la
superficie (m). Para este caso, el ancho y el area

b=3+y

Ac = 3y +

y2
2

109. Suponga que esta dise

nando un tanque esferico
para almacenar agua para un poblado peque
no en
un pas en desarrollo. El volumen de lquido que
puede contener se calcula con
V = h2

3R h
3

donde V es el volumen en m3 , h es la profundidad

del agua en el tanque en m, y R es el radio del
tanque en m. Si R = 3m, a que profundidad debe
llenarse el tanque de modo que contenga 30m3 ?
110. La concentracion de saturacion de oxgeno disuelto
en agua dulce se calcula con la ecuacion (APHA,
1992)
1.575702 105
Ta
7
1.243800 1010
6.642308 10
+

2
Ta
Ta3
8.621949 1011

Ta4

ln Osf = 139.34411 +

donde Osf es la concentracion de saturaci

on de
oxgeno disuelto en agua dulce a 1 atm (mg/L) y Ta
es la temperatura absoluta (K). Recuerde el lector
que Ta = T + 273.15, donde T es la temperatura en
grados cent`grados. De acuerdo con esta ecuaci
on,
la saturacion disminuye con el incremento de la
temperatura. Para aguas naturales comunes en
climas templados, la ecuacion se usa para determinar que la concentracion de oxgeno vara de 14.621
mg/L a 0 grado centgrados a 6.413 mg/L a 40 grados centgrados. Dado un valor de concentraci
on de
oxgeno, puede emplearse esta formula y el metodo
de Newton para resolver para la termperatura en
grados centgrados.
111. El balance de masa de un contaminante en un lago
bien mezclado se expresa as:
V

dc
= W Qc kV c
dt

Dados los valores de parametros V = 1 106 m3 ,

Q = 1 105 m3 /anno y W = 1 106 g/anno, y
k = 0.25m0.5 /anno, use el metodo de Newton para
resolver para la concentracion de estado estable.
Emplee un valor inicial c = 4g/m3 y = 0.5.

Compiladores: R. Soto-Villalobos, J.A. Briones-Carrillo, B.E. De Leon-Jasso

Universidad Aut
onoma de Nuevo Le
on

UANL-FCT

112. The Van Der Waals equation of state for a vapor is


a
P 2 (v b) = RT
v
2

where P is the pressure (N/m ), v is the specific

volume (m3 /kg), T is the temperature (K), R is
the gas constant (J/kg K), and a and b are empirical constants. Consider water vapor, for which
R = 461.495J/kg K, a = 1703.28P a (m3 /kg)3 ,
and b = 0.00169099(m3/kg). The equation can be
rearranged into the form (show it!!)
P v 3 (P b + RT )v 2 + av ab = 0
Calculate the specific volume v for P = 10000kP a
and T = 800K. Use the ideal gas law, P v = RT ,
to obtain the initial guess (or guesses).
113. When an incompressible fluid flows steadily
through a round pipe, the pressure drop due to
the effects of wall friction is given by the empirical
formula:
 
L
2
P = 0.5f V
D
where P is the pressure drop, is the density,
V is the velocity, L is the pipe length, D is the
pipe diameter, and f is the DArcy friction coefficient. Several empirical formulas exist for the friction coefficient f as a function of the dimensionless
Reynolds number, Re = DV /, where is the
viscosity. For flow in the turbulent regime between
completely smooth pipe surfaces and wholly rough
pipe surfaces, Colebrook (1939) developed the following empirical equation for the friction coefficient
f:


1
/D
2.51
= 2 log10
+
3.7
f 1/2
Ref 1/2
where  is the pipe surface roughness. Develop a
procedure to determine f for specified values of
/D and Re. Use the approximation proposed by
Genereaux (1939) to determine the initial approximation(s):
f = 0.16Re0.16

Facultad de Ciencias de la Tierra

Solve for f for a pipe having /D = 0.001 for

Re = 10n , for n = 4, 5, 6, 7.
114. Consider quasi-one-dimensional isentropic flow of a
perfect gas through a variable-area channel. The
relationship between the Mach number M and
the flow area A, derived by Zucrow and Hoffman
(1976), is given by
1
A
= =
A
M

2
+1

1 2
1+
M
2

(+1)/2(1)

where A is the choking area (i.e., the area where

M = 1) and is the specific heat ratio of the flowing gas. For each value of , two values of M exist,
one less than unity (i.e., subsonic flow) and one
greater than unity (i.e., supersonic flow). a) Calculate both values of M for  = 10.0 and = 1.4
by Newtons method. For the subsonic root, let
M1 = 0.2. For the supersonic root, let M1 = 5.0.
b) For the subsonic root, let M1 = 0.4. For the
supersonic root, M1 = 5.0.
115. La ley de los gases ideales esta dada por
pV = nRT
donde p es la presion absoluta, V es el volumen,
n es el n
umero de moles, R es la constante universal de los gases y T es la temperatura absoluta.
Aunque esta ecuacion se utiliza ampliamente por
los ingenieros y cientficos, solo es exacta en un
rango limitado de presion y temperatura. Adem
as,
la ecuacion es apropiada solamente para algunos
gases. Una ecuacion de estado alternativa para los
gases esta dada por:

a
p + 2 (v b) = RT
v
conocida como la ecuacion de Van Der Waals,
donde v = V /n es el volumen molar, a y b son
constantes empricas que dependen del gas que se
analiza. del chapra pagina 200

Compiladores: R. Soto-Villalobos, J.A. Briones-Carrillo, B.E. De Leon-Jasso