Problemas de Valor Inicial-Unidad 5 - Calculo Numerico

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 22

UNIVERSIDAD TÉCNICA DE MANABÍ

Facultad de Ciencias Matemáticas, Físicas y Químicas

Carrera:

Ingeniera industrial

TEMA:

PROBLEMAS DE VALOR INICIAL

MATERIA:

CÁLCULO NUMÉRICO

DOCENTE:

ING. GABRIEL MOREIRA

Estudiante:

Leonardo Juliano Pachay Mero.

NIVEL:

PARALELO:

“C”
1. Resuelva el siguiente problema de valor inicial en el intervalo de t = 0 a 2, donde

𝒚(𝟎) = 𝟏. Muestre todos sus resultados en la misma gráfica.

𝑑𝑦
= 𝑦𝑡 3 − 1.1𝑦
𝑑𝑡

a. Analíticamente.

𝑑𝑦
= 𝑦𝑡 3 − 1.1𝑦
𝑑𝑡

𝑑𝑦
= 𝑦(𝑡 3 − 1.1)
𝑑𝑡

𝑑𝑦
= (𝑡 3 − 1.1)𝑑𝑡
𝑦

𝑑𝑦
∫ = ∫(𝑡 3 − 1.1)𝑑𝑡
𝑦

𝑡4
𝑙𝑛(𝑦) = − 1.1𝑡 + 𝐶
4

𝑡4
𝑦= 𝐶𝑒 4 −1.1𝑡

04
1= 𝐶𝑒 4 −1.1(0)

𝐶=1

𝑡4
𝑦= 𝑒 4 −1.1𝑡

b. Método de Euler con 𝒉 = 𝟎. 𝟓 y 𝒉 = 𝟎. 𝟐𝟓

𝑡1 = 𝑡0 + ℎ

𝑦1 = 𝑦0 + ℎ ∗ 𝑓(𝑡0 , 𝑦0 )
Euler con ℎ = 0.5

Iteración 𝑓(𝑥0 , 𝑦0 ) t y
0 1
1 -1.1 0.5 0.45
2 -0.439 1 0.2306
3 -0.023 1.5 0.2191
4 0.4984 2 0.4683

Euler con ℎ = 0.25

Iteración 𝑓(𝑥0 , 𝑦0 ) t y
0 1
1 -1.1 0.25 0.725
2 -0.786 0.5 0.5284
3 -0.515 0.75 0.3996
4 -0.271 1 0.3318
5 -0.033 1.25 0.3235
6 0.276 1.5 0.3925
7 0.893 1.75 0.6157
8 2.622 2 1.2713

c. Método de RK2 con 𝒉 = 𝟎. 𝟓


𝑦𝑖+1 = 𝑦1 + (𝑘1 + 𝑘2 )
2

𝑘1 = 𝑓(𝑡𝑖 , 𝑦𝑖 )

𝑘2 = 𝑓(𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘1 ℎ)

𝑓(𝑡𝑖 , 𝑦𝑖 ) = 𝑦𝑡 3 − 1.1𝑦

Iteración 𝑘1 𝑘2 t y
0 1
1 -1.1 -0.4388 0.5 0.6153
2 -0.5999 -0.0315 1 0.4574
3 -0.0457 0.9885 1.5 0.6931
4 1.5768 10.222 2 3.4071
d. Método RK4 con 𝒉 = 𝟎. 𝟐𝟓

1
𝑦𝑖+1 = 𝑦1 + ℎ(𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 )
6

𝑘1 = 𝑓(𝑡𝑖 , 𝑦𝑖 )

1 1
𝑘2 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘1 ℎ)
2 2

1 1
𝑘3 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘2 ℎ)
2 2

𝑘4 = 𝑓(𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘3 ℎ)

𝑓(𝑡𝑖 , 𝑦𝑖 ) = 𝑦𝑡 3 − 1.1𝑦

Iteración 𝑘1 𝑘2 𝑘3 𝑘4 t y
0 1
1 -1.1 -0.9471 -0.9681 -0.8219 0.25 0.7603
2 -0.8244 -0.6883 -0.7061 -0.5692 0.5 0.5860
3 -0.5714 -0.4404 -0.4544 -0.3203 0.75 0.4743
4 -0.3216 -0.1867 -0.1939 -0.0426 1 0.4274
5 -0.0427 0.1367 0.1439 0.3953 1.25 0.4655
6 0.3971 0.7725 0.8429 1.5384 1.5 0.6808
7 1.5488 2.7902 3.2854 6.3982 1.75 1.5182
8 6.4666 12.7768 17.1086 39.9879 2 5.9442

e. Compare sus resultados.

PROBLEMA 1 EDO
7

4 Método Analítico
Método Euler con h=0,5
y

3 Método Euler con h=0,25


Método RK2 con h=0,5
2
Método RK4 con h=0,25
1

0
0 0,5 1 1,5 2 2,5
t
Al analizar la gráfica podemos llegar a la conclusión de que el método que más acerca

sus resultados a los datos reales de la ecuación diferencial es el Método RK4 con ℎ =

0.25, dando una exactitud casi igual al de la gráfica de la función, en cambio los otros

métodos, aunque al principio no había mucha diferencia poco a poco se fueron

distanciando más de la solución real.

2. Resuelva de t = 0 a 3, con 𝒉 = 𝟎. 𝟏, con los métodos de Euler y RK2.

𝑑𝑦
= 𝑦𝑠𝑒𝑛3 (𝑡) 𝑦(0) = 1
𝑑𝑡

a. Analíticamente.

𝑑𝑦
= 𝑦𝑠𝑒𝑛3 (𝑡)
𝑑𝑡

𝑑𝑦
= 𝑠𝑒𝑛3 (𝑡)𝑑𝑡
𝑦

𝑑𝑦
∫ = ∫ 𝑠𝑒𝑛3 (𝑡)𝑑𝑡
𝑦

𝑙𝑛(𝑦) = ∫ 𝑠𝑒𝑛2 (𝑡)𝑠𝑒𝑛(𝑡)𝑑𝑡 = ∫(1 − 𝑐𝑜𝑠 2 (𝑡))𝑠𝑒𝑛(𝑡)𝑑𝑡

𝑙𝑛(𝑦) = ∫(𝑠𝑒𝑛(𝑡) − 𝑐𝑜𝑠 2 (𝑡)𝑠𝑒𝑛(𝑡))𝑑𝑡 = ∫ 𝑠𝑒𝑛(𝑡)𝑑𝑡 − ∫ 𝑐𝑜𝑠 2 (𝑡)𝑠𝑒𝑛(𝑡)𝑑𝑡

𝑐𝑜𝑠 3 (𝑡)
𝑙𝑛(𝑦) = −𝑐𝑜𝑠(𝑡) + ∫ 𝑐𝑜𝑠 2 (𝑡)(−𝑠𝑒𝑛(𝑡))𝑑𝑡 = −𝑐𝑜𝑠(𝑡) + +𝐶
3

𝑐𝑜𝑠 3 (𝑡)
𝑙𝑛(𝑦) = −𝑐𝑜𝑠(𝑡) + +𝐶
3

𝑐𝑜𝑠3 (𝑡)
𝑦 = 𝐶𝑒 −𝑐𝑜𝑠(𝑡)+ 3
𝑐𝑜𝑠3 (0)
1 = 𝐶𝑒 −𝑐𝑜𝑠(0)+ 3

2
𝐶 = 𝑒3

2 𝑐𝑜𝑠3 (𝑡)
−𝑐𝑜𝑠(𝑡)+
𝑦= 𝑒3 ∗𝑒 3

𝑐𝑜𝑠3 (𝑡) 2
𝑦 = 𝑒 −𝑐𝑜𝑠(𝑡)+ 3
+
3

b. Con los métodos propuestos.

Método de Euler:

𝑡1 = 𝑡0 + ℎ

𝑦1 = 𝑦0 + ℎ ∗ 𝑓(𝑡0 , 𝑦0 )

𝑓(𝑡0 , 𝑦0 ) = 𝑦𝑠𝑒𝑛3 (𝑡)

Iteración 𝑓(𝑥0 , 𝑦0 ) t y
0 1
1 0 0.1 1
2 5.31657e-09 0.2 1.000000001
3 4.25324e-08 0.3 1.000000005
4 1.43546e-07 0.4 1.000000019
5 3.40253e-07 0.5 1.000000053
6 6.64547e-07 0.6 1.00000012
7 1.14832e-06 0.7 1.000000234
8 1.82345e-06 0.8 1.000000417
9 2.72182e-06 0.9 1.000000689
10 3.87531e-06 1 1.000001077
11 5.31577e-06 1.1 1.000001608
12 7.07507e-06 1.2 1.000002316
13 9.18505e-06 1.3 1.000003234
14 1.16776e-05 1.4 1.000004402
15 1.45844e-05 1.5 1.00000586
16 1.79374e-05 1.6 1.000007654
17 2.17684e-05 1.7 1.000009831
18 2.61091e-05 1.8 1.000012442
19 3.09914e-05 1.9 1.000015541
20 3.64469e-05 2 1.000019186
21 4.25075e-05 2.1 1.000023436
22 4.92049e-05 2.2 1.000028357
23 5.65708e-05 2.3 1.000034014
24 6.46369e-05 2.4 1.000040478
25 7.34349e-05 2.5 1.000047821
26 8.29964e-05 2.6 1.000056121
27 9.33532e-05 2.7 1.000065456
28 0.000104537 2.8 1.00007591
29 0.000116579 2.9 1.000087568
30 0.000129511 3 1.000100519

Método de RK2:


𝑦𝑖+1 = 𝑦1 + (𝑘1 + 𝑘2 )
2

𝑘1 = 𝑓(𝑡𝑖 , 𝑦𝑖 )

𝑘2 = 𝑓(𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘1 ℎ)

𝑓(𝑡0 , 𝑦0 ) = 𝑦𝑠𝑒𝑛3 (𝑡)

Iteración 𝑘1 𝑘2 t y
0 1
1 0 0.000995011 0.1 1.000049751
2 5.31683e-09 0.00784177 0.2 1.000441839
3 4.25511e-08 0.025819831 0.3 1.001732833
4 1.43794e-07 0.059156317 0.4 1.004690656
5 3.41849e-07 0.1107123 0.5 1.010226288
6 6.71343e-07 0.181860895 0.6 1.019319366
7 1.1705e-06 0.272526199 0.7 1.032945735
8 1.88352e-06 0.381313303 0.8 1.052011494
9 2.86339e-06 0.505649686 0.9 1.077294122
10 4.17484e-06 0.641877119 1 1.109388186
11 5.89725e-06 0.785271907 1.1 1.148652077
12 8.12678e-06 0.930017638 1.2 1.195153365
13 1.09775e-05 1.069197348 1.3 1.248613781
14 1.45807e-05 1.194901361 1.4 1.308359578
15 1.90816e-05 1.298553707 1.5 1.373288218
16 2.46331e-05 1.371534729 1.6 1.441866186
17 3.13868e-05 1.40611425 1.7 1.512173467
18 3.94811e-05 1.396612047 1.8 1.582006044
19 4.90279e-05 1.340590462 1.9 1.649038018
20 6.01014e-05 1.239795734 2 1.71103081
21 7.27303e-05 1.100545769 2.1 1.766061735
22 8.68969e-05 0.933345193 2.2 1.81273334
23 0.000102545 0.751688515 2.3 1.850322893
24 0.000119595 0.570236971 2.4 1.878840721
25 0.000137967 0.402740195 2.5 1.898984629
26 0.000157601 0.260144429 2.6 1.911999731
27 0.000178481 0.149256791 2.7 1.919471494
28 0.000200642 0.072156347 2.8 1.923089344
29 0.000224175 0.026336418 2.9 1.924417373
30 0.000249212 0.005408423 3 1.924700255

Comparación de los datos obtenidos:

PROBLEMA 2 EDO
4

3,5

2,5

2 Método Analítico
y

1,5 Método de Euler


Método de RK2
1

0,5

0
0 0,5 1 1,5 2 2,5 3 3,5
t

Como se puede observar en la representación gráfica, el método que más se acerca al

valor real es el Método de RK2, no obstante, aún hay una distancia significativa entre la

gráfica real y la obtenida.


3. Use los métodos de Euler y RK4 para resolver:

𝑑𝑦
= −2𝑦 + 5𝑒 −𝑡
𝑑𝑡

𝑑𝑧 𝑦𝑧 2
=−
𝑑𝑦 2

En el rango de x = 0 a 0.4, con un tamaño de paso de 0.1, con 𝒚(𝟎) = 𝟐 y 𝒛(𝟎) = 𝟒

𝒅𝒚
Euler para = −𝟐𝒚 + 𝟓𝒆−𝒕
𝒅𝒕

𝑡1 = 𝑡0 + ℎ

𝑦1 = 𝑦0 + ℎ ∗ 𝑓(𝑡0 , 𝑦0 )

𝑓(𝑡0 , 𝑦0 ) = −2𝑦 + 5𝑒 −𝑡

Iteración 𝑓(𝑥0 , 𝑦0 ) t y
0 2
1 1 0.1 2.1
2 0.32418709 0.2 2.132418709
3 -0.171183653 0.3 2.115300344
4 -0.526509584 0.4 2.062649385

𝒅𝒚
RK4 para = −𝟐𝒚 + 𝟓𝒆−𝒕
𝒅𝒕

1
𝑦𝑖+1 = 𝑦1 + ℎ(𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 )
6

𝑘1 = 𝑓(𝑡𝑖 , 𝑦𝑖 )

1 1
𝑘2 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘1 ℎ)
2 2

1 1
𝑘3 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘2 ℎ)
2 2

𝑘4 = 𝑓(𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘3 ℎ)
𝑓(𝑡0 , 𝑦0 ) = −2𝑦 + 5𝑒 −𝑡

Iteración 𝑘1 𝑘2 𝑘3 𝑘4 t y
0 2
1 1 0.656147123 0.69053241 0.386080608 0.1 2.067990661
2 0.388205768 0.128737983 0.154684761 -0.073264509 0.2 2.082687107
3 -0.071720449 -0.264198254 -0.244950473 -0.412293016 0.3 2.057648592
4 -0.41120608 -0.550736127 -0.536783122 -0.656340329 0.4 2.00360551

Comparación de los datos obtenidos:

PROBLEMA 3 EDO 1
2,14

2,12

2,1

2,08

2,06
y

Método de Euler
2,04 Método de RK4

2,02

1,98
0 0,05 0,1 0,15 0,2 0,25 0,3 0,35 0,4 0,45
t

Las gráficas obtenidas son muy similares en cuanto a sus datos, no obstante, presentan

una diferencia notable en cuanto al desplazamiento de los datos en el eje de las “y”.

𝒅𝒛 𝒚𝒛𝟐
Euler para 𝒅𝒚 = − 𝟐

𝑡1 = 𝑡0 + ℎ

𝑦1 = 𝑦0 + ℎ ∗ 𝑓(𝑡0 , 𝑦0 )

𝑦𝑧 2
𝑓(𝑡0 , 𝑦0 ) = −
2
Iteración 𝑓(𝑥0 , 𝑦0 ) y z
0 2
1 0 0.1 2
2 -0.2 0.2 1.98
3 -0.39204 0.3 1.940796
4 -0.565003367 0.4 1.884295663

𝒅𝒛 𝒚𝒛𝟐
RK4 para𝒅𝒚 = − 𝟐

1
𝑦𝑖+1 = 𝑦1 + ℎ(𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 )
6

𝑘1 = 𝑓(𝑡𝑖 , 𝑦𝑖 )

1 1
𝑘2 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘1 ℎ)
2 2

1 1
𝑘3 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘2 ℎ)
2 2

𝑘4 = 𝑓(𝑡𝑖 + ℎ, 𝑦𝑖 + 𝑘3 ℎ)

𝑦𝑧 2
𝑓(𝑡0 , 𝑦0 ) = −
2

Iteración 𝑘1 𝑘2 𝑘3 𝑘4 y z
0 2
1 0 -0.1 -0.099500625 -0.198014938 0.1 1.99004973
2 -0.198014896 -0.29407425 -0.292649392 -0.3844677 0.2 1.960784232
3 -0.384467481 -0.471207321 -0.469104538 -0.549436926 0.3 1.91387543
4 -0.549437874 -0.622740703 -0.620323176 -0.685864583 0.4 1.851851593

Comparación de los datos obtenidos:


PROBLEMA 3 EDO 2
2,02
2
1,98
1,96
1,94
z

1,92 Método de Euler

1,9 Método de RK4

1,88
1,86
1,84
0 0,1 0,2 0,3 0,4 0,5
y

En este caso, ambos métodos acercan sus resultados, distanciándose por pocos decimales

y estando muy cerca sus valores.

4. La siguiente es una ecuación diferencial de segundo orden con valor inicial:

𝑑2𝑥 𝑑𝑥
2
+ (5𝑥) + (𝑥 + 7)𝑠𝑒𝑛(𝜔𝑡) = 0
𝑑𝑡 𝑑𝑡

Donde

𝑑𝑥
(0) = 1.5 𝑦 𝑥(0) = 6
𝑑𝑡

Observe que 𝜔 = 1. Descomponga la ecuación en dos ecuaciones diferenciales de

primer orden. Después de la descomposición, resuelva el sistema de t = 0 a 15, y

grafique sus resultados con el método que desee.

Descomposición:

𝑑2𝑥 𝑑𝑥
2
+ (5𝑥) + (𝑥 + 7)𝑠𝑒𝑛(𝑡) = 0
𝑑𝑡 𝑑𝑡

𝑑2𝑥 𝑑𝑥
= −(5𝑥) − (𝑥 + 7)𝑠𝑒𝑛(𝑡)
𝑑𝑡 2 𝑑𝑡
𝑑𝑥
𝑦1 = 𝑥 𝑦2 =
𝑑𝑡

𝑑𝑦1 𝑑𝑥 𝑑𝑦2 𝑑 2 𝑥
= = 2
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡

dy1 dy2
= y2 = −5y1 y2 − (y1 + 7)sen(t)
dt dt

Se utilizará el método de Euler para resolver el sistema de ecuaciones diferenciales, se

aplican los valores iniciales:

dy1
= y2 y1 (0) = 6
{ dt →{
dy2 y2 (0) = 1.5
= −5y1 y2 − (y1 + 7)sen(t)
dt

Para este método se utilizará un paso de 0.05, debido a que esta ecuación se trata de una

ecuación no lineal de segundo orden con términos trigonométricos, las cuales no tienen

soluciones analíticas generales sencillas, por ende, si se utilizaba un paso de mayor

distancia esto provocaba que los cálculos tiendan a llevar a un resultado impreciso o

inestable, por eso se optó por este tipo de paso, además de que hay que tener en cuenta

también el intervalo establecido de 0 a 15.

Iteración t 𝑥 = 𝑦1 𝑦2
0 6 1.5
1 0.05 6.075 -0.75
2 0.1 6.0375 0.356388618
3 0.15 6.055319431 -0.246614361
4 0.2 6.042988713 0.029169694
5 0.25 6.044447198 -0.144460431
6 0.3 6.037224176 -0.087526963
7 0.35 6.032847828 -0.148060148
8 0.4 6.02544482 -0.14820081
9 0.45 6.01803478 -0.178574216
10 0.5 6.009106069 -0.193027578
11 0.55 5.99945469 -0.214891665
12 0.6 5.988710107 -0.232315911
13 0.65 5.977094312 -0.25119662
14 0.7 5.964534481 -0.268518201
15 0.75 5.95110857 -0.285720805
16 0.8 5.93682253 -0.302030802
17 0.85 5.92172099 -0.317770407
18 0.9 5.90583247 -0.332725274
19 0.95 5.889196206 -0.346944636
20 1 5.871848974 -0.36035198
21 1.05 5.853831375 -0.372933251
22 1.1 5.835184713 -0.384646753
23 1.15 5.815952375 -0.395466094
24 1.2 5.79617907 -0.405360061
25 1.25 5.775911067 -0.414302138
26 1.3 5.755195961 -0.422266217
27 1.35 5.73408265 -0.429228684
28 1.4 5.712621215 -0.435167591
29 1.45 5.690862836 -0.440063145
30 1.5 5.668859679 -0.443897616
31 1.55 5.646664798 -0.446655493
32 1.6 5.624332023 -0.448323538
33 1.65 5.601915846 -0.448890883
34 1.7 5.579471302 -0.448349105
35 1.75 5.557053847 -0.446692315
36 1.8 5.534719231 -0.443917229
37 1.85 5.51252337 -0.440023254
38 1.9 5.490522207 -0.435012558
39 1.95 5.468771579 -0.428890143
40 2 5.447327072 -0.421663911
41 2.05 5.426243876 -0.413344725
42 2.1 5.40557664 -0.403946464
43 2.15 5.385379317 -0.393486069
44 2.2 5.365705013 -0.381983583
45 2.25 5.346605834 -0.369462178
46 2.3 5.328132725 -0.355948172
47 2.35 5.310335317 -0.341471037
48 2.4 5.293261765 -0.326063387
49 2.45 5.276958596 -0.309760957
50 2.5 5.261470548 -0.292602563
51 2.55 5.24684042 -0.27463005
52 2.6 5.233108917 -0.255888213
53 2.65 5.220314506 -0.236424712
54 2.7 5.208493271 -0.216289957
55 2.75 5.197678773 -0.19553698
56 2.8 5.187901924 -0.174221286
57 2.85 5.17919086 -0.152400686
58 2.9 5.171570825 -0.130135105
59 2.95 5.16506407 -0.107486385
60 3 5.159689751 -0.084518054
61 3.05 5.155463848 -0.061295096
62 3.1 5.152399093 -0.037883691
63 3.15 5.150504909 -0.014350957
64 3.2 5.149787361 0.009235327
65 3.25 5.150249127 0.032807006
66 3.3 5.151889478 0.056295834
67 3.35 5.154704269 0.079633767
68 3.4 5.158685958 0.102753251
69 3.45 5.16382362 0.125587513
70 3.5 5.170102996 0.148070843
71 3.55 5.177506538 0.170138866
72 3.6 5.186013481 0.191728813
73 3.65 5.195599922 0.212779765
74 3.7 5.20623891 0.233232893
75 3.75 5.217900555 0.253031678
76 3.8 5.230552139 0.272122112
77 3.85 5.244158244 0.29045288
78 3.9 5.258680888 0.30797552
79 3.95 5.274079664 0.324644568
80 4 5.290311893 0.340417672
81 4.05 5.307332776 0.355255693
82 4.1 5.325095561 0.369122783
83 4.15 5.3435517 0.381986439
84 4.2 5.362651022 0.393817545
85 4.25 5.382341899 0.404590386
86 4.3 5.402571419 0.414282651
87 4.35 5.423285551 0.422875422
88 4.4 5.444429322 0.430353139
89 4.45 5.465946979 0.436703564
90 4.5 5.487782158 0.441917721
91 4.55 5.509878044 0.445989833
92 4.6 5.532177535 0.448917252
93 4.65 5.554623398 0.450700372
94 4.7 5.577158416 0.451342545
95 4.75 5.599725544 0.450849987
96 4.8 5.622268043 0.449231684
97 4.85 5.644729627 0.446499284
98 4.9 5.667054591 0.442667006
99 4.95 5.689187942 0.437751528
100 5 5.711075518 0.431771885
101 5.05 5.732664112 0.424749368
102 5.1 5.753901581 0.416707417
103 5.15 5.774736952 0.407671518
104 5.2 5.795120527 0.397669103
105 5.25 5.815003983 0.386729448
106 5.3 5.834340455 0.374883576
107 5.35 5.853084634 0.362164158
108 5.4 5.871192842 0.348605419
109 5.45 5.888623113 0.334243046
110 5.5 5.905335265 0.319114091
111 5.55 5.921290969 0.303256889
112 5.6 5.936453814 0.286710965
113 5.65 5.950789362 0.26951695
114 5.7 5.96426521 0.251716493
115 5.75 5.976851034 0.233352184
116 5.8 5.988518644 0.214467467
117 5.85 5.999242017 0.195106564
118 5.9 6.008997345 0.175314393
119 5.95 6.017763065 0.15513649
120 6 6.025519889 0.134618935
121 6.05 6.032250836 0.113808274
122 6.1 6.03794125 0.092751445
123 6.15 6.042578822 0.071495703
124 6.2 6.046153607 0.050088545
125 6.25 6.048658034 0.028577641
126 6.3 6.050086916 0.007010759
127 6.35 6.050437454 -0.014564309
128 6.4 6.049709239 -0.036099816
129 6.45 6.047904248 -0.05754813
130 6.5 6.045026842 -0.07886181
131 6.55 6.041083751 -0.099993671
132 6.6 6.036084067 -0.120896861
133 6.65 6.030039224 -0.141524927
134 6.7 6.022962978 -0.161831886
135 6.75 6.014871384 -0.181772298
136 6.8 6.005782769 -0.201301336
137 6.85 5.995717702 -0.220374857
138 6.9 5.984698959 -0.238949472
139 6.95 5.972751486 -0.256982618
140 7 5.959902355 -0.274432631
141 7.05 5.946180723 -0.291258818
142 7.1 5.931617782 -0.307421528
143 7.15 5.916246706 -0.322882227
144 7.2 5.900102594 -0.337603573
145 7.25 5.883222416 -0.351549487
146 7.3 5.865644941 -0.364685233
147 7.35 5.84741068 -0.37697749
148 7.4 5.828561805 -0.388394434
149 7.45 5.809142084 -0.398905811
150 7.5 5.789196793 -0.408483019
151 7.55 5.768772642 -0.417099188
152 7.6 5.747917683 -0.42472926
153 7.65 5.72668122 -0.431350069
154 7.7 5.705113716 -0.436940429
155 7.75 5.683266695 -0.441481209
156 7.8 5.661192634 -0.444955424
157 7.85 5.638944863 -0.447348314
158 7.9 5.616577447 -0.448647429
159 7.95 5.594145076 -0.448842712
160 8 5.57170294 -0.447926582
161 8.05 5.549306611 -0.445894017
162 8.1 5.52701191 -0.442742629
163 8.15 5.504874779 -0.438472744
164 8.2 5.482951142 -0.433087474
165 8.25 5.461296768 -0.426592788
166 8.3 5.439967129 -0.418997578
167 8.35 5.41901725 -0.410313713
168 8.4 5.398501564 -0.400556098
169 8.45 5.378473759 -0.389742713
170 8.5 5.358986624 -0.377894648
171 8.55 5.340091891 -0.365036134
172 8.6 5.321840084 -0.351194549
173 8.65 5.304280357 -0.336400422
174 8.7 5.287460336 -0.320687421
175 8.75 5.271425965 -0.304092322
176 8.8 5.256221349 -0.286654969
177 8.85 5.2418886 -0.268418206
178 8.9 5.22846769 -0.249427807
179 8.95 5.2159963 -0.229732367
180 9 5.204509681 -0.209383198
181 9.05 5.194040521 -0.18843418
182 9.1 5.184618812 -0.166941616
183 9.15 5.176271732 -0.144964051
184 9.2 5.169023529 -0.122562081
185 9.25 5.162895425 -0.099798141
186 9.3 5.157905518 -0.076736281
187 9.35 5.154068704 -0.053441915
188 9.4 5.151396608 -0.029981574
189 9.45 5.149897529 -0.00642263
190 9.5 5.149576398 0.017166977
191 9.55 5.150434747 0.040719026
192 9.6 5.152470698 0.064165304
193 9.65 5.155678963 0.087437897
194 9.7 5.160050858 0.110469478
195 9.75 5.165574332 0.133193599
196 9.8 5.172234012 0.155544965
197 9.85 5.18001126 0.177459711
198 9.9 5.188884246 0.19887566
199 9.95 5.198828029 0.219732568
200 10 5.209814657 0.239972356
201 10.05 5.221813275 0.259539328
202 10.1 5.234790241 0.278380359
203 10.15 5.248709259 0.296445076
204 10.2 5.263531513 0.313686008
205 10.25 5.279215814 0.330058725
206 10.3 5.29571875 0.345521944
207 10.35 5.312994847 0.360037625
208 10.4 5.330996728 0.373571038
209 10.45 5.34967528 0.386090817
210 10.5 5.368979821 0.397568985
211 10.55 5.38885827 0.407980976
212 10.6 5.409257319 0.41730562
213 10.65 5.4301226 0.425525133
214 10.7 5.451398857 0.432625078
215 10.75 5.473030111 0.438594321
216 10.8 5.494959827 0.44342497
217 10.85 5.517131075 0.44711231
218 10.9 5.539486691 0.449654725
219 10.95 5.561969427 0.451053616
220 11 5.584522108 0.451313309
221 11.05 5.607087773 0.450440964
222 11.1 5.629609822 0.448446473
223 11.15 5.652032145 0.445342363
224 11.2 5.674299263 0.44114369
225 11.25 5.696356448 0.43586794
226 11.3 5.718149845 0.429534918
227 11.35 5.739626591 0.422166651
228 11.4 5.760734923 0.413787278
229 11.45 5.781424287 0.404422953
230 11.5 5.801645435 0.394101741
231 11.55 5.821350522 0.382853516
232 11.6 5.840493198 0.370709868
233 11.65 5.859028691 0.357704001
234 11.7 5.876913891 0.343870644
235 11.75 5.894107423 0.329245953
236 11.8 5.910569721 0.313867423
237 11.85 5.926263092 0.297773799
238 11.9 5.941151782 0.281004989
239 11.95 5.955202032 0.263601976
240 12 5.96838213 0.245606736
241 12.05 5.980662467 0.227062155
242 12.1 5.992015575 0.208011949
243 12.15 6.002416172 0.188500583
244 12.2 6.011841201 0.168573191
245 12.25 6.020269861 0.148275504
246 12.3 6.027683636 0.127653765
247 12.35 6.034066324 0.106754662
248 12.4 6.039404058 0.085625248
249 12.45 6.04368532 0.064312868
250 12.5 6.046900963 0.042865087
251 12.55 6.049044218 0.021329614
252 12.6 6.050110698 -0.000245764
253 12.65 6.05009841 -0.021813263
254 12.7 6.049007747 -0.043325166
255 12.75 6.046841489 -0.064733897
256 12.8 6.043604794 -0.085992092
257 12.85 6.039305189 -0.10705267
258 12.9 6.033952556 -0.127868905
259 12.95 6.02755911 -0.148394493
260 13 6.020139386 -0.168583627
261 13.05 6.011710204 -0.188391063
262 13.1 6.002290651 -0.207772196
263 13.15 5.991902041 -0.226683129
264 13.2 5.980567885 -0.24508074
265 13.25 5.968313848 -0.262922762
266 13.3 5.95516771 -0.280167848
267 13.35 5.941159318 -0.296775646
268 13.4 5.926320535 -0.312706873
269 13.45 5.910685192 -0.327923386
270 13.5 5.894289022 -0.342388259
271 13.55 5.877169609 -0.356065855
272 13.6 5.859366317 -0.368921906
273 13.65 5.840920221 -0.380923584
274 13.7 5.821874042 -0.392039583
275 13.75 5.802272063 -0.402240195
276 13.8 5.782160053 -0.411497391
277 13.85 5.761585184 -0.419784901
278 13.9 5.740595939 -0.427078293
279 13.95 5.719242024 -0.433355056
280 14 5.697574271 -0.438594685
281 14.05 5.675644537 -0.442778761
282 14.1 5.653505599 -0.445891036
283 14.15 5.631211047 -0.447917516
284 14.2 5.608815171 -0.448846548
285 14.25 5.586372844 -0.448668896
286 14.3 5.563939399 -0.44737783
287 14.35 5.541570508 -0.444969204
288 14.4 5.519322047 -0.441441533
289 14.45 5.497249971 -0.436796071
290 14.5 5.475410167 -0.431036884
291 14.55 5.453858323 -0.424170914
292 14.6 5.432649777 -0.416208048
293 14.65 5.411839375 -0.407161171
294 14.7 5.391481316 -0.397046216
295 14.75 5.371629006 -0.385882205
296 14.8 5.352334895 -0.373691283
297 14.85 5.333650331 -0.360498736
298 14.9 5.315625394 -0.346333003
299 14.95 5.298308744 -0.331225672
300 15 5.281747461 -0.315211459

Análisis gráfico de los datos obtenidos:


Método de Euler
6,2

5,8
x=y1

5,6

5,4

5,2

5
0 2 4 6 8 10 12 14 16
t
Aparentemente la gráfica que nos muestran los datos obtenidos se asemeja a una gráfica

tipo senoidal normal, no obstante, su similitud no indica que es una gráfica senoidal como

tal, debido a que al inicio de la gráfica sus datos están un poco dispersos y luego comienza

a tomar forma senoidal.

También podría gustarte