UNIVERSIDAD NACIONAL HERMILIO VALDIZAN - HUANUCO

FACULTAD DE INGENIERIA INDUSTRIAL Y DE SISTEMAS

ESCUELA PROFESIONAL DE INGENIERIA INDUSTRIAL

CUADERNO VIRTUAL

ALUMNO:
❖ LEYVA HUAMAN ANTHONY CESAR

CURSO:
❖ METODOS NUMERICOS

HUANUCO – PERÚ
2020
 CAPÍTULO 3
3.1 Convierta los números que se presentan base 2 a base 10:
a)101101;

x (base 10) =1*2^0 + 0*2^1 + 1*2^2 + 1*2^3 + 0*2^4 + 1*2^5

x (base 10) =45

b) 110.011

x (base 10) =1*2^2 + 1*2^1 + 0*2^0 + 0*2^0 + 1*2^-1 + 1*2^-2

x (base 10) =6.375

c) 0.01101.

x (base 10) =0*2^0 + 0*2^0 + 1*2^-1 + 1*2^-2 + 0*2^-3 + 1*2^-4

x (base 10) =0.40625

3.2 Convierta los siguientes números base 8 a base 10:

a) 71 263

x (base 10) =3*8^0 + 6*8^1 + 2*8^2 + 1*8^3 + 7*8^4

x (base 10) =29363

b) 3.147

x (base 10) =3*8^0 + 1*8^-1 + 4*8^-2 + 7*8^-3

x (base 10) =3.201171875

3.3 La derivada de f(x) = 1/(1 – 3x2) está dada por (6*x)/((1-3x^2)^2), ¿Esperaría
dificultades para evaluar esta función para x = 0.577? Inténtelo con aritmética de
3 y 4 dígitos con corte.

X X+0.1111 X+0.1111+0.1111 …
0.577 0.6881 0.7992 …
FX = 6x/(1-3x^2) ^2 2352910.793 23.3552637 5.712974066 …

Respuesta: se observa que al evaluar para 0.577 se obtiene 2352910.8, pero al manejar 4
decimales esta cambia repentinamente dando como resultado un error en escala diminutiva.

3.4 Evalúe el polinomio y = x3 – 5x2 + 6x + 0.55

a) En x = 1.37. Utilice aritmética de 3 dígitos con corte. Evalúe el error relativo
porcentual.

x x+0.0001 …
 1.37 1.3701 1.3711 1.3721 …
y = x3 – 5x2 + 6x + 0.55 1.2624 1.26190396 1.25693916 1.25196636 …
4% 39% 40% …
Error% Error% Error% Error% …

b) Repita el inciso a) pero exprese a y como y = ((x – 5) *x + 6) *x + 0.55 y evalúe el

error y compárelo con el inciso a).

y = ((x – 5)x + x x+0.0001

6)x
1.37 1.3701 1.3702 1.3703 1.3704 1.3705 1.3706
0.7124 0.71190396 0.71140784 0.71091164 0.71041536 0.709919 0.70942256
7% 7% 7% 7% 7% 7%
Error%
dato
constante =
0.55

De acuerdo con lo que podemos comparar para datos elevados en el inciso b) que sería en
avance iteracional, la variación en el error es constante, mientras en la b) aumenta
agigantadamente.

3.5 Utilice aritmética de 5 dígitos con corte para determinar las raíces de la
ecuación siguiente.
F(X)=x2 – 5 000.002x + 10

Calcule los errores relativos porcentuales de sus resultados.

x x+1.5
99999 100000.5 100002 100003.5 ..
x2 – 5 000.002x + 10 9499804811 9500097310 9500389814 9500682322 …
0% 0% 0% …
Error% Error% Error% Error% …

3. 6 La función exponencial y = e ^x está dada por la serie infinita

Evalúe esta función para x = 10 y x = –10; esté atento al problema del error de redondeo.

y=e^-x y=e ^x x
0.367879441 2.718281828 1
0.135335283 7.389056099 2
0.049787068 20.08553692 3
0.018315639 54.59815003 4
0.006737947 148.4131591 5
0.002478752 403.4287935 6
0.000911882 1096.633158 7
0.000335463 2980.957987 8
0.00012341 8103.083928 9
4.53999E-05 22026.46579 10

CAPÍTULO 4
4.1 La serie infinita se utiliza para aproximar e^x:

IT RESULTADO ET(%) EA(%)

1 1.000000000000 18.127%
2 1.000000000000 18.127% 0.0000000%
3 0.500000000000 59.063% -100.0000000%
4 0.166666666667 86.354% -200.0000000%
5 0.041666666667 96.589% -300.0000000%
6 0.008333333333 99.318% -400.0000000%
7 0.001388888889 99.886% -500.0000000%
8 0.000198412698 99.984% -600.0000000%
9 0.000024801587 99.998% -700.0000000%
10 0.000002755732 100.000% -800.0000000%
11 0.000000275573 100.000% -900.0000000%
12 0.000000025052 100.000% -1000.0000000%
13 0.000000002088 100.000% -1100.0000000%
14 0.000000000161 100.000% -1200.0000000%
15 0.000000000011 100.000% -1300.0000000%
Valor verdadero para e^2 = 1.221402758 en x=1

4.2 La expansión en serie de Maclaurin para cos x es

Iniciando con el primer término cos x = 1, agregue los términos uno a uno para
estimar cos(p/3). Después de que agregue cada uno de los términos, calcule los
errores relativos porcentuales verdaderos y aproximados. Use una calculadora
para determinar el valor exacto. Agregue términos hasta que el valor absoluto
del error aproximado se encuentre dentro de cierto criterio de error,
considerando dos cifras significativas.

IT RESULTADO ET(%) EA(%)

1 1.000000000000 -100.000%
2 0.500000000000 0.000% -100.0000000%
3 0.541666666667 -8.333% 7.6923077%
4 0.540277777778 -8.056% -0.2570694%
5 0.540302579365 -8.061% 0.0045903%
6 0.540302303792 -8.060% -0.0000510%
7 0.540302305880 -8.060% 0.0000004%
8 0.540302305868 -8.060% 0.0000000%
9 0.540302305868 -8.060% 0.0000000%
10 0.540302305868 -8.060% 0.0000000%
11 0.540302305868 -8.060% 0.0000000%
12 0.540302305868 -8.060% 0.0000000%
13 0.540302305868 -8.060% 0.0000000%
14 0.540302305868 -8.060% 0.0000000%
15 0.540302305868 -8.060% 0.0000000%
Valor verdadero para cos (pi/3) = 0.5 en x=1
4.3 Repita los cálculos del problema 4.2, pero ahora usando la expansión de la
serie de Maclaurin para sen x, para evaluar el sen(p/3).

IT RESULTADO ET(%) EA(%)

1 1.000000000000 -100.000%
2 0.833333333333 -66.667% -20.0000000%
3 0.841666666667 -68.333% 0.9900990%
4 0.841468253968 -68.294% -0.0235793%
5 0.841471009700 -68.294% 0.0003275%
6 0.841470984648 -68.294% -0.0000030%
7 0.841470984809 -68.294% 0.0000000%
8 0.841470984808 -68.294% 0.0000000%
9 0.841470984808 -68.294% 0.0000000%
10 0.841470984808 -68.294% 0.0000000%
11 0.841470984808 -68.294% 0.0000000%
12 0.841470984808 -68.294% 0.0000000%
13 0.841470984808 -68.294% 0.0000000%
14 0.841470984808 -68.294% 0.0000000%
15 0.841470984808 -68.294% 0.0000000%
Valor verdadero para sen (pi/3) = 0.866025404 en x=1

CAPÍTULO 5
5.1 Determine las raíces reales de f (x) = −0.5x2 + 2.5x + 4.5:
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 xl = 5 y xu = 10. Calcule el error
estimado ea y el error verdadero e para cada iteración.
5.2 Determine las raíces reales de f (x) = 5x3 − 5x2 + 6x − 2:
a) Gráficamente.

b) Utilizando el método de bisección para localizar la raíz. Use los valores iniciales xl =
0 y xu = 1 iterando hasta que el error estimado ea se encuentre debajo de es = 10%.

5.6 Determine la raíz real de ln (x2) = 0.7:

a) Gráficamente.
b) Empleando tres iteraciones en el método de bisección con los valores iniciales xl =
0.5 y xu = 2.

c) Usando tres iteraciones del método de la falsa posición, con los mismos valores
iniciales de b).
5.3 Determine las raíces reales de f (x) = −25 + 82x − 90x2 + 44x3 − 8x4 + 0.7x5:
a) Gráficamente.

b) Usando el método de bisección para localizar la raíz más grande con e = 10%.
Utilice como valores iniciales xl = 0.5 y xu = 1.0.
c) Realice el mismo cálculo que en b), pero con el método de la falsa posición y e =
0.2%.

CAPÍTULO 6
6.1 Utilice la iteración simple de punto fijo para localizar la raíz de f(x) = sen (x)–
x
Haga una elección inicial de x0 = 0.5 e itere hasta que E ABS ≤ 0.01%.
 n x0 x1 error E%
1 0.5 0.479425539
2 0.47942554 0.461269555 0.0393609 4%
3 0.46126956 0.445085337 0.03636206 4%
4 0.44508534 0.430534905 0.03379617 3%
5 0.4305349 0.417356953 0.03157477 3%
6 0.41735695 0.405345694 0.02963214 3%
7 0.40534569 0.394336465 0.02791836 3%
8 0.39433647 0.384195664 0.02639489 3%
9 0.38419566 0.374813558 0.02503139 3%
10 0.37481356 0.366099037 0.02380373 2%
11 0.36609904 0.357975715 0.02269238 2%
12 0.35797572 0.350378991 0.02168145 2%
13 0.35037899 0.343253797 0.0207578 2%
14 0.3432538 0.336552853 0.01991052 2%
15 0.33655285 0.330235304 0.01913045 2%
16 0.3302353 0.324265627 0.01840984 2%
17 0.32426563 0.31861277 0.01774209 2%
18 0.31861277 0.313249451 0.01712156 2%
19 0.31324945 0.308151582 0.01654338 2%
20 0.30815158 0.303297805 0.01600333 2%
21 0.30329781 0.298669108 0.01549774 2%
22 0.29866911 0.294248496 0.0150234 2%
23 0.2942485 0.290020727 0.01457747 1%
24 0.29002073 0.285972087 0.01415747 1%
… … … … …
89 0.17209249 0.171244303 0.00495308 0%

6.2 Utilice los métodos de a) Newton-Raphson, para determinar una raíz de

f(x) = -0.9x2 + 1.7x+ 2.5 con el uso de x0 = 5. Haga el cálculo hasta que sea
menor que es = 0.01%.
IT XI FXI F'XI FXI / F'XI EA ER ERP
0 5 -11.5 -7.3 1.57534247
1 3.42465753 -2.2335335 - 0.50030054 1.57534247 0.46 46%
4.46438356
2 2.924357 - - 0.06321002 0.50030054 0.17108053 17%
0.22527057 3.56384259
3 2.86114698 - - 0.00104229 0.06321002 0.02209255 2%
0.00359596 3.45006456
4 2.86010469 -9.7773E-07 - 2.8355E-07 0.00104229 0.00036442 0%
3.44818844
5 2.86010441 -7.1942E-14 - 2.0864E-14 2.8355E-07 9.9139E-08 0%
3.44818793

6.3 Determine las raíces reales de f(x) = -1 + 5.5x - 4x2 + 0.5x3: a) con el método
de Newton-Raphson dentro de ERP = 0.01%.

IT XI FXI F'XI FXI / F'XI EA ER ERP

0 10 154 75.5 2.0397351
1 7.9602649 41.5225334 36.8666067 1.12629116 2.0397351 0.2562396 26%
2 6.83397374 9.35827837 20.8830056 0.4481289 1.12629116 0.16480765 16%
3 6.38584483 1.21031834 15.5817627 0.07767532 0.4481289 0.07017535 7%
4 6.30816952 0.03342491 14.7241479 0.00227007 0.07767532 0.01231345 1%
5 6.30589944 2.8142E-05 14.6993561 1.9145E-06 0.00227007 0.00035999 0%

6.4 Emplee el método de Newton-Raphson para determinar una raíz real de

f(x) = –1 + 5.5x – 4x2 + 0.5x3 con el uso de valores iniciales de:
a) 4.52
IT XI FXI F'XI FXI / F'XI EA ER ERP
0 4.52 -11.688896 -0.0144 811.728889
1 -807.208889 -265593868 983842.457 -269.955689 811.728889 1.00559954 101%
2 -537.2532 -78694170.8 437265.028 -179.969048 269.955689 0.50247386 50%
3 -357.284152 -23316586.1 194341.721 -119.977254 179.969048 0.50371405 50%
4 -237.306898 -6908481.7 86375.8012 -79.9816801 119.977254 0.50557845 51%
5 -157.325218 -2046867.07 38390.9382 -53.3164117 79.9816801 0.50838436 51%
6 -104.008806 -606419.264 17064.3182 -35.5372689 53.3164117 0.5126144 51%
7 -68.4715376 -179640.316 7585.79949 -23.6811316 35.5372689 0.5190079 52%
8 -44.790406 -53200.888 3373.09395 -15.7721335 23.6811316 0.52870991 53%
9 -29.0182725 -15746.4064 1500.73639 -10.4924532 15.7721335 0.5435242 54%
10 -18.5258193 -4654.80193 668.515527 -6.96289278 10.4924532 0.56636919 57%
11 -11.5629265 -1372.39016 298.555318 -4.59677011 6.96289278 0.60217392 60%
12 -6.96615643 -402.447707 134.020255 -3.0028872 4.59677011 0.65987179 66%
13 -3.96326923 -116.754525 60.7674084 -1.9213346 3.0028872 0.75767934 76%
14 -2.04193463 -33.1655487 28.0897226 -1.18070047 1.9213346 0.94093835 94%
15 -0.86123415 -9.02308406 13.5024596 -0.66825484 1.18070047 1.37094014 137%
16 -0.19297931 -2.21394364 7.099696 -0.3118364 0.66825484 3.46283155 346%
17 0.11885709 -0.40195449 4.5703338 -0.08794861 0.3118364 2.62362475 262%
18 0.2068057 -0.02922066 3.90970732 -0.00747387 0.08794861 0.42527169 43%
19 0.21427957 -0.0002059 3.85463704 -5.3416E-05 0.00747387 0.03487908 3%
20 0.21433299 -1.0496E-08 3.85424405 -2.7232E-09 5.3416E-05 0.00024922 0%
 21 0.21433299 -7.2858E-17 3.85424403 -1.8903E-17 2.7232E-09 1.2705E-08 0%

b) 4.54.
IT XI FXI F'XI FXI / F'XI EA ER ERP
0 4.54 -11.688068 0.0974 -120.000698
1 124.540698 904479.551 22274.7527 40.6055935 120.000698 0.96354605 96%
2 83.9351047 267945.914 9901.67186 27.060674 40.6055935 0.48377367 48%
3 56.8744307 79358.891 4402.55585 18.0256409 27.060674 0.47579683 48%
4 38.8487897 23491.6041 1958.55238 11.9943711 18.0256409 0.46399492 46%
5 26.8544186 6945.22366 872.404349 7.96101449 11.9943711 0.44664423 45%
6 18.8934041 2047.17241 389.793846 5.25193619 7.96101449 0.42136475 42%
7 13.6414679 598.937461 175.502727 3.41269603 5.25193619 0.38499788 38%
8 10.2287719 171.853967 80.6114864 2.13187939 3.41269603 0.33363693 33%
9 8.0968925 46.7090288 39.0643622 1.19569413 2.13187939 0.263296 26%
10 6.90119836 10.7905316 21.7302214 0.49656796 1.19569413 0.17325891 17%
11 6.40463041 1.50500275 15.7918927 0.09530224 0.49656796 0.07753265 8%
12 6.30932817 0.05049239 14.7368076 0.00342628 0.09530224 0.01510497 2%
13 6.30590189 6.4124E-05 14.6993829 4.3623E-06 0.00342628 0.00054334 0%
14 6.30589753 1.0388E-10 14.6993352 7.0671E-12 4.3623E-06 6.9179E-07 0%

6.5 Determine la raíz real más pequeña de f(x) = -12 - 21x + 18x^2- 2.4x^3, realice
gráfico.
x f(x)=12+21*x g(x)=19*x^2-2.4*x^3
1 33.0 16.6
1.904 52.0 52.3
2.808 71.0 96.7
3.712 90.0 139.0
4.616 108.9 168.8
5.52 127.9 175.3
6.424 146.9 147.8
7.328 165.9 75.9
8.232 184.9 -51.3
9.136 203.9 -244.3
10.04 222.8 -513.7
10.944 241.8 -870.2
11.848 260.8 -1324.5
 6.6 Localice la primera raíz positiva de f(x) = sen x + cos (1 + x2) – 1 donde x está
en radianes. Para localizar la raíz, use cuatro iteraciones del método de la
secante con valores iniciales de:
a) xi – 1 = 1.0 y xi= 3.0

i xi xi-x(i-1) fxi f(xi-1) E abs

0 1
1 3 2 -1.69795152 -0.57467585 200%
2 -0.02321428 -3.02321428 -20.2310775 55 302%
Respuesta 3 0.78978594 0.81300022 -10.3619786 -20.2310775 81%
pedida
Respuesta 4 1.64338879 0.85360285 6.273685229 -10.3619786 85%
con los
valores
iniciales
5 1.32147584 -0.32191295 -0.98495369 6.27368523 32%
6 1.36515749 0.04368165 -0.07693257 -0.98495369 4%
Respuesta 7 1.36885844 0.00370095 0.001061895 -0.07693257 0%
optima

b) xi – 1 = 1.5 y xi = 2.5

i xi xi-x(i-1) fxi f(xi-1) E abs

0 1.5
1 2.5 1 0.166396317 -0.99663469 100%
2 2.35692873 -0.14307127 27.77251908 33.125 14%
3 1.61457215 -0.74235658 5.568344539 27.7725191 74%
Respuesta 4 1.42840458 -0.18616757 1.279155685 5.56834454 19%
con los
valores
iniciales
 5 1.37288423 -0.05552036 0.086091876 1.27915569 6%
Respuesta 6 1.36887786 -0.00400637 0.001471498 0.08609188 0%
optima

c) xi – 1 = 1.5 y xi = 2.25.

i xi xi-x(i-1) fxi f(xi-1) E abs

0 1.5
1 2.25 0.75 0.753820863 -0.99663469 75%
2 1.92701799 -0.32298201 13.85276205 24.015625 32%
3 1.48676873 -0.44024926 2.575124305 13.8527621 44%
Respuesta con 4 1.38624267 -0.10052607 0.369663312 2.57512431 10%
los valores
iniciales
5 1.36939322 -0.01684945 0.012345682 0.36966331 2%
Respuesta 6 1.36881105 -0.00058217 6.21514E-05 0.01234568 0%
optima

CAPÍTULO 7
7.1 Sea f(x)= x5- 3.5x4 + 2.75x3 + 2.125x2 - 3.875x + 1.25

Tenemos que almacenar los coeficientes del polinomio en el vector fila p= [1 -3.5 2.75 2.125 -3.875
1.25];

Mediante la función plolyval, podemos calcular el valor del polinomio cuando proporcionamos el
valor de x

código en Matlab:
7.2 Consideremos el siguiente mecanismo formado por tres barras articuladas.
Conocidas las longitudes de las barras r2, r3 y r4 y la distancia r1 entre la primera
articulación y la última queremos relacionar los ángulos α y φ. Como vemos en el
dibujo, el vector →r1=→r2+→r3+→r4r→1=r→2+r→3+r→4 es la suma de los otros
tres.

r2cosφ+r3cosβ−r4cosα=r1
r2sinφ=r3sinβ+r4sinα
Eliminamos el ángulo β, de la forma en q se muestra:

Conocidos r1, r2, r3 y r4, dado el ángulo α se obtiene el ángulo φ

Resultado 120.00 y 86.1015

7.3 Deflación polinomial planteamiento del problema, divida el polinomio de

segundo grado f(x) = (x – 4) (x + 6) = x2 + 2x – 24 entre el factor x – 4.

Solución

Usando el método propuesto en el pseudocódigo anterior, los parámetros son n = 2,

a0= –24, a1= 2, a2 = 1 y t = 4. Estos valores se usan para calcular

r = a2 = 1 a2 = 0

El loop o ciclo se itera después desde i = 2 – 1 = 1 hasta 0. Para i = 1,

 s = a1 = 2 → a1 = r = 1 → r = s + rt = 2 + 1(4) = 6 → Para i = 0 → s = a0 = –24 →
a0 = r = 6 → r = –24 + 6(4) = 0

resultando el cociente a0 + a1x = 6 + x, con un residuo de cero.

CAPÍTULO 8
8.1 Método de la bisección sea f(x)= (2.0+(12.02/(x^2)))*(x-0.08407)-0.082054*375,
con un intervalo de [10, 20] y un error de 0.001.
IT XA XB XR F(XA) F(XR) F(XA)*F(XR) EA(%) ET(%)
1 10 20 15 - - 1.37959567 1400.0000000
9.74649521 0.14154787
2 15 20 17.5 - 4.74516748 -0.67166836 14.2857143 1650.0000000
0.14154787
3 15 17.5 16.25 - 2.29747548 -0.32520277 7.6923077 1525.0000000
0.14154787
4 15 16.25 15.625 - 1.0767509 -0.1524118 4.0000000 1462.5000000
0.14154787
5 15 15.625 15.3125 - 0.46727983 -0.06614247 2.0408163 1431.2500000
0.14154787
6 15 15.3125 15.15625 - 0.16278308 -0.0230416 1.0309278 1415.6250000
0.14154787
7 15 15.15625 15.078125 - 0.01059656 -0.00149992 0.5181347 1407.8125000
0.14154787
8 15 15.078125 15.0390625 - - 0.00926869 0.2597403 1403.9062500
0.14154787 0.06548096
9 15.0390625 15.078125 15.0585938 - - 0.00179703 0.1297017 1405.8593750
0.06548096 0.02744352
10 15.0585938 15.078125 15.0683594 - - 0.00023118 0.0648088 1406.8359375
0.02744352 0.00842381
11 15.0683594 15.078125 15.0732422 - 0.00108629 -9.1507E-06 0.0323939 1407.3242188
0.00842381
12 15.0683594 15.0732422 15.0708008 - - 3.0905E-05 0.0161996 1407.0800781
0.00842381 0.00366878
13 15.0708008 15.0732422 15.0720215 - - 4.7373E-06 0.0080991 1407.2021484
0.00366878 0.00129125
14 15.0720215 15.0732422 15.0726318 - - 1.3232E-07 0.0040494 1407.2631836
0.00129125 0.00010248
15 15.0726318 15.0732422 15.072937 - 0.00049191 -5.0409E-08 0.0020247 1407.2937012
0.00010248
16 15.0726318 15.072937 15.0727844 - 0.00019472 -1.9954E-08 0.0010123 1407.2784424
0.00010248
17 15.0726318 15.0727844 15.0727081 - 4.6119E-05 -4.7261E-09 0.0005062 1407.2708130
0.00010248

8.2 Método de la bisección sea f(x)= ((R+1)/(R(1+0.1R))-Ln((1+0.1R)/(0.1R)) , con

un intervalo de [0.01, 2] y un error de 0.001.
IT XA XB XR F(XA) F(XR) F(XA)*F(XR) EA(%) ET(%)
1 0.01 2 1.005 93.9903461 - -54.5639535 0.5000000
0.58052721
2 0.01 1.005 0.5075 93.9903461 - -19.1151416 98.0295567 49.2500000
0.20337346
3 0.01 0.5075 0.25875 93.9903461 1.06201023 99.8187092 96.1352657 74.1250000
4 0.25875 0.5075 0.383125 1.06201023 0.1773294 0.18832563 32.4632953 61.6875000
5 0.383125 0.5075 0.4453125 0.1773294 - -0.00849191 13.9649123 55.4687500
0.04788776
6 0.383125 0.4453125 0.41421875 0.1773294 0.05385329 0.00954977 7.5066013 58.5781250
 7 0.41421875 0.4453125 0.42976563 0.05385329 0.00058547 3.153E-05 3.6175241 57.0234375
8 0.42976563 0.4453125 0.43753906 0.00058547 - -1.4178E-05 1.7766271 56.2460938
0.02421568
9 0.42976563 0.43753906 0.43365234 0.00058547 - -7.0025E-06 0.8962753 56.6347656
0.01196038
10 0.42976563 0.43365234 0.43170898 0.00058547 - -3.3514E-06 0.4501550 56.8291016
0.00572431
11 0.42976563 0.43170898 0.4307373 0.00058547 -0.0025787 -1.5098E-06 0.2255852 56.9262695
12 0.42976563 0.4307373 0.43025146 0.00058547 - -5.8485E-07 0.1129200 56.9748535
0.00099894
13 0.42976563 0.43025146 0.43000854 0.00058547 - -1.2138E-07 0.0564919 56.9991455
0.00020732
14 0.42976563 0.43000854 0.42988708 0.00058547 0.00018893 1.1061E-07 0.0282539 57.0112915
15 0.42988708 0.43000854 0.42994781 0.00018893 -9.2292E- -1.7437E-09 0.0141250 57.0052185
06
16 0.42988708 0.42994781 0.42991745 0.00018893 8.9842E-05 1.6974E-08 0.0070630 57.0082550
17 0.42991745 0.42994781 0.42993263 8.9842E-05 4.0304E-05 3.621E-09 0.0035314 57.0067368
18 0.42993263 0.42994781 0.42994022 4.0304E-05 1.5537E-05 6.2621E-10 0.0017657 57.0059776
19 0.42994022 0.42994781 0.42994402 1.5537E-05 3.1538E-06 4.9E-11 0.0008828 57.0055981

8.3 Método de la bisección sea f(x)= 0.05-((x/(1-x))*(raiz(6/(2+x)))), con un

intervalo de [0.01, 0.05] y un error de 0.001.
IT XA XB XR F(XA) F(XR) F(XA)*F(XR) EA(%) ET(%)
1 0.01 0.05 0.03 0.03254811 - -0.00010322 97.0000000
0.00317128
2 0.01 0.03 0.02 0.03254811 0.01482745 0.00048261 50.0000000 98.0000000
3 0.02 0.03 0.025 0.01482745 0.00586344 8.694E-05 20.0000000 97.5000000
4 0.025 0.03 0.0275 0.00586344 0.00135499 7.9449E-06 9.0909091 97.2500000
5 0.0275 0.03 0.02875 0.00135499 - -1.2275E-06 4.3478261 97.1250000
0.00090591
6 0.0275 0.02875 0.028125 0.00135499 0.0002251 3.0501E-07 2.2222222 97.1875000
7 0.028125 0.02875 0.0284375 0.0002251 - -7.6594E-08 1.0989011 97.1562500
0.00034026
8 0.028125 0.0284375 0.02828125 0.0002251 -5.7546E- -1.2954E-08 0.5524862 97.1718750
05
9 0.028125 0.02828125 0.02820313 0.0002251 8.3786E-05 1.886E-08 0.2770083 97.1796875
10 0.02820313 0.02828125 0.02824219 8.3786E-05 1.3122E-05 1.0995E-09 0.1383126 97.1757813
11 0.02824219 0.02828125 0.02826172 1.3122E-05 -2.2211E- -2.9146E-10 0.0691085 97.1738281
05
12 0.02824219 0.02826172 0.02825195 1.3122E-05 -4.5444E- -5.9633E-11 0.0345662 97.1748047
06
13 0.02824219 0.02825195 0.02824707 1.3122E-05 4.289E-06 5.6282E-11 0.0172861 97.1752930
14 0.02824707 0.02825195 0.02824951 4.289E-06 -1.2767E- -5.4756E-13 0.0086423 97.1750488
07
15 0.02824707 0.02824951 0.02824829 4.289E-06 2.0807E-06 8.924E-12 0.0043213 97.1751709
16 0.02824829 0.02824951 0.0282489 2.0807E-06 9.765E-07 2.0318E-12 0.0021606 97.1751099
17 0.0282489 0.02824951 0.02824921 9.765E-07 4.2442E-07 4.1445E-13 0.0010803 97.1750793
18 0.02824921 0.02824951 0.02824936 4.2442E-07 1.4838E-07 6.2974E-14 0.0005401 97.1750641

8.4 Método de NEWTON RAPHSON sea f(x)= 10(1-(e^(-0.04*x)))+4*(e^(-0.04*x))-

0.93(10), con un error de 0.001.
f '(x) = 6/(25*(e^(0.04*x)))

IT Xi f(Xi) f'(Xi) f(Xi)/f'(Xi) Ea(%)

0 1 -5.06473663 0.230589465 -21.9643019
1 22.9643019 -1.69453101 0.095781241 -17.6916795 95.6454152
2 40.6559814 -0.48000681 0.047200272 -10.1695771 43.5155638
3 50.8255586 -0.08563508 0.031425403 -2.72502713 20.0087858
4 53.5505857 -0.00450211 0.028180084 -0.15976194 5.08869716
5 53.7103477 -1.4355E-05 0.028000574 -0.00051266 0.29745096
6 53.7108603 -1.4718E-10 0.028 -5.2564E-09 0.00095448

8.5 Método de SECANTE sea f(x)= ((4+x)/((42-2x)^2)*(30-x))-0.015, con un error de

0.0001.

IT XI F(XI) F(XI-1) XI-(XI-1) EA(%)

0 15
1 20 0.585 -0.0062037 5
2 15.0524667 -0.00599159 0.585 -4.94753328 32.8685881
3 15.1026258 -0.00578265 -0.00599159 0.05015911 0.33212179
4 16.4908198 0.00364984 -0.00578265 1.38819401 8.41798059
5 15.9536672 -0.00105405 0.00364984 -0.53715269 3.36695434
6 16.0740328 -0.00014864 -0.00105405 0.12036569 0.74882071
7 16.0937929 7.2629E-06 -0.00014864 0.01976008 0.12278078
8 16.0928724 -4.7339E-08 7.2629E-06 -0.00092055 0.00572021
9 16.0928783 -1.4964E-11 -4.7339E-08 5.9612E-06 3.7043E-05

8.6 Método de SECANTE sea f(x)= ((4+x)/((42-2x)^2)*(30-x))-0.015, con un error de

0.0001

IT XI F(XI) F(XI-1) XI-(XI-1) EA(%)

0 15
1 20 0.585 -0.0062037 5
2 15.0524667 -0.00599159 0.585 -4.94753328 32.8685881
3 15.1026258 -0.00578265 -0.00599159 0.05015911 0.33212179
4 16.4908198 0.00364984 -0.00578265 1.38819401 8.41798059
5 15.9536672 -0.00105405 0.00364984 -0.53715269 3.36695434
6 16.0740328 -0.00014864 -0.00105405 0.12036569 0.74882071
7 16.0937929 7.2629E-06 -0.00014864 0.01976008 0.12278078
8 16.0928724 -4.7339E-08 7.2629E-06 -0.00092055 0.00572021
9 16.0928783 -1.4964E-11 -4.7339E-08 5.9612E-06 3.7043E-05

CAPÍTULO 9
9.1 Método de Gauss-Jordan, Con la técnica de Gauss-Jordan resuelva el sistema
del ejemplo 9.5:
3x1 – 6x2 – x3 = -38
3x1 -1x2 +7x3 = –34
8x1 +x2 -2x3 = -20
Código y desarrollo en Matlab:
9.2 Utilice la regla de Cramer para resolver
0.3x1 + 0.52x2 + x3 = –0.01
0.5x1 + x2 + 1.9x3 = 0.67
0.1x1 + 0.3x2 + 0.5x3 = –0.44
Código y desarrollo en Matlab:
9.3 Método de Cramer

Planteamiento del problema Con la técnica de Cramer resuelva el sistema del

Sea el sistema

3x1−x2=5

−2x1+x2+x3=0

2x1−x2+4x3=15

Código y desarrollo en Matlab:

 CAPÍTULO 10
10.1 Código y Desarrollo en MATLAB:

%grafica funcion
clc
syms x1 x2 x3;
F1=10*x1+2*x2-1*x3-27;
F2=-3*x1-6*x2+2*x3+61.5;
F3=1*x1+1*x2+5*x3+21.5;
%calcular las dterminates
A=[ 10 2 -1 ; -3 -6 2 ;1 1 5];
disp(A);
D=det(A);
disp(D);
x1=[ 27 -1 2 ; -61.5 2 -6 ;-21.5 5 1];
Dx=det(x1);
x1=Dx/D;
disp ('El valor de X1 es: ');
disp(x1);
x2=[ 10 -1 27 ; -3 2 -61.5 ;1 5 -21.5];
Dy=det(x2);
x2=Dy/D;
disp ('El valor de x2 es: ');
disp(x2);
x3=[ 10 2 27 ; -3 -6 -61.5 ;1 1 -21.5];
Dz=det(x3);
x3=Dz/D;
disp ('El valor de x3 es: ');
disp(x3)
10.2 Código y Desarrollo en MATLAB:

clc
syms x1 x2 x3;
F1=8*x1+4*x2-1*x3-11;
F2=-2*x1+5*x2+1*x3-4;
F3=2*x1-1*x2+6*x3-7;
%calcular las dterminates
A=[ 8 4 -1 ; -2 5 1 ;2 -1 6];
disp(A);
D=det(A);
disp(D);
x1=[ 11 -1 4 ; 4 1 5 ;7 6 -1];
Dx=det(x1);
x1=Dx/D;
disp ('El valor de X1 es: ');
disp(x1);
x2=[ 8 -1 11 ; -2 1 4 ;2 6 7];
Dy=det(x2);
x2=Dy/D;
disp ('El valor de x2 es: ');
disp(x2);
x3=[ 8 4 11 ; -2 5 4 ;2 -1 7];
Dz=det(x3);
x3=Dz/D;
disp ('El valor de x3 es: ');
disp(x3)

10.3 Código y Desarrollo en MATLAB:

clc
syms x1 x2 x3;
F1=2*x1-6*x2-1*x3+38;
F2=-3*x1-1*x2+7*x3+34;
F3=-8*x1+1*x2-2*x3+20;
%calcular las dterminates
A=[ 2 -6 -1 ; -3 -1 7 ;-8 1 -2];
disp(A);
D=det(A);
disp(D);
x1=[ -38 -1 -6 ; -34 7 -1 ;-20 -2 1];
Dx=det(x1);
x1=Dx/D;
disp ('El valor de X1 es: ');
disp(x1);
x2=[ 2 -1 -38 ; -3 7 -34 ;-8 -2 -20];
Dy=det(x2);
x2=Dy/D;
disp ('El valor de x2 es: ');
disp(x2);
x3=[ 2 -6 -38 ; -3 -1 -34 ;-8 1 -20];
Dz=det(x3);
x3=Dz/D;
disp ('El valor de x3 es: ');
disp(x3)

10.4 Código y Desarrollo en MATLAB:

%grafica funcion
clc
syms x1 x2 x3;
F1=3*x1-2*x2+1*x3+10;
F2=2*x1+6*x2-4*x3-44;
F3=-1*x1-2*x2+5*x3+26;
%calcular las dterminates
A=[ 3 -2 1 ; 2 6 -4 ;1 -2 5];
disp(A);
D=det(A);
disp(D);
x1=[ -10 1 -2 ; 44 -4 6 ;-26 5 -2];
Dx=det(x1);
x1=Dx/D;
disp ('El valor de X1 es: ');
disp(x1);
x2=[ 3 1 -10 ; 2 -4 44 ;-1 5 -26];
Dy=det(x2);
x2=Dy/D;
disp ('El valor de x2 es: ');
disp(x2);
x3=[ 3 -2 -10 ; 2 6 44 ;-1 -2 -26];
Dz=det(x3);
x3=Dz/D;
disp ('El valor de x3 es: ');
disp(x3)

CAPÍTULO 11
11.1 Emplee el método de Gauss-Seidel para resolver el sistema siguiente hasta
que el error relativo porcentual esté por debajo de ԑ= 5%.
10x1 + 2x2 – x3 = 27
–3x1 – 6x2 + 2x3 = –61.5
x1 + x2 + 5x3 = –21.5
Código y Desarrollo en MATLAB:

11.2 Utilice el método de Gauss-Seidel, para resolver el sistema siguiente para una
tolerancia de ԑ = 5%. Si es necesario, reacomode las ecuaciones para lograr
convergencia.
–3x1 + x2 – 12x3 = 50
6x1 – x2 – x3 = 3
6x1 + 9x2 + x3 = 40
Código y Desarrollo en MATLAB:
11.3 Use el método de Gauss-Seidel, para resolver el sistema siguiente para una
tolerancia de ԑ = 5%. Si es necesario, reacomode las ecuaciones para lograr
convergencia.
2x1 – 6x2 – x3 = –38
–3x1 – x2 + 7x3 = –34
–8x1 + x2 – 2x3 = –20
Código y Desarrollo en MATLAB:
11.4 Determinar concentraciones (las c están en g/m3) en una serie de reactores
acoplados como función de la cantidad de masa de entrada a cada uno de ellos
(los lados derechos están en g/d).
15c1 – 3c2 – c3 = 3 800
–3c1 + 18c2 – 6c3 = 1 200
–4c1 – c2 + 12c3 = 2 350
Código y Desarrollo en MATLAB:
 CAPÍTULO 12
12.1 un circuito con resistores para resolverse usando ecuaciones algebraicas
lineales simultaneas.

SOLUCIÓN:
la regla de la corriente de Kirchhoff se aplica a cada nodo para obtener:
i12 + i52 + i32 = 0
i65 – i52 – i54 = 0
i43 – i32 = 0
i54 – i43 = 0
La aplicación de la regla de voltajes en cada una de las mallas da:
–i54R54 – i43R43 – i32R32 + i52R52 = 0
–i65R65 – i52R52 + i12R12 – 200 = 0
o, sustituyendo el valor de las resistencias de la figura 1 y pasando las constantes al
lado derecho,
–l5i54 – 5i43 – l0i32 + l0i52 = 0
–20i65 – 10i52 + 5i12 = 200
 Por lo tanto, el problema consiste en la solución del siguiente conjunto de seis
ecuaciones con seis corrientes como incógnitas:

Aunque no es práctico resolverlo a mano, este sistema se resuelve de manera sencilla

con un método de eliminación. Si se procede de esta forma, la solución es:
i12 = 6.1538 i52 = –4.6154 i32 = –1.5385
i65 = –6.1538 i54 = –1.5385 i43 = –1.5385
12.2 Las fuerzas (F) representan ya sea la tensión o la compresión sobre los
componentes de la armadura. Las reacciones externas (H2, V2 y V3) son fuerzas
que caracterizan cómo interactúa dicha estructura con la superficie de soporte. El
apoyo fijo en el nodo 2 puede transmitir fuerzas horizontales y verticales a la
superficie, mientras que el apoyo móvil en el nodo 3 transmite sólo fuerzas
verticales. Se observa que el efecto de la carga externa de 1 000 lb se distribuye
entre los componentes de la armadura.
SOLUCIÓN:
Este tipo de estructura se puede describir como un conjunto de ecuaciones algebraicas
lineales acopladas. Los diagramas de fuerza de cuerpo libre para cada nodo se
muestran en la figura 2. La suma de las fuerzas en ambas direcciones, vertical y
horizontal, deben ser cero en cada nodo, ya que el sistema está en reposo. Por lo tanto,
para el nodo 1,
ΣFH = 0 = –F1 cos 30° + F3 cos 60° + F1, h
ΣFV = 0 = –F1 sen 30° – F3 sen 60° + F1, v
PARA EL NODO 2,
ΣFH = 0 = F2 + F1 cos 30° + F2, h + H2
ΣFV = 0 = F1 sen 30° + F2, v + V2
PARA EL NODO 3,
ΣFH = 0 = –F2 – F3 cos 60° + F3, h
Solución del siguiente sistema de seis ecuaciones con seis incógnitas:

Sin embargo, como este problema es un caso de estudio ideal, para demostrar la utilidad
de la matriz inversa se utiliza la descomposición LU para calcular:
F1 = –500 F2 = 433 F3 = –866
H2 = 0 V2 = 250 V3 = 750
la matriz inversa es:
Debido a que las fuerzas externas no tienen efecto sobre la descomposición LU, no se
necesita aplicar el método una y otra vez para estudiar el efecto de diferentes fuerzas

Si la fuerza del viento se puede idealizar como dos fuerzas puntuales de 1 000 libras
sobre
los nodos 1 y 2 (figura a), el vector del lado derecho es:
{F}T = ⎣–1 000 0 1 000 0 0 0⎦
que se utiliza para calcular:
F1 = –866 F2 = 250 F3 = –500
H2 = –2 000 V2 = –433 V3 = 433
Para un viento de la derecha (figura 12.2b), F1, h = –1 000, F3, h = –1 000, y todas las
demás fuerzas externas son cero, con lo cual resulta:
F1 = –866 F2 = –1 250 F3 = 500
H2 = 2 000 V2 = 433 V3 = –433
12.3 Use la descomposición LU con el propósito de obtener los desplazamientos
y generar la inversa de [K]. Si m1 = 2 kg, m2 = 3 kg, m3 = 2.5 kg, y todas las k = 10
kg/s2.
SOLUCIÓN:
Sustituyendo los parámetros del modelo se obtiene:
La descomposición LU se utiliza con el objetivo de obtener x1 = 7.35, x2 = 10.045 y x3
= 12.495. Estos desplazamientos se utilizaron para construir la figura 12.11b. La inversa
de la matriz de rigidez calculada es:

Cada elemento de la matriz k–1 ji nos indica el desplazamiento de la masa i debido a

una fuerza unitaria impuesta sobre la masa j. Así, los valores 0.1 en la columna 1 nos
indican que una carga unitaria hacia abajo en la primera masa desplazará todas las
masas 0.1 m hacia abajo. Los otros elementos se interpretan en forma similar. Por lo
tanto, la inversa de la matriz de rigidez proporciona una síntesis de cómo los
componentes del sistema responden a fuerzas que se aplican en forma externa.

EJERCICIOS PROPUESTOS
12.3 Se muestra un problema donde donde se aprecia que hay cinco reactores
interconectados o acoplados, se necesitan cinco ecuaciones de balance de masa para
caracterizar el sistema.

ENTRADAS = SALIDAS
HALLANDO ECUACIONES SUSTITUYENDO
1 Q01*C01+Q31*C3=Q15*C1+Q12*C1 6C1-C3=50
2 Q12*C1=Q25*C2+Q24*C2+Q23C3 -3C1+3C2=0
3 Q03*C03+Q23*C2=Q31*C3+Q34*C3 -C2+9C3=160
4 Q54*C5+Q24*C2+Q34*C3=Q44*C4 -C2-8C3+11C4-2C5=0
5 Q15*C1+Q25*C2=Q55*C5+Q54*C5 -3C1-C2+4C5=0
 A X B

6 0 -1 0 0 C1 50
-3 3 0 0 0 C2 0
0 -1 9 0 0 C3 160
0 -1 -8 11 -2 C4 0
-3 -1 0 0 4 C5 0

C
[A]-1 mg/m³
0.169811 0.006289 0.018868 0.000000 0.000000 11.51 c1
0.169811 0.339623 0.018868 0.000000 0.000000 11.51 c2
0.018868 0.037736 0.113208 0.000000 0.000000 19.06 c3
0.060034 0.074614 0.087479 0.090909 0.045455 17.00 c4
0.169811 0.089623 0.018868 0.000000 0.250000 11.51 c5

12.4 Vuelva a calcular las concentraciones para los cinco reactores que se muestran en
la figura 12.3, si los flujos cambian
como sigue:
 ENTRADAS = SALIDAS
HALLANDO ECUACIONES SUSTITUYENDO
1 Q01*C01+Q31*C3=Q15*C1+Q12*C1 8C1-3C3=50
2 Q12*C1=Q25*C2+Q24*C2+Q23C3 -2C1+2C2=0
3 Q03*C03+Q23*C2=Q31*C3+Q34*C3 -2C2+10C3=160
4 Q54*C5+Q24*C2+Q34*C3=Q44*C4 -7C3+10C4-3C5=0
5 Q15*C1+Q25*C2=Q55*C5+Q54*C5 -4C1-2C2+6C5=0

A X B
8 0 -3 0 0 C1 50
-2 2 0 0 0 C2 0
0 -2 10 0 0 C3 160
0 0 -7 10 -3 C4 0
-4 -2 0 0 6 C5 0

C
[A]-1 mg/m³
0.1351 0.0405 0.0405 0.0000 0.0000 13.24 c1
0.1351 0.5405 0.0405 0.0000 0.0000 13.24 c2
0.0270 0.1081 0.1081 0.0000 0.0000 18.65 c3
0.0595 0.1378 0.0878 0.1000 0.0500 17.03 c4
0.1351 0.2072 0.0405 0.0000 0.1667 13.24 c5

12.6 En la figura P12.6 se muestran tres reactores conectados por tubos. Como se
indica, la tasa de transferencia de productos químicos a través de cada tubo es igual a
la tasa de flujo (Q, en unidades de metros cúbicos por segundo) multiplicada por la
concentración del reactor desde el que se origina el flujo (c, en unidades de miligramos
por metro cúbico). Si el sistema se encuentra en estado estacionario (estable), la
transferencia de entra da a cada reactor balanceará la de salida. Desarrolle las
ecuaciones del balance de masa para los reactores y resuelva las tres ecuaciones
algebraicas lineales simultáneas para sus concentraciones.
 ENTRADAS = SALIDAS
HALLANDO ECUACIONES SUSTITUYENDO
1 400+Q21*C2=Q13*C1+Q12*C1 120C1-20C2=400
2 Q12*C1=Q21*C2+Q23*C2 -80C1+80C2=0
3 Q13*C1+Q23*C2+200=Q33*C3 -40C1-60C2+120C3=200

A X B
120 -20 0 C1 400
-80 80 0 C2 0
-40 -60 120 C3 200

C
[A]-1 mg/m³
0.01000 0.00250 0.00000 4.00 c1
0.01000 0.01500 0.00000 4.00 c2
0.00833 0.00833 0.00833 5.00 c3

12.7 Con el empleo del mismo enfoque que en la sección 12.1, determine la
concentración de cloruro en cada uno de los Grandes Lagos con el uso de la información
que se muestra en la figura P12.7.
 ENTRADAS = SALIDAS
HALLANDO ECUACIONES SUSTITUYENDO
superior 180=QSH*CS 67CS=180
michigan 710=QMH*CM 36CM=710
hurón QSH*CS+QMH*CM+740=QHE*CH -67CS-36CM+161CH=740
erie QHE*CH+3850=QEO*CE -161CH+182CE=3850
ontario QEO*CE+4720=QOO*CO -182CE+212CO=4720

A X B
67 0 0 0 0 SUP. 180
0 36 0 0 0 MICH. 710
-67 -36 161 0 0 HUR. 740
0 0 -161 182 0 ERIE 3850
0 0 0 -182 212 ONT. 4720

C
mg/m
[A]-1 ³
0.014925 0.000000 0.000000 0.000000 0.000000
4 0 0 0 0 2.687 SUP.
0.000000 0.027777 0.000000 0.000000 0.000000
0 8 0 0 0 19.722 MICH.
0.006211 0.006211 0.006211 0.000000 0.000000
2 2 2 0 0 10.124 HUR.
0.005494 0.005494 0.005494 0.005494 0.000000
5 5 5 5 0 30.110 ERIE
0.004717 0.004717 0.004717 0.004717 0.004717
0 0 0 0 0 48.113 ONT.
 CAPÍTULO 13
13.1 Dada la fórmula f(x) = –x2 + 8x – 12
a) Determine en forma analítica (esto es, por medio de derivación)
el valor máximo y el correspondiente de x para esta función.
𝑑𝑓
= −2𝑥 + 8 = 0
𝑑𝑥
2𝑥 = 8
𝑥=4
𝑑 2𝑓
= −2
(𝑑𝑥 )2
Debido a que la segunda derivada es negativa, es en ese punto donde existe un
máximo (4)

b) Verifique que la ecuación (13.7) produce los mismos resultados

con base en los valores iniciales de x0 = 0, x1 = 2 y x2 = 6.

X0=0, x1= 2 y x2=6

Calculamos el valor de la función en esos puntos. F(x)=x2-8x+12

F(x0) = (0)2-8(0) +12=12

F(x1) = (2)2-8(2) +12=0
F(x2) = (6)2-8(6) +12=0
Luego reemplazamos en la ecuación 13.7 y nos da:
12(4−36)+0(36−0)+0(0−4)
X3=2(12)(2−6)+2(0)(6−0)+2(0)(0−2)

X3=4
Con esta ecuación se llega al mismo resultado que es 4.
13.2 Dada la función f(x) = –1.5x6 – 2x4 + 12x

a) Grafique la función.
b) Utilice métodos analíticos para probar que la función es cóncava
para todos los valores de x.

f’(x)=-9x5 -8x3 +12

f’’(x)=-45x4-24x2
como la segunda derivada siempre dará negativo para cualquier valor de x, se
concluye que la función es cóncava.

c) Derive la función y después use algún método de localización

de raíces para resolver cuál es el máximo f(x) y el valor correspondiente
de x.
f´(x)= -9x5 -8x3 +12=0
al resolver esta ecuación nos sale que el valor de x = 0.916915
f´´(x)= -45x4-24x2
como ya vimos, esta ecuación será negativa, lo que quiere decir que hay un máximo
en el punto x = 0.916915.
f (0.916915) = -1.5(0.916915)6-2(0.916915)4+12(0.916915)
f (0.916915) = 8.6979
es decir que el máximo valor de la función es 8.6979 en el punto 0.916915.

13.3 Considere la función siguiente: f(x) = –x4 – 2x3 – 8x2 – 5x

Use los métodos analítico y gráfico para demostrar que la función
tiene un máximo para algún valor de x en el rango –2 <=x <= 1.

Método analítico:
f ´(x)=-4x3-6x2-16x-5=0
el valor de x = -0.34725
f ´´(x)=-12x2-12x-16
dado que la segunda derivada va a ser negativa hay un máximo en el punto
-0.34725
Método gráfico:
En el grafico se nota que el máximo de la función esta entre -0.6 y -0.2, pero más
exacto el valor es de -0.34725 como ya se demostró en el método analítico.

13.4 Considere la función que sigue: f(x) = 3 + 6x + 5x2 + 3x3 + 4x4

Localice el mínimo por medio de encontrar la raíz de la derivada de
dicha función. Utilice el método de bisección con valores iniciales
de xl = –2 y xu = 1
f ´(x) = 6+10x+9x2 +16x3 = 0
utilizando el método de la bisección visto en los capítulos anteriores nos sale un valor
aproximado para esta ecuación de x = -0.58668
f ´(x) = 10+18x+48x2
dado que la segunda derivada siempre dará positivo para cualquier valor de x, habrá
un mínimo en el punto -0.58668
f (-0.58668) = 3+6(-0.58668) +5(-0.58668)2+3(-0.58668)3+4(-0.58668)4
f (-0.58668) =1.0694

Gráfica de la función.

CAPÍTULO 14
14.1 Encuentre la derivada direccional de f (x, y) = x2 + 2y2, si x = 2 y y = 2, en la
dirección de h = 2i + 3j.
(d f(x,y))/dx=2x
(d f(x,y))/dx=2(2)=4
(d f(x,y))/dy=4y
(d f(x,y))/dy=4(2)=8
Entonces el gradiente es 4i+8j y nos pie en la dirección de h= 2i+3j.
Sabemos que la derivada direccional de la función en dirección “v” es igual al producto
escalar del gradiente con “v”
 (4,8). (2,3)
8+24
32
Entonces la tasa de cambio en función a “h” es de 32.
14.2 Repita el ejemplo 14.2 para la función siguiente, en el punto (0.8, 1.2).
f (x, y) = 2xy + 1.5y – 1.25x2 – 2y2 + 5
En el ejemplo 14.2 se halla el gradiente de la función para hallar la máxima pendiente:
Determinamos la elevación:
f (0.8,1.2) = 2(0.8) (1.2) +1.5(1.2)-1.25(0.8)2-2(1.2)2+5
f (0.8,1.2) = 5.04
evaluamos las derivadas parciales:
(d f(x,y))/dx=2y-2.5x=2(1.2)-2.5(0.8)=0.4
(d f(x,y))/dy=2x+1.5-4y=2(0.8)+1.5-4(1.2)=-1.7
Y el gradiente es 0.4i -1.7j, además la pendiente es la magnitud del gradiente

√ (〖 (0.4) 〗^2+〖 (-1.7) 〗^2) =1.7464

Esto quiere decir que en el primer paso se ha ganado 1.7464 de aumento de elevación
por unidad de distancia recorrida a lo largo de la trayectoria con la mayor pendiente.

14.3 Dada f (x, y) = 2.25xy + 1.75y – 1.5x2 – 2y2 construya y resuelva un sistema
de ecuaciones algebraicas lineales que maximice f(x). Observe que esto se logra
por medio de igualar a cero las derivadas parciales de f con respecto tanto de x
como de y.

𝑑 𝑓(𝑥, 𝑦)
= 2.25𝑦 − 3𝑥 = 0
𝑑𝑥
𝑑 𝑓(𝑥, 𝑦)
= 2.25𝑥 + 1.75 − 4𝑦 = 0
𝑑𝑦
Resolviendo esta ecuación el valor de x=0.56756 e y=0.75675
Usamos el hessiano
 𝑑 2 𝑓(𝑥, 𝑦)
= −3
(𝑑𝑥)2
𝑑 2 𝑓(𝑥, 𝑦)
= −4
(𝑑𝑦)2
𝑑 2 𝑓(𝑥, 𝑦)
= 2.25
𝑑𝑥𝑑𝑦
Con esto el hessiano sale 6.938, se concluye que existe un mínimo local en el punto
x=0.56756 e y=0.75675
14.4 Desarrolle una ecuación unidimensional en la dirección del gradiente de
presión en el punto (4, 2). La función de presión es f(x, y) = 6x2y – 9y2 – 8x2
Evaluamos las derivadas parciales de la ecuación en los puntos (4,2)
(d f(x,y))/dx=12xy-16x=12(4)(2)-16(4)=32
(d f(x,y))/dy=6x^2-18y=6(4)^2-18(2)=60
El gradiente es 32i+60j, luego hallamos la ecuación unidimensional con la siguiente
formula.

f (4+32h,2+60h) =6(4+32h)2(2+60h)-9(2+60h)2-8(4+32h)2
f (4+32h,2+60h) =368640h3+63856h2+4624h+28
entonces la función unidimensional en la dirección del gradiente será:
g(h)= 368640h3+63856h2+4624h+28
MATLAB OTROS

PROBLEMA 11
clear all;

clc

%F=x.^3+4*x.^3-3*x

syms x

f=input('INGRESE LA FUNCION: ');

integral=int(f,x);

disp('LA INTEGRAL ES: ');

disp(integral);

ezsurf(f);

pause(4);

hold on;

ezsurf(integral);

hold on;

PROBLEMA 12

clear all;

clc

syms x y
f=(x.^3+2*x.^2-y+3)

integral=int(int(f,y),(x));

disp('LA INTEGRAL ES: ');

disp(integral);

figure

ezsurf(f);

figure

%grafica de la funcion integrada

ezsurf(integral);

figure

pause(2);

%grafica de ambas funciones

ezsurf(f)

hold on;

ezsurf(integral);

PROBLEMA 13

clear all;

clc

syms x y
%f=x.^2+y.^2+1

f=input('INGRESE LA FUNCION: ');

integral=int(int(f,y),(x));

disp('LA INTEGRAL ES: ');

disp(integral);

figure

ezsurf(f);

figure

%grafica de la funcion integrada

ezsurf(integral);

figure

pause(2);

%grafica de ambas funciones

ezsurf(f)

hold on;

ezsurf(integral);

PROBLEMA 14

%======integral doble definida=====

clear all;
clc

syms x y

%f=x.^2+x-y.^6+y

f=input('INGRESE LA FUNCION: ');

integral=int(int(f,y),(x));

disp('LA INTEGRAL ES: ');

disp(integral);

figure

ezsurf(f);

figure

%grafica de la funcion integrada

ezsurf(integral);

figure

pause(2);

%grafica de ambas funciones

ezsurf(f)

hold on;

ezsurf(integral);

PROBLEMA 15

%======integral definida=====

clear all;
clc

syms x y

%f=(x+y.^3*sin(x-2))

%intervalo 0 1

f=input('INGRESE LA FUNCION: ');

a=input('INGRESE EL INTERVALO INFERIOR: ');

b=input('INGRESE EL INTERVALO SUPERIOR: ');

integral=int(f,x,a,b);

disp('LA INTEGRAL ES: ');

disp(integral);

figure

ezsurf(f);

figure

%grafica de la funcion integrada

ezsurf(integral);

figure

pause(2);

%grafica de ambas funciones

ezsurf(f)

hold on;

ezsurf(integral);

PROBLEMA 17

clear all;

clc

syms x y
%f=(2*x+sec(x)-x*y.^2)

%intervalo 0 1 1 1

f=input('INGRESE LA FUNCION: ');

a=input('DESDE|X|: ');

b=input('HASTA|X|: ');

a1=input('DESDE|Y|: ');

b1=input('HASTA|Y|: ');

F=int(int(f,x,a,b),y,a1,b1);

disp('LA INTEGRAL ES: ');

disp(F);

figure

ezsurf(f);

figure

%grafica de la funcion integrada

ezsurf(F);

figure

pause(2);

%grafica de ambas funciones

ezsurf(f)

hold on;

ezsurf(F);

PROBLEMA 18

clear all;

clc

syms x y
%f=(2*(x.^2*y.^2)+2*y)

%intervalo 1 2 0 1

f=input('INGRESE LA FUNCION: ');

a=input('DESDE|X|: ');

b=input('HASTA|X|: ');

a1=input('DESDE|Y|: ');

b1=input('HASTA|Y|: ');

F=int(int(f,x,a,b),y,a1,b1);

disp('LA INTEGRAL ES: ');

disp(F);

figure

ezsurf(f);

figure

%grafica de la funcion integrada

ezsurf(F);

figure

pause(2);

%grafica de ambas funciones

ezsurf(f)

hold on;

ezsurf(F);
 EJERCICIO 19

clear all

clc

syms x y z

%f=(sin(x)+cos(6*y)-14*z)

f=input('DIGITE LA FUNCION:')

integral=int(int(int(f,x),y),z);

disp('LA INTEGRAL ES: ');

disp(integral);

pause(7);

EJERCICIO 20

clear all

clc

syms x y z;

%f=(sin(2*x)^3-5*cos(2*y)-sin(2*z)-24)

f=input('DIGITE LA FUNCION:');

a=input('DESDE|X|:');

b=input('HASTA|X|:');

a1=input('DESDE|Y|:');

b1=input('HASTA|Y|:');

a2=input('DESDE|Z|:');

b2=input('HASTA|Z|:');

integral=int(int(int(f,x,a,b),y,a1,b1),z,a2,b2);

disp('LA INTEGRAL ES: ');

disp(integral);

hold on;
 CAPITULO 15
PROBLEMAS PROPÙESTOS

x1: cantidad de gas de la calidad regular a producir

x2: cantidad de gas de la calidad premium a producir

Max(Z)= 150x1+175x2

7x1 + 11x2 ≤ 77 Restricción respecto a la materia prima

10x1 + 8x2 ≤ 80 Restricción de tiempo de producción
x1 ≤ 9 Restricción de almacenamiento
x2 ≤ 6 Restricción de almacenamiento

regular premium total disponible

Producido 4.888889 3.888889

materia prima 7 11 77 77
tiempo 10 8 80 80
almacen. gas regular 4.888889 9
almacen. Gas premium 3.888889 6

ganacia por unidad 150 175

ganacia 733.333333 680.555556 1413.889
 Max(Z)= 150x1+175x2

(0,6) 1050
(1.57;6) 1286
(4.889;3.889) 1414
(0;8) 1400

Método Gráfico
12
10
8
(0;6)
6
(1.57;6)
4 (4.889;3.889)
2 y = -0.6364x + 7

0 (0;8)
0 2 4 6 8 10 12
x1: cantidad de gas de la calidad regular a producir
x2: cantidad de gas de la calidad premium a producir
x3: cantidad de gas crudo a producir

Max(Z)= 150x1+175x2+250x3

Sujeto a:
7x1 + 11x2 + 15x3 ≤ Restricción respecto a la materia
154 prima
10x1 + 8x2 + 12x3 ≤
80 Restricción de tiempo de producción
x1 ≤ 9 Restricción de almacenamiento
x2 ≤ 6 Restricción de almacenamiento
x3 ≤ 5 Restricción de almacenamiento

regular premium crudo total disponible

Producido 0 6 2.6666667

materia prima 7 11 15 106 154

tiempo 10 8 12 80 80
almacen. gas regular 0 9
almacen. Gas premium 6 6
almacen. Gas crudo 2.66666667 5
ganacia por unidad 150 175 250
ganacia 0 1050 666.666667 1716.66667
GRACIAS