Universidad Andrés Bello Curso: ICI 2204 - Métodos matemáticos

Facultad de Ingenierı́a Semestre: 1-2018

Ingenierı́a Civil Industrial Profesores: Diego Beneventti - Samantha Reid
Juan Carlos Velazquez

Control recuperativo (b)

Plazo: Domingo 24 de Junio hasta las 11:59
Instrucciones:

 Esta evaluación se podrá realizar con un máximo de 2 integrantes.

 Se debera de anexar cada uno de los calculos que se realicen para la resolución del problmea.
 Utilice matlab en caso de ser necesario para encontrar la solución.

Nombre alumno:

1. Problema 1

Mediante el método de la secante calcule la raı́z de la siguiente función tomando en cuenta como
puntos iniciales x0 = −2 y x1 = −1 y un error de 0.0001

f (x) = e−x (x2 + 5x + 2) + 1

iter xi−1 xi f (xi−1 ) f (xi ) xi+1 Error

1 -2 -1 -28.5562244 -4.436563657 -0.816060279
2 -1 -0.816060279 -4.436563657 -2.198648037 -0.635348057 -0.284430275
3 -0.816060279 -0.635348057 -2.198648037 -0.45931395 -0.587626549 -0.081210605
4 -0.635348057 -0.587626549 -0.45931395 -0.066919168 -0.579488104 -0.014044198
5 -0.587626549 -0.579488104 -0.066919168 -0.002586656 -0.579160876 -0.000565002
6 -0.579488104 -0.579160876 -0.002586656 -1.54772E-05 -0.579158907 -3.40104E-06

Cuadro 1: Tabla de resultados

2. Problema 2

Dada la siguiente ecuación apliqué el método de falsa posición para poder encontrar su solución
utilizando un intervalo de a = 3 y b = 4, tomando en cuenta un error del 0.0001
f (x) = e−x (3,2sen (x) − 0,5cos (x))
 Iter a b x error
1 3 4 3.551212931
2 3 3.551212931 3.368301562 5.43037
3 3 3.368301562 3.314887572 1.61134
4 3 3.314887572 3.301127814 0.41682
5 3 3.301127814 3.297706885 0.10374
6 3 3.297706885 3.296864057 0.02556
7 3 3.296864057 3.296656873 0.00628
8 3 3.296656873 3.296605971 0.00154
9 3 3.296605971 3.296593467 0.00038

Cuadro 2: Tabla de resultados

3. Problema 3

Considere el siguiente sistema de ecuaciones, encuentre los valores de x,y y z para los cuales las
funciones se acercan a cero con un error de 0.01.
Los puntos iniciales son (0,1, 0,1, −0,1)

1
3x − cos(yz) − =0
2
x2 − 81(y + 0,1)2 + sen(z) + 1,06 = 0
10π − 3
e−xy + 20z + =0
3

(a) Despeje
1 1
x = cos(yz) +
3 6
1 p
y= x + sen(z) + 1,06 − 0,1
9
1 10π − 3
z = − e−xy −
20 60
(b) Solución

iter xi+1 yi+1 zi+1 f (xi+1 ) f (yi+1 ) f (zi+1 )

1 0.499983333 0.00944115 -0.523102492 -3.78043E-05 -0.159753394 0.005240872
2 0.499995935 2.55023E-05 -0.523364535 -1.21952E-05 -0.000214383 0.004696543
3 0.5 1.22713E-05 -0.523599362 -6.84293E-11 -0.000199315 6.61535E-06
4 0.5 -3.13741E-08 -0.523599693 -2.06416E-11 -2.8643E-07 6.1513E-06
5 0.5 -4.90549E-08 -0.523600001 0 -2.66352E-07 8.84043E-09

Cuadro 3: Tabla de resultados

4. Problema 4

Se recolectaron los siguientes datos para una sección transversal de un rı́o (y = distancia de una
ribera, H = profundidad y U = velocidad):

y, m 0 1 3 5 7 8 9 10
H, m 0 1 1.5 3 3.5 3.2 2 0
U , m/s 0 0.1 0.12 0.2 0.25 0.3 0.15 0

Use integración numérica para calcular:

a) El área de la sección transversal

b) El caudal

Observe que el área de la sección transversal (Ac ) y el caudal (Q) se puede calcular como

Z y
Ac = H(y)dy
0

Z y
Q= H(y)U (y)dy
0

Respuesta

a) el área de la sección transversal

Z y
Ac = H(y)dy
0
Z 10
= H(y)dy
0

Aplicando Simpson 3/8:

3h 2a + b a + 2b
I= [f (a) + 3 ∗ f ( ) + 3 ∗ f( ) + f (b)]
8 3 3
3 ∗ 3,333
= [f (0) + 3 ∗ f (3) + 3 ∗ f (7) + f (10)]
8
3 ∗ 3,333
= (0 + 3 ∗ 1,5 + 3 ∗ 3,5 + 0) = 18,75
8

b) El caudal tiene la siguiente ecuación:

Z y
Q= H(y)U (y)dy
0
Z 10
= H(y)U (y)dy
0
 Aplicando simpson 3/8 tenemos:

3h 2a + b a + 2b
I= [f (a) + 3 ∗ f ( ) + 3 ∗ f( ) + f (b)]
8 3 3
3 ∗ 3,333
= [f (0) + 3 ∗ f (3) + 3 ∗ f (7) + f (10)]
8
3 ∗ 3,333
= ((0 ∗ 0,1) + 3(1,5 ∗ 0,12) + 3(3,5 ∗ 0,3) + 0) = 3,95
8

5. Problema 5

La función f (x) = 2e−1,5x se puede utilizar para generar la tabla siguiente de datos espaciados en
forma desigual:

x 0 0.05 0.15 0.25 0.35 0.475 0.6

f (x) 2 1.8555 1.5970 1.3746 1.1831 0.9808 0.8131

Evalúe la integral de a = 0 a b = 0,6, con el uso de :

a) la regla del trapecio

b) una combinación de las reglas de trapecio y de Simpson; emple la regla de Simpson siempre
que sea posible a fin de obtener la exactitud más alta.

c) Calcular el error relativo porcentual para cada una de ellas.

Respuesta

a) Al ocupar Trapecio tenemos:

h
I= [f (x0 ) + f (x1 )]
2
0,6
= [2 + 0,8131]
2
= 0,84393

b) Para este ejercicio, tenemos que analizar los intervalos y ver cuales puntos presentan la misma
distancia entre ellos. Analizando esto, ocupar los siguientes métodos:

 Trapecio: 0 - 0.05, con h = 0.05

 Trapecio: 0.05 - 0.15, con h = 0.1

 Simpson 13 : 0.15 - 0.35, con h = 0.1

 Simpson 13 : 0.35 - 0.6, con h = 0.125

Y se calcula:
  1◦ Trapecio:
0,05
I= [2 + 1,8555] = 0,096388
2
 2◦ Trapecio:
0,1
I= [1,8555 + 1,5970] = 0,172625
2
 1◦ Simpson 13 :
0,1
I= [1,5970 + 4 ∗ 1,3746 + 1,1831] = 0,27595
3
 1◦ Simpson 13 :

0,125
I= [1,1831 + 4 ∗ 0,9808 + 0,8131] = 0,246642
3

Sumando todos los elementos, da como resultado I = 0,791605

c) Calculamos el valor real de la integral, el cual es 0.79124

Analizamos para ambos casos:

 Para (a) tenemos un valor t = 6,65917 %

 Para (b) tenemos un valor t = 0,04613 %

6. Problema 6

Considere una función f (x) tal que f (2) = 1,5713, f (3) = 1,5719, f (5) = 1,5738 y f (6) = 1,5751.
Calcule f (4) usando un polinomio de interpolación de segundo grado y uno de tercer grado. Redondee
los resultados finales a cuatro lugares decimales. ¿Existe alguna ventaja aquı́ al usar un polinomio
de tercer grado?

Respuesta

f2 (x) = b0 + b1 (x − x0 ) + b2 (x − x0 ) (x − x1 )
Entonces lo que se tiene que calcular son los b sub indices.

b0 = f (x0 )

f (x1 )−f (x0 )

b1 = x1 −x0
f (x2 )−f (x1 ) f (x )−f (x )
x2 −x1
− x1 −x 0
b2 = x2 −x0
1 0

f2 (4) = b0 + b1 (x − x0 ) + b2 (x − x0 ) (x − x1 )

b0 = f (5) = 1,5738

f (6)−f (5)
b1 = 6−5
f (3)−f (5) f (5)−f (3)
6−5
− 5−3
b2 = 6−5 = −0,00285
 1,5751−1,5738
b1 = 6−5 = 0,0013

f2 (4) = 1,5738 + 0,0013(4 − 5) + −0,00285(4 − 5)(4 − 6) = 1,5668

7. Problema 7

La función siguiente muestra regiones tanto inclinadas en una región de x relativamente corta:

1 1
f (x) = + −6
(x − 0,3) + 0,01 (x − 0,9)2 + 0,04
2

Determine el valor de la integral definida de la función entre x 0 y 1, resuelva el problema mediante

las siguientes metodologı́as.

 Euler

 Taylor

 Kutta 2 y 4

Respuesta

Debido a que no se presenta el tamaño de paso, y considerando como supuesto que f (0) = 0, tenemos
múltiples resultados dependiendo de estos elementos. A continuación, daremos algunos resultados:

(a) Para Euler, con tamaño de paso 0.1, tenemos:

Iteración xi yi h yi+1
1 0 0 0.1 0.51764706
2 0.1 0.51764706 0.1 2.06470588
3 0.2 2.06470588 0.1 6.65338513
4 0.3 6.65338513 0.1 16.3033851
5 0.4 16.3033851 0.1 21.0482127
6 0.5 21.0482127 0.1 22.9482127
7 0.6 22.9482127 0.1 24.1174435
8 0.7 24.1174435 0.1 25.3556788
9 0.8 25.3556788 0.1 27.1402942
10 0.9 27.1402942 0.1 29.3105644
(b) Para Taylor, con tamaño de paso 0.1, tenemos:

Iteración x y h yi+1
1 0 0 0.1 0.80519031
2 0.1 0.80519031 0.1 3.1349481
3 0.2 3.1349481 0.1 10.1987074
4 0.3 10.1987074 0.1 19.8112074
5 0.4 19.8112074 0.1 21.996582
6 0.5 21.996582 0.1 22.996582
7 0.6 22.996582 0.1 23.688298
8 0.7 23.688298 0.1 24.475625
9 0.8 24.475625 0.1 25.7862758
10 0.9 25.7862758 0.1 27.9127185

(c) Para RK2, con tamaño de paso 0.1, tenemos:

Iteración x y h yi+1’ xi+1 yi+1

1 0 0 0.1 0.51764706 0.1 1.032352941
2 0.1 1.03235294 0.1 2.57941176 0.2 4.100221976
3 0.2 4.10022198 0.1 8.68890122 0.3 11.2195616
4 0.3 11.2195616 0.1 20.8695616 0.4 18.41697539
5 0.4 18.4169754 0.1 23.161803 0.5 21.73938918
6 0.5 21.7393892 0.1 23.6393892 0.6 23.27400457
7 0.6 23.2740046 0.1 24.4432353 0.7 24.4777376
8 0.7 24.4777376 0.1 25.7159729 0.8 25.98916294
9 0.8 25.9891629 0.1 27.7737783 0.9 27.96660577
10 0.9 27.9666058 0.1 30.136876 1 29.8517409

(d) Para RK4, con tamaño de paso 0.1, tenemos:

Iteración x y h x+1 y+1 k1 k2 k3 k4 x y

1 0 0 0.1 0.1 0.95108957 5.176470588 9.10457886 9.10457886 15.4705882 0.1 0.91045789
2 0.1 0.95108957 0.1 0.2 4.10227806 37.46858987 26.4289818 26.4289818 45.8867925 0.2 3.59398775
3 0.2 4.10227806 0.1 0.3 11.5230412 44.09713845 76.1621622 76.1621622 96.5 0.3 11.7184943
4 0.3 11.5230412 0.1 0.4 19.0086406 94.00885826 76.919708 76.919708 47.4482759 0.4 19.215012
5 0.4 19.0086406 0.1 0.5 21.9848876 44.00304913 28.8929421 28.8929421 19 0.5 21.8979348
6 0.5 21.9848876 0.1 0.6 23.3429268 14.00224915 13.9469496 13.9469496 11.6923077 0.6 23.3795825
7 0.6 23.3429268 0.1 0.7 24.3695503 4.00198521 11.3032674 11.3032674 12.3823529 0.7 24.4732536
8 0.7 24.3695503 0.1 0.8 25.6454478 -0.11583171 14.7058824 14.7058824 17.8461538 0.8 25.8401385
9 0.8 25.6454478 0.1 0.9 27.3532501 -2.15221317 20.7294118 20.7294118 21.7027027 0.9 27.718389
10 0.9 27.3532501 0.1 1 28.8877557 -3.29586835 19.8415505 19.8415505 16 1 29.3374051