RAICES DE ECUACIONES NO LINEALES

1. Determine las raíces reales de f  x   0.4x2  2.2x  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 x  10 . Calcule el error estimado ε 0 y el error
verdadero ε1 para cada iteración.
2. Determine las raíces reales de f  x   2  7 x  5x2  6 x3 :
a) Gráficamente
b) Utilizando el método de bisección para localizar la raíz más pequeña. Use los valores
iniciales x1  0 y x 0  1 iterando hasta que el error estimado ε0 se encuentre debajo de
 s  10%
3. Determine las raíces reales de f  x   26  82.3x  88x2  45.4x3  9x4  0.65x5
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 x 0  1.0
c) Realice el mismo cálculo en b), pero con el método de la falsa posición y ε s  0.1%
4. Calcule las raíces reales de f  x   11  22x  17 x2  2.5x3 :
a)Gráficamente
b) Empleando el método de la falsa posición con un valor de  s correspondiente a tres cifras
significativas para determinar la raíz más pequeña.
5. Localice la primera raíz no trivial de senx  x 2 , donde x esta en radianes. Use una técnica
gráfica y bisección con un intervalo inicial de 0.5 a 1. Haga el cálculo hasta que  a sea menor
que  s  2% . Realice también una prueba de error sustituyendo la respuesta final en la
ecuación original.
6. Determine la raíz real de: lnx 2  0.7
a) Gráficamente
b) Empleando tres iteraciones en el método de bisección con los valores iniciales
x1  0.5 y xa  2
c) Usando tres iteraciones del método de la falsa posición, con los mismos valores iniciales de
b)
0.9  0.4 x
7. Determine la raíz real de f  x   :
x
a) Analíticamente b) Gráficamente.
c) Empleando tres iteraciones en el método de la falsa posición. Con valores iniciales de 1 a 3.
Calcule el error aproximado  a y el error verdadero t en cada iteración
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
9. Encuentre la raíz positiva más pequeña de la función (x esta en radianes) x 2 sen 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
 a sea menor que  s  1% . Compruebe su respuesta final sustituyéndola en la función
original.
10. Calcule la raíz real positiva de f  x   x2  8x3  36x2  462x  1010 utilizando el método de
la falsa posición. Use una grafica para escoger el valor inicial y realice el cálculo con
 s  1.0%.
11. Determine la raíz real de x3.3  79 :
a) Analíticamente
b) Con el método de falsa posición para  s  0.1% . Use como valores iniciales de 3.0 a 4.0.
gm
12. La velocidad de caída de un paracaidista está dada por: v (1  e  ct / m )
c
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% .
13. La concentración de saturación del oxigeno disuelto en agua se calcula con la ecuación (APHA,
1992)
1.575701x105 6.642308 x107 1.243800 x1010 8.621949 x1011
ln osf  139.34411    
Ta Ta2 Ta3 Ta4
Donde osf =Concentracion de saturación de oxigeno disuelto en agua a 1 atm (mg/lt) y Ta =
temperatura absoluta (K). Recuerde 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 típicas en climas templados. La ecuación sirve para determinar rangos de
concentración de oxigeno desde 14.621 mg/L hasta °0 a 6.949 mg/L hasta 35°C. Dado un valor
de concentración de oxigeno, esta fórmula y el método de bisección son útiles para resolver
la temperatura en °C.
a) Si los valores iniciales se fijan en 0 y 35°C. ¿Cuántas iteraciones del método de bisección
se requieren para determinar la temperatura con un error absoluto de 0.05°C?
b) Con base en a), desarrolle y pruebe un programa de bisección para determinar T como
función de una concentración de oxigeno dada. Pruebe el programa para osf  8,10 y 14
mg/L. Compruebe sus resultados.
14. Con el método de iteración simple de punto fijo localice la raíz de f  x   sen( x  x) . Use
un valor inicial de x o  0.5 y haga iteraciones hasta que  a  0.01% .
15. Utilice a) la iteración de punto fijo y b) el método de Newton-Raphson para determinar la raíz
de f  x   0.9x2  1.7 x  2.5 usando x0  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.
16. Determine las raíces reales de f  x   2.0  6x  4x2  0.5x3 :
a) Gráficamente b) usando el método de Newton-Raphson que cumpla con  s  0.01%
17. Emplee el método de Newton-Raphson para determinar la raíz real de
f  x   2.0  6 x  4 x2  0.5x3 , 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.
18. Determine la menor raíz real de f  x   11  22x  17 x2  25x3 : a) gráficamente y b) usando
el método de la secante para un valor de  s con tres cifras significativas.
 
19. Localice la primera raíz positiva de: f  x   senx  cos 1  x 2  1 ,donde x esta en radianes.
Use cuatro iteraciones con el método de La secante con valores iniciales de
a) xi 1  1.0 y xo  3.0, y
 b) xi 1  1.5 y xi  2.5 . Para localizar la raíz
c) Use el método grafico para verificar los resultados.
20. Calcule la raíz real de x3.3  79 , con el método de la secante modificado que cumpla con
 s  0.1% . Intente diferentes valores de  y analice los resultados.
21. Determine la mayor raíz real de f  x   x3  6x2  11x  6.1 :
a) Gráficamente b) con el método de Newton-Raphson (tres iteraciones, x i  3.5 ) c)
Utilizando el método de la secante (tres iteraciones, x i 1  2.5 y x i  3.5 d) utilizando el
método de la secante modificado (tres iteraciones x i  3.5, δ  0.02
22. Determine la menor raíz positiva de f  x   7e  x sen  x   1: a) gráficamente b) con el método
de Newton-Raphson (tres iteraciones , x i  0.3 ) c) utilizando secante (tres iteraciones,
x i 1  0.5 y x i  0.4)) d) Usando el método de la secante modificado (cinco iteraciones
x i  0.5, δ  0.03 ).
23. La función x 3  2x 2  5x  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 x0  0.2 .

24. El balance de masa para un contaminante bien mezclado en un lago se escribe así
dc
v  w  Q. c  K .v. c
dt
Dados los valores de los parámetros
v  1x106 m3 , Q  1x105 m3 / año, w  1x106 m3 / año , k  0.2m0.5 / g 0.5 / años. Use el
método de la secante modificado para determinar la concentración en estado
estacionario. Emplee como valores iniciales c  4 g / m3 y   0.5 . Realice dos iteraciones y
determine el error relativo porcentual después de la segunda iteración.
25. En el problema anterior la raíz puede localizarse con el método de iteración de punto fijo así:
W  Qc 2 W  Kv c
c( ) o así: c  solo una ecuación funciona siempre para un valor
Kv Q
inicial de c  1 Seleccione la correcta y demuestre porque siempre funciona.

26. El método de “divide y promedia”, es un antiguo método para aproximar la raíz cuadrada de
xa/x
cualquier número positivo a. se puede formular como x  . Demuestre que esta
2
fórmula es equivalente al algoritmo de Newton-Raphson para resolver x  a
27. a) Aplique el método de Newton-Raphson a la función f  x   tanh( x2  9) para evaluar su
raíz real ya conocida en x  3 . Con valor inicial use x0  3.1 y realice un mínimo de 4
iteraciones b) ¿El método muestra convergencia hacia su raíz real? Dibuja una grafica con
los resultados para cada iteración efectuada.
28. El polinomio f  x   0.0074x2  0.28x3  3.355x2  12.183x  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 x0  16.15 .
Explique sus resultados.
29. Utilice el método de la secante con la función que describe un círculo ( x  1) 2  ( y  2) 2  16
para encontrar una raíz real positiva. Tome como valor inicial xi  3 y xi 1  0.5 . Aproxime la
 solución del primer y cuarto cuadrante. Al encontrar la solución para f ( x) en el cuarto
cuadrante asegúrese de tomar el valor negativo de la raíz cuadrada. ¿Por qué la solución
diverge?
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:
Volumen de fluido que regresa a la entrada
R
volumen que sale del sistema
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 autocatalítica (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 
R(1  X Af ) R 1  R(1  X Af ) 
Donde X Af 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
X Af  0.9

Alimentación Reactor de flujo Producto

Tipo tapón

Reciclaje
Representación esquemática de un reactor de flujo tipo tapón con recirculación
31. En un proceso de ingeniería química el vapor de agua (H 2 O) se calienta a temperaturas lo
suficientemente altas para que una porción significativa del agua se disocie, o se rompa, para
formar oxigeno (O 2 ) e hidrogeno (H 2 ) :
1
H 2 O H 2  O2
2
Si se asume que esta es la única reacción que se lleva a cabo la fracción molar x de H 2 O que
x 2Pt
se disocia se representa por: K 
1 x 2  x
Donde k= la constante de equilibrio de la reacción y Pt  la presión total de la mezcla. Si
Pt  3 atm y k  0.05 , determine el valor de x que satisfaga la ecuación anterior.
32. La siguiente ecuación permite calcular la concentración de un químico en un reactor donde se
tiene una mezcla completa:
C  Cent (1  e0.04t )  Co e0.04t
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 C ent .
33. Una reacción química reversible : 2A  B C , se caracteriza por la relación de equilibrio
C
K  2c
Ca C b

Donde la nomenclatura Ci representa la concentración del componente i. Suponga que se

define una variable x que representa el número de moles de C producido. La conservación
de la masa se utiliza para reformular la relación de equilibrio como:
(Cc,0  x)
K
(Ca,0  2x) 2 (Cb,0  x)
Donde el subíndice o indica la concentración inicial de cada component.e.
Si K  0.015.Cao  42. Cbo  30 y Cco  4 , calcule x. Use bisección y falsa posición con base en
el método grafico y considere una tolerancia al error de 0.5%
34. Las siguientes reacciones químicas se llevan a cabo en un sistema cerrado
2A  B C
A D C
C Cc
En equilibrio estas pueden caracterizarse por K1  2 c , K 2 
Ca Cb Ca Cd
Donde la nomenclatura Ci representa la concentración del componente i. Si x1 y x 2 son el
número de moles de C que se producen debido a la primera y segunda reacciones,
respectivamente, emplee un método para reformular las relaciones de equilibrio en términos de
las concentraciones iniciales de los componentes (revisar el ejercicio 33). Después, use el método
de Newton Raphson para resolver el par de ecuaciones simultaneas no lineales para x1 y x 2 si
K1  4 x104 , K 2  3.7x102 , Cao  50, Cbo  20 , Cco  5 y Cdo  10 Utilice el método grafico
para proponer los valores iniciales.
RT a
35. La ecuación de estado de Redlich-Kwong está dada por: p  
v  b v(v  b) T
donde R= la constante universal de los gases =0.518KJ/Kg.°K. T=temperatura absoluta (°K).
p=presión absoluta (KPa) y v=volumen de un Kg de gas ( m3 / Kg ) .
Los parámetros a y b se calculan mediante:
R 2 Tc2.5 T
a  0.427 , b  0.0866R c
Pc Pc
Donde Pc  4600 kPa y Tc  191  K ,como ingeniero, se le pide determinar la cantidad de
combustible metano Pc=4600 kPa, Tc=191K , que se puede almacenar en un tanque de 3m3 a
una temperatura de -40°C con una presión de 65000 kPa. Emplee el método de localización de
raíces de su elección para calcular v y luego determine la masa de metano contenida en el
tanque.
36. El volumen V de un líquido contenido en un tanque horizontal cilíndrico de radio r y longitud L
está relacionado con la profundidad del líquido h por:
 rh 2
V   r 2 cos 1    (r  h) 2rh  h  L
  r  
Determine h para r =2 m , L=5 m, y V= 8m 3 . Observe que si usted utiliza un lenguaje de
programación o herramienta de software, el arco coseno se puede calcular como:
π x
cos 1 x   tan 1 ( )
2 1  x2
37. El volumen V del líquido contenido en un tanque esférico de radio r está relacionado con la
πK 2 (3r  h)
profundidad “h” del líquido por: V  . Determine h para r=1 m y V =0.5 m3 .
3
Use secante y método de falsa posición.
38. Para el tanque esférico del problema anterior es posible desarrollar las siguientes formulas
h3  (3V / ) V
para el método de punto fijo: h  , h  3 3(rh 2  )
3r 
3
Si r = 1 m y V=0.5 m determine si cualquiera de las dos alturas es estable y el rango de
valores iniciales para los que sí son estables.
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 , G 0 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.
PρDp ε 3 (1  ε)
 150  1.75
G o L (1  ε)
2 Dp G o
( )
μ
Dados los siguientes valores para los parámetros encuentre fracción vacía  del lecho
Dp G o PρDp
 1000 ;  10
μ Do2 L
H 2  0.5O2  H 2 O
41. Las siguientes tres reacciones tienen lugar en un gasificador de carbón. C  0.5O2  CO
C O2  CO2
Las expresiones de equilibrio se dan abajo .
Si P  20 atm, K1  106 , K 2  108 , K3  1014 , ntotal  3.38 mol; e1 , e2 y e3 representan el número
de moles de H 2 formadas en la reacción 1, número de moles de CO formadas en la reacción
2 y número de moles de CO2 formadas en la reacción 3, respectivamente. Usando -0.95, 1.95
y 0.01 como valores iniciales para e1 , e2 y e3 respectivamente, establezca un esquema
iteractivo y encuentre e1 , e2 y e3
(1  e1)( ntotal  0.5 e1  0.5) 0.5
K1 
(2  0.5e1  0.5 e2  e3 ) 0.5 (0.38  (e1 ) P 0.5
e2 P 0.5
K2 
(2  0.5e1  0.5e2  e3 ) 0.5 ( ntotal  0.5e1  0.5) 0.5
e3
K3 
(2  0.5e1  0.5e2  e3 ) 0.5
42. La operación de un reactor de flujo tipo tapón con densidad constante, para la producción de
una sustancia, mediante una reacción enzimática se describe con la ecuación dad abajo,
donde V es el volumen del reactor, F la velocidad de flujo de reactivo C. Cent y Csal son,
respectivamente, las concentraciones del reactivo que entran y salen del reactor, K y K max
son constantes para un reactor de 100L, con una concentración de entrada Cent  0.2 M y
una velocidad de flujo de entrada de 80 L/s, K max  102 s 1 y K  0.1 M, determine la
concentración C a la salida del reactor.
 Csal
V K 1
   dC
F Cent K max C K max
43. El desplazamiento de una estructura está definido por la siguiente ecuación para una
vibración amortiguada: y  8e kt cos wt
Donde k  0.5 y w  3 a) utilice el método grafico para realizar una primera estimación del
tiempo requerido para que el desplazamiento disminuya a 4 b) Use el método de Newton-
Raphson para calcular la raíz con ε s  0.01% c) con el método de la secante determine la raíz
al ε s  0.01%
44. La fórmula que define la fuerza por unidad de area : P/A, que causa un máximo esfuerzo σ m
en una columna que tiene una relación de esbeltez L / k es:
P σm

A 1  (ec / k 2 )sec 0.5 P / (EA)(L / k) 
 
Si E  200000MPa , ec / k  0.2 y σ m  250 MPa . Calcule 𝑃/𝐴 para L / k  100. Recuerde
2

que secx=1/cosx.
45. Un cable catenario es aquel que cuelga entre dos puntos que no están en la misma línea
vertical. Como se muestra en la figura no está sujeto a otras cargas más que a su propio peso.
Así que su peso w( N / m) actúa como una carga uniforme por unidad de longitud a lo largo
del cable. Un diagrama de cuerpo libre de una sección AB se representa en la figura inferior,
donde TA y TB son las fuerzas de tensión en sus extremos. Considerándolos equilibrios de
fuerzas horizontal y vertical, se puede obtener el siguiente modelo de ecuación diferencial:

2
d2 y w  dy 
 1  
dx 2
TA  dx 
Puede emplearse el cálculo para resolver esta ecuación para la altura “y” del cable como
función de la distancia “x”
T  w  T e x  e x
y  A cosh  x   y0  A donde cosh  x  
w  TA  w 2
Calcule TA , si w=10, y0=5, y=15, x=50

46. En la figura que se muestra se tiene una viga uniforme sujeta a una carga distribuida
creciente linealmente. La ecuación para calcular la curva elástica resultante es :
w0
y ( x 5  2 L2 x 3  L4 x )
120 EIL
Utilice el método de la bisección para determinar el punto de máxima deflexión (es decir el
valor de x donde dy / dx  0 . Después sustituya este valor en la ecuación dada para
 determinar el valor de máxima deflexión. Use los siguientes valores de los parámetros en
sus cálculos: L  450 cm , E  50000 kN / cm2 , I  30000cm 4 y w0  1.75kN / cm

47. En ingeniería ambiental, la siguiente ecuación sirve para calcular el nivel de oxigeno en un
rio aguas abajo desde una descarga de aguas residuales:
C  10  20(e0.2 x  e0.75 x )
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 esta dentro de
los 2 km de la descarga) Determine una respuesta con un 1% de error. Calcule la distancia
aguas abajo a la cual el oxigeno se encuentra al minimo. ¿Cuál es la concentración en dicha
ubicacion?
48. La concentración de bacterias contaminantes “ c ” en un lago decrece de acuerdo con la
relación c  70e 1.5t  25e 0.075t
.
Determine el tiempo requerido para que la concentración de bacterias se reduzca a 9
usando
a) el método grafico y
b) el método de Newton-Raphson considere valor inicial t=10 y tolerancia 0.5%
49. En ingeniería marítima, la ecuación de una ola estacionaria reflejada en un puerto está dada
  2x   2tw  x 
por   16 , t  12 , v  48 , h  h0  sen   cos    e  Encuentre el valor
       
positivo más bajo de x si: h  0.4h0

 2 
1  R  
50. La carga en un circuito RLC serie esta dada por q  t   q0 e Rt /(2 L ) cos  
 LC  2 L  
t
 
-6
suponga q/q0 = 0.01, t = 0.05 s, L = 5H y C = 10 F. Encuentre el valor de la Resistencia R
usando el método de Newton.
51. La carga en un circuito RLC paralelo la impedancia esta dada por
2
1 1  1 
   wC  -6
 suponga R = 225, L = 0.5H y C = 0.610 F y Z = 75 Ohms.
Z R 2  wL 
Encuentre el valor de la frecuencia w usando el método de Newton y secante.
52. El valor acumulado de una cuenta de ahorros puede calcularse con la ecuación de anualidad
vencida : A = P[(1 + i )n - 1 ] / i
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?
53. Si se compra una pieza de un equipo que cuesta $20000 al contado y en pagos de $4000 al año
durante 6 años, ¿Qué tasa de interés se está pagando? La fórmula que la relaciona el valor
presente P, los pagos anuales A el número de años “n” y la tasa de interés es:
i (1  i ) n
A P
(1  i ) n  1

54. Supongamos que la oscilación de una estructura, dotada de un sistema de amortiguación,

t
ante un movimiento oscilatorio, viene dada por la función y  10e 2
cos 2t , ¿ en que
instante t la posición de la estructura es : y(t) =4?. Utilice el punto fijo para resolver el
problema con tolerancia=0.001
log x
55. Encuentre el máximo de la función f ( x)  , x  0 . Utilice bisección con tolerancia de
x 1
0,001
1
56. La función de densidad de Cauchy tiene la forma f ( x)  , x  R , donde "  "

1   x   
2

es un parámetro
a) Genere 50 valores usando   1 y graficarlo usando Matlab.
b) Trate los datos que conseguió en el paso anterior como si fuesen simples observaciones de
la función de distribución de Cauchy con un "  " desconocido . Graficar la función de

 
n
l ()  n ln    ln 1   xi    ,  R .
2
probabilidad logarítmica asociada:
i 1
c) Calcular el máximo valor de la función anterior usa el método de falsa posición para hallar un
estimador de "  " con tolerancia de 0.0001 (para hacer esto se requiere de la derivada :
n   xi
l '()   ,  R

i 1 1   x   2
i 
y graficarla).

57. La función de distribución de Weibull es de la forma: F ( x)   
1  exp    x  , si x  0
,
 0, en otro caso
donde  ,  son parámetros positivos.
a) Genere 50 números aleatorios de la función de distribución de Weibull con   1.8 ,   1

b) Adicione 3 números más a los ya calculados anteriormente. Trate a los 53 datos como si
fuesen observaciones de la función de Weibull con un  desconocido , pero mantenga   1 .
Plotee la función de probabilidad y la de probabilidad logarítmica (dependiente de  ).
Considere las funciones de probabilidad:
1
 n   n 
L()   n   xk  .exp    xk  
 k 1   k 1 
 n n
y la función de densidad de probabilidad logarítmica l ()  n .ln   (  1)  ln xk   xk 
k 1 k 1
c) Use el método de Newton para hallar el máximo estimador de probabilidad  . Para hacer
esto, necesitamos resolver la ecuación l '()  0 , para los puntos estacionarios. Considere que
la primera y segunda derivada de l () son:
n n n
l ' ( )  .   ln xk   xk  ln xk
 k 1 k 1
n n
.   xk   ln xk 
2
l '' ( )  
 2
k 1
(Use el criterio de la segunda derivada para probar que es un máximo).

58.Muchos campos de la ingeniería requieren estimaciones exactas de la población, por ejemplo,

los ingenieros en transporte pueden considerar necesario determinar por separado la
tendencia del crecimiento demográfico de una ciudad y de los suburbios adyacentes. Si la
población del area urbana disminuye con el tiempo de acuerdo con: P  t   Pmax e
 kt
 P.min

Ps.max
Mientras que la población suburbana crece de acuerdo con: Ps  t  
1   Ps.max / P0  1 e  K s t

Donde P.max , k , Ps.max , P0 y K s son parámetros obtenidos en forma empírica. Determine el

tiempo y los valores correspondientes de P (t ) , Ps (t ) cuando la población en la ciudad
sea 20% mayor que la suburbana. Los valores de los parámetros son Pμmax  75000 ,
K μ  0.045 / año , Pμmin  100000 personas , Ps.max  300000 personas, P0  10000
personas, K s  0.08 / año K s = 0.075/𝑎ñ𝑜. Para obtener las soluciones, use a) un
método grafico, b) el método de la falsa posición y c) el método de la secante modificada.

59. Una corriente alterna en un circuito eléctrico se describe mediante I  9et sen(2 t ), donde “t”
esta en segundos. Determine todos los valores de t tales que I  3.5 . Use punto fijo y falsa
posición.

60. La resistividad  de un silicon revestido depende de la carga “q”en un electron, la densidad del
electron “n”y la movilidad del electrón ”  ”. La densidad del electrón está dada en
términos de la densidad de revestimiento “N” y la densidad portadora intrínseca “ ni ”.La
movilidad del electrón está definida por la temperatura ”T”, la temperatura de referencia
“ T0 ” y la movilidad de referencia “  0 ”. Las ecuaciones necesarias para calcular la resistividad
1 1 T 2.42
son:   , donde n  ( N  N 2  4ni 2 ) y   0 ( )
qn 2 T0

Determine N, dados T0  300 K , T  1000K , 0  1350cm2 Vs 1 , q  1.7 x1019 C 

, y ni  6.21x109 cm3 una resistividad deseada   6.5 x106 Vs cm / C Use los métodos a) de
bisección y b) de la secante
61.Una carga total Q se distribuye en forma uniforme alrededor de un conductor con forma de anillo
circular con radio “a”. Una carga q se localiza a una distancia x del centro del anillo (ver figura).
La fuerza ejercida sobre la carga por el anillo está dada por:

1 qQx
F
4e0 ( x  a 2 )3/2
2

Donde e0  8.85 x1012 C 2 / ( Nm2 ). Encuentre la distancia “x” donde la fuerza es de 1 N si “q”
y ”Q” son 2 x105 C para un anillo de 0.9m

62. En la figura de la parte inferior se muestra un circuito con un resistor, un inductor y un capacitor
en paralelo. Las reglas de Kirchhoff sirven para expresar la impedancia del sistema como:

1 1 1 2
 2
 ( wC  )
Z R wL

Donde Z  impedancia  Ω  y w=la frecuencia angular. Encuentre  para que la impedancia

resultante sea de 100Ω, usando los métodos de la bisección y de la falsa posición con valores
iniciales de 1 y 1000 para los siguientes parámetros: R  225 Ω , C  0.6 x106 F y L  0.5H

63. Para el flujo de fluidos en una tubería, la fricción se describe mediante un número adimensional:
el factor de fricción de Fanning f, el cual depende de varios parámetros relacionados con el fluido
y el tamaño de la tubería, que se pueden representar por otra cantidad adimensional: el número
de Reynolds Re. Una fórmula que predice f dado Re es la ecuación de von Karman:

1
f

 4log10 Re f  0.4 
Los valores típicos del número de Reynolds para un flujo turbulento van de 10000 a 500000, y
los del factor de fricción de Fanning van de 0.001 a 0.01. Calcule f si el usuario introduce valores
entre 2500 y 1000000. (use bisección y falsa posición)

64. Los sistemas mecánicos reales llegan a involucrar la deflexión de resortes no lineales. En la figura
P8.34, una masa m se suelta desde una distancia h sobre un resorte no lineal. La fuerza de
resistencia F del resorte está dada por: F  (k1d  k2 d 3/2 )
 2k2 d 5/2 1
Con la conservación de la energía se demuestra que: 0  k1 d 2  mgd  mgh
5 2

Encuentre 𝑑 dados los siguientes valores de los parámetros: .k1  40000 g / s 2 ,

k2  40 g / ( s 2 m0.5 ) , m  95 g , g  9.8m / s 2 y h  0.43m .

65.Los ingenieros mecánicos, así como los demás ingenieros, usan la termodinámica en su trabajo.
El siguiente polinomio se puede usar para relacionar el calor específico a presión cero
del aire seco, c p kJ / (kg K ), con la temperatura.

c p  0.99403  1.671x104 T  9.7215 x108 T 2  9.5838 x1011 T 3  1.9520 x1014 T 4

Determine la temperatura que corresponda a un calor especifico de 1.2kJ / (kg K ) .

66.Los ingenieros en aeronáutica suelen calcular las trayectorias de proyectiles como los cohetes.
Un problema relacionado con dicho tema es la descripción de la trayectoria de una pelota. La
trayectoria de una pelota lanzada por un jugador que se encuentra en el jardín derecho está
definida por las coordenadas (x,y) como se presenta en la fig P8.36. La trayectoria se puede
modelar como:

g
y   tan0  x  x 2  1.8
2v02 cos 2 0

Calcule el angulo inicial 0 , si la velocidad inicial es 20 m/s y la distancia “x” al cátcher es de 40 m.

Observe que la pelota sale de la mano del lanzador con una elevación y0=1.8m y el cátcher la recibe
a 1m. Exprese el resultado final en grados (considere g=9.8m/s2)

67. La velocidad hacia arriba de un cohete se calcula con la siguiente fórmula:

 m0
v  ln  gt
m0 qt

Donde v  velocidad hacia arriba,   la velocidad con que el combustible sale del cohete, m0 
masa inicial del cohete en el tiempo, t  0 , q  el consumo de combustible y la aceleración hacia
abajo debido a la gravedad (considere la gravedad constante = 9.8m/s 2. Si   2200 m / s ,
m0  160000 kg y q  2680 kg / s , calcule el tiempo en que v  1000 m / s (Sugerencia: “t”está
entre 10 y 50 segundos). Determine el resultado dentro del 1% del valor verdadero.
Compruebe su respuesta.

68. En un tanque mezclador entran dos líquidos a diferente temperatura y salen a la misma
temperatura. La capacidad calorífica del fluido A está dada por:

c p  3.381  1.804 x102 T  4.300 x106 T 2

Y la capacidad calorífica del fluido B está dada por: c p  8.592  1.290 x101 T  4.078 x105 T 2

Donde las unidades de 𝑐𝑝 son cal / mol K , y las de T esta en unidades de K Recuerde que:
T2
ΔH   c p dT
T1

El fluido A entra al mezclador a 400°C y el B entra al mezclador a 600 °C. La cantidad de A que
entra al mezclador es el doble de la cantidad de B. ¿A qué temperatura salen los dos fluidos del
mezclador?

69. Un compresor está operando con una razón de compresión Rc , de 3.0 (la presión del gas que
sale es tres veces mayor que la presión que el gas que entra). Los requerimientos de energía del
compresor Hp, se determina mediante la ecuación dada abajo. Suponiendo que los
requerimientos de energía del compresor son exactamente igual a zRT1 / MW , encuentre
la eficiencia poli trópica n del compresor. El parámetro z es la compresibilidad del gas bajo las
condiciones de operación del compresor, R es la constante de los gases, T1 es la temperatura
del gas a la entrada del compresor y MW es el peso molecular del gas.

zRT1 n  n 
n 1
Hp   Rc  1
MW n  1  


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


(x  4)  (y 4)  5
2 2
a)  use Newton

 x 2  y 2  16


 y   x  x  0.75
2
b)  use x=y=1.2 y analize los resultados.(Use Newton)

 y  5 xy  x 2
  y  x 2 1
c)  use secante
 y  2cos x

71. Con ayuda de funciones de singularidad es posible calcular el esfuerzo cortante a lo largo de
una viga , expresándola en la siguiente forma:

1 1 0
V (x)  20[ x  0  x  5 ]  15 x  8  57 . Por definición , la función de singularidad se
expresa del modo siguiente:

n (x  a) n , x  a
xa  , halle los puntos en los que el esfuerzo cortante es cero.
 0 , x  a

72.La forma general para un campo tensorial de tres dimensiones es la siguiente:

  xx  xy  xz 
 
  xy  yy  yz 
  xz  yz  zz 


En la que los términos en la diagonal principal representan esfuerzos a la tensión o a la

compresión y los términos fuera de la diagonal representan los esfuerzos cortantes. Un campo
tensorial (en MPa) esta dado por la siguiente matriz:

10 14 25
14 7 15 
 
 25 15 16 

Para resolver cuales son los esfuerzos principales, es necesario construir la matriz siguiente (en
MPa):

10   14 25 
 14 7 15 

 25 15 16   

1 , 2 , 3 se obtienen con la ecuación : 3  I 2  II   III  0 , donde:

I   xx   yy   zz
II   xx  yy   xx  xx   yy  zz   2 xy   2 xz   2 yz
III   xx  yy  zz   xx 2 yz   yy  2 xz   zz  2 xy  2 xy  xz  yz

I,II y III se conocen como las invariantes de esfuerzos. Estime 1 , 2 , 3 , usando bisección y
secante.
  Subir el archivo digital al aula virtual (resuelto en Matlab ) el domingo 05 de mayo e
impreso en hora de clase el lunes 06 de mayo

N° Ejercicios asignados
APELLIDOS Y NOMBRES
1 AGUERO RODRIGUEZ, MARIA FERNANDA 1 29 60

2 BERNABE RUIZ, DANIEL 3 31 61

3 CHAPI SUYO, ISAIAS ALDO 41 33 63

4 DURAN RIVERA, MARY CECILIA 11 51 65

5 ESPINOZA REYNAFARJE, PATRICK LUIS 4 34 62

6 FIGUEROA VALERIANO, DAN JAROD 5 36 66

7 GABRIEL RODRIGUEZ, JAKELINE STHEFANY 8 32 64

8 HUAMAN UMERES, CRISTIAN MARCELO 7 37 68

9 LOBATON TARAZONA, GRECIA ISABEL 9 35 70

10 QUISPE AVENDAÑO, KELVIN MISAEL 10 38 67

11 RAMIREZ ALBERTO, ANGIE LIZBETH 12 55 69

12 RAMIREZ CABELLO, MARIA GRAZIA 18 2 1

13 RAMOS MORALES, JOICE YORELI 6 42 2

14 SALVADOR TORIBIO, PAUL KENNEDY 14 40 3

15 SALVATIERRA CAPCHA, KATERIN 13 57 4

16 SILVERA GARCIA, KATTYA ROSSCELYN 11 44 5

17 URIARTE BARRAZA, KATHERIN LOURDES 46 15 6

18 VERDE ROBLES, ANYELA BRILLYTH 17 60 7

19 MORAN PISCO JULIO CESAR 16 39 8

20 ALIAGA ZEGARRA, OLIMBER OCTAVIO 19 49 9

21 BALTAZAR MARIANO, ALEXIS DIOWIL 22 50 10

22 HUAMAN CARDENAS, THEO ALESSANDRO 20 52 11

23 LIJARZA RAFAILE, YEISON YENCO 23 48 12

24 LLANA ANAYA, YEIDI MAYLIN 21 51 13

25 MACHCO HUAMAN, VIDAL FERNANDO 56 53 14

26 MARAVI PEREZ, ANTHONY MANUEL 25 45 16

27 PALMA MAGARIÑO, CESAR RENATO 27 47 18

28 PORTILLA CARRASCO, ANATOLLY EULER 28 54 15

29 SERAFIN FAUSTINO, ESTER 26 57 19

30 VARGAS FLORES, FRANCIS STARLYN 41 56 20

31 VELA JARAMILLO, SUMMER NAIN 30 59 17