Trabajo de métodos numéricos: Solución a los problemas

propuestos
Franz Gracian , Daniel Hernandez, Andres Hernandez y Mario Martinez

Sección 4.3
1b. Procedemos a utilizar el método del trapecio para aproximar la integral
Z 0,5
2
dx
0 x−4
como
b
b−a
Z
f (x)dx ' (f (a) + f (b))
a 2
entonces definimos a = 0 y b = 0,5 como
2
f (x) =
x−4
procedemos a buscar las imágenes, por lo tanto
2 1
f (a) = f (0) = =−
0−4 2
y
2 4
f (b) = f (0, 5) = =−
0, 5 − 4 7
reemplazando en la formula tenemos que
Z 0,5
2 0, 5 1 4 15
dx ' (− − ) = − und2
0 x − 4 2 2 7 56
6c. Procedemos a utilizar el método de simpson para aproximar la integral
Z 1,3
[sin2 (x) − 2xsinx + 1]dx
0,75

como
b
b−a
Z
a+b
f (x)dx ' (f (a) + 4f ( ) + f (b))
a 6 2
a+b
sea a = 0, 75 , b = 1, 3 y 2 = 1, 025 como

f (x) = sin2 (x) − 2xsinx + 1

procedemos a a calcular las imágenes

f (a) = f (0, 75) = sin2 (0, 75) − 2(0, 75)sin(0, 75) + 1 = 0, 9805369442

f (b) = f (1, 3) = sin2 (1, 3) − 2(1, 3)sin(1, 3) + 1 = 0, 9415276478

y
a+b
f( ) = f (1, 025) = sin2 (1, 025) − 2(1, 025)sin(1, 025) + 1 = 0, 9636482298
2
procedemos a reemplazar en la formula
Z 1,3
[sin2 (x) − 2xsin(x) + 1]dx
0,75

1
 1, 3 − 0, 75
' ((0, 9805369442) + 4(0, 9636482298) + (0, 9415276478)) = 0, 5295269385 und2
6
Ejercicio 5.2
1b. Procedemos a utilizar el método de Euler para aproximar la solución de la ecuación diferencial
dada.
y 0 = 1 + (t − y)2 , 2 6 t 6 3, y(2) = 1 h = 0, 5
dado que
dy
= f (t, y), a 6 t 6 b, y(a) = α
dx
ti = a + ih
y
wi+1 = wi + h[1 + (ti − wi )2 ]
sea a = 2 , b = 3 y α = 1 entonces
ti = 2 + 0, 5i
e iterando desde
i = 1, 2
obtenemos los siguientes resultados:
cuando i = 0
t0 = 2 w1 = 1
cuando i = 1

t1 = 2, 5 w2 = w1 + h[1 + (t1 − w1 )2 ] = 1 + 0, 5[1 + (2 − 1)2 ] = 2

cuando i = 2

t2 = 3 w3 = w2 + h[1 + (t2 − w2 )2 ] = 2 + 0, 5[1 + (2, 5 − 2)2 ] = 2, 625

por lo que podemos llevar estos datos obtenidos en una tabla

i ti wi+1
0 2 1
1 2, 5 2
2 3 2, 625

Ejercicio 5.4
1b. Procedemos a utilizar el método de Taylor de orden 4 para aproximar la solución de la
ecuación diferencial.

y 0 = 1 + (t − y)2 ; 2 6 t 6 3; y(2) = 1 ; h = 0, 5

de donde
f (t, y(t)) = 1 + (t − y)2
por lo que tenemos

f 0 (t, y(t)) = 2(t − y)(1 − y 0 ) = 2(t − y)(1 − 1 − (t − y)2 ) = −2(t − y)3

f 00 (t, y(t)) = −6(t − y)2 (1 − y 0 ) = −6(t − y)2 (1 − 1 − (t − y)2 ) = 6(t − y)4

f 000 (t, y(t)) = 24(t − y)3 (1 − 1 − (t − y)2 ) = −24(t − y)5
como sabemos,

h 0 h2 h3
T (4) (ti , wi ) = f (ti , wi ) + f (ti , wi ) + f 00 (ti , wi ) + f 000 (ti , wi )
2 6 24
procedemos a reemplazar las derivadas encontradas , por lo que tenemos

h h2 h3
T (4) (ti , wi ) = 1 + (ti − wi )2 + (−2)(ti − wi )3 + (6)(ti − wi )4 + (−24)(ti − wi )5
2 6 24

2
 = 1 + (ti − wi )2 − h(ti − wi )3 + h2 (ti − wi )4 − h3 (ti − wi )5
sabemos que la ecuación recursiva se denota por
wi+1 = wi + hT (n) (ti , wi )
en nuestro caso tenemos lo siguiente
wi+1 = wi + h[1 + (ti − wi )2 − h(ti − wi )3 + h2 (ti − wi )4 − h3 (ti − wi )5 ]
si i = 0
t0 = 2 y w1 = 1
si i = 1
t1 = 2, 5 y w2 = w1 + h[1 + (t1 − w1 )2 − h(t1 − w1 )3 + h2 (t1 − w1 )4 − h3 (t1 − w1 )5 ]
t1 − w1 = 2, 5 − 1 = 1, 5
993
w2 = 1 + 0, 5[1 + (1, 5)2 − 0, 5(1, 5)3 + 0, 52 (1, 5)4 − 0, 53 (1, 5)5 ] = = 1, 939453125
512
si i = 2
t2 = 3 y w3 = w2 + h[1 + (t2 − w2 )2 − h(t2 − w2 )3 + h2 (t2 − w2 )4 − h3 (t2 − w2 )5 ]
993 543
t2 − w2 = 3 − = = 1, 060546875
512 512
993 543 2 543 3 543 4 543 5
w3 = + 0, 5[1 + ( ) − 0, 5( ) + 0, 52 ( ) − 0, 53 ( ) ] = 2, 777898337
512 512 512 512 512
con estos datos obtenidos procedemos a representarlos en una tabla para compararlos con
los valores exactos provenientes de la solución analı́tica de la ecuación diferencial, la cual es
1
y(t) = t +
1−t
i ti wi+1 Exacto yi = y(ti )
0 2 1 1
1 2, 5 1, 939453125 1, 833333333
2 3 2, 777898337 2, 5
Ejercicio 5.9
[Link] a realizar el metodo de Runge Kutta para aproximar las soluciones del sistema
de ecuaciones propuesto. Nuestro sistema de ecuaciones esta dado por:

u01 = −4u1 − 2u2 + cos(t) + 4sin(t)

u02 = 3u1 + u2 − 3sin(t)
con las condiciones
u1 (0) = 0 u2 (0) = −1 h = 0, 1 0 6 t 6 2
Cabe denotar que conocemos la solución analı́tica de las mismas, las cuales son
U1 (t) = 2e−t − 2e−2t + sin(t) ; U2 (t) = −3e−t + 2e−2t
el método descrito está definido en nuestro caso de la siguiente forma sea
du1 du2
= f (t, u1 , u2 ) y = g(t, u1 , u2 )
dt dt
tenemos que
{i+1} {i} 1 {i+1} {i} 1
u1 = u1 + (K1 + 2K2 + 2K3 + K4 ) y u2 = u2 + (L1 + 2L2 + 2L3 + L4 )
6 6
donde
K1 = hf (ti , u1 (ti ), u2 (ti )) L1 = hg(ti , u1 (ti ), u2 (ti ))
h K1 L1 h K1 L1
K2 = hf (ti + , ui (ti ) + , u2 (ti ) + ) L2 = hg(ti + , ui (ti ) + , u2 (ti ) + )
2 2 2 2 2 2
h K2 L2 h K2 L2
K3 = hf (ti + , u1 (ti ) + , u2 (ti ) + ) L3 = hg(ti + , u1 (ti ) + , u2 (ti ) + )
2 2 2 2 2 2
K4 = hf (ti + h, u1 (ti ) + K3 , u2 (ti ) + L3 ) L4 = hg(ti + h, u1 (ti ) + K3 , u2 (ti ) + L3 )
A continuación se muestran los resultados de las iteraciones y de los valores exactos.

3
 t U1 U2 u1 u2 k1 k2 k3 K4 L1 L2 L3 L4
0 0 -1 0 -1 0,3 0,269866694 0,27339273 0,24572938 -0,1 -0,07499375 -0,07826343 -0,05575855
0,1 0,27204675 -1,07705075 0,27204137 -1,07704549 0,246026332 0,22164386 0,22459831 0,20213123 -0,05604216 -0,03682174 -0,03951808 -0,02226525
0,2 0,49549074 -1,11555217 0,49548169 -1,11554333 0,202390381 0,182541803 0,18505025 0,16666655 -0,02251063 -0,00789799 -0,01014464 0,00293472
0,3 0,67953338 -1,12483139 0,67952186 -1,12482019 0,166897025 0,150600435 0,15276403 0,13756731 0,00271848 0,01367568 0,01177905 0,02155615
0,4 0,83140051 -1,11230221 0,83138741 -1,11228951 0,137776376 0,124242411 0,12614217 0,11341343 0,02136177 0,02943216 0,02780558 0,03498282
0,5 0,95672798 -1,0838331 0,9567139 -1,0838195 0,113606812 0,102203866 0,10390418 0,09306982 0,03480456 0,0406073 0,039187 0,04432943
0,6 1,05987732 -1,04404648 1,05986269 -1,0440324 0,093251957 0,083477674 0,08502926 0,07563478 0,04416282 0,04819558 0,04693107 0,05049214
0,7 1,14419437 -0,99656198 1,14417945 -0,99654769 0,075809051 0,067266962 0,06870921 0,06039773 0,05033376 0,05299548 0,05184726 0,05418973
0,8 1,21222098 -0,94419386 1,21220597 -0,94417953 0,060566624 0,052946967 0,05430994 0,04680448 0,05403701 0,05564656 0,05458409 0,05599716
0,9 1,26586845 -0,8891112 1,26585346 -0,88909695 0,046969766 0,040033736 0,04133981 0,03442876 0,05584827 0,05665957 0,05565973 0,05637296

4
1 1,3065593 -0,83296776 1,3065444 -0,83295364 0,034591594 0,02815838 0,02942358 0,02294873 0,05622666 0,05644106 0,0554868 0,0556815
1,1 1,33534321 -0,77700693 1,33532844 -0,776993 0,02310978 0,017045885 0,01828103 0,0121277 0,05553702 0,05531337 0,0543926 0,05421108
1,2 1,3529916 -0,72214673 1,35297699 -0,72213299 0,012287211 0,006497636 0,00770927 0,00179813 0,05406807 0,0535309 0,0526356 0,05218868
1,3 1,36007462 -0,66904822 1,36006018 -0,6690347 0,001956023 -0,00362304 -0,00243185 -0,00815155 0,05204713 0,05129334 0,05041879 0,04979199
1,4 1,35702353 -0,61817077 1,3570093 -0,61815747 -0,00799562 -0,01340284 -0,01223183 -0,01778789 0,04965212 0,04875641 0,04790054 0,04715905
1,5 1,34418117 -0,56981634 1,34416716 -0,56980329 -0,01763449 -0,02288845 -0,02173959 -0,02714157 0,04702132 0,04604058 0,04520345 0,04439621
1,6 1,32184223 -0,52416515 1,32182847 -0,52415236 -0,02699143 -0,03209483 -0,03097184 -0,036215 0,04426122 0,04323814 0,04242148 0,04158446
1,7 1,29028532 -0,48130403 1,29027184 -0,48129154 -0,03606895 -0,04101216 -0,03992012 -0,04498869 0,04145296 0,04041892 0,03962574 0,03878465
1,8 1,24979796 -0,44124922 1,24978481 -0,44123705 -0,04484767 -0,04961167 -0,04855661 -0,05342662 0,03865745 0,0376349 0,03686917 0,03604164
1,9 1,20069578 -0,40396431 1,200683 -0,40395251 -0,05329162 -0,05785053 -0,05683919 -0,06148081 0,03591962 0,03492397 0,03419035 0,0333877
2 1,14333672 -0,36937457 1,14332436 -0,36936318 -0,06135282 -0,06567604 -0,06471551 -0,06909507 0,03327176 0,03231294 0,03161652 0,03084518