Objetivos

General

Resolver ecuaciones diferenciales ordinarias mediante métodos numéricos.

Específicas

• Comparar valores de y(x) mediante el método de Euler, para obtener una

aproximación a los valores dados.

• Comparar valores de y(x) mediante el método de Euler mejorado, para obtener

una aproximación a los valores dados.

• Comparar valores de y(x) mediante el método de RK4, para obtener una

aproximación a los valores dados.
Marco Teórico

(Mecánicos, 2012).
(Varios., s.f.).
(Desconocido).
Solución y Resultados
En los siguientes problemas construir una tabla para comparar los valores indicados
de 𝑦(𝑥) mediante el método de Euler, el método de Euler mejorado y el método RK4.
Para obtener una aproximación a los valores indicados. Calcule redondeando a cuatro
cifras decimales. Primero use h = 0.1 y luego h = 0.05.

*Los valores aproximados se encuentran en las tablas en color rojo*

1. 𝑦 ′ = 𝑥 + 𝑦 2 𝑐𝑜𝑛 𝑦(0) = 0 𝑦(0.1), 𝑦(0.2), 𝑦(0.3), 𝑦(0.5)

1.1. Método de Euler.

1.1.1. Con h = 0.1

a) y (0.1)
x y f(x,y)
0 0 0
0.1 0 0.1

b) y (0.2)
x y f(x,y)
0 0 0
0.1 0 0.1
0.2 0.01 0.2001

c) y (0.3)
x y f(x,y)
0 0 0
0.1 0 0.1
0.2 0.01 0.2001
0.3 0.0300 0.3009006

d) y (0.5)
x y f(x,y)
0 0 0
0.1 0 0.1
0.2 0.01 0.2001
0.3 0.03001 0.3009006
0.4 0.06010006 0.403612017
0.5 0.1005 0.510092465
 Imagen 1. Primeras cuatro iteraciones con h = 0.1. Método de Euler.

1.1.2. Con h = 0.05

a) y (0.1)

x y f(x,y)
0 0 0
0.05 0 0.05
0.1 0.0025 0.10000625

b) y (0.2)

x y f(x,y)
0 0 0
0.05 0 0.05
0.1 0.0025 0.10000625
0.15 0.007500313 0.150056255
0.2 0.0150 0.200225094

c) y (0.3)

x y f(x,y)
0 0 0
0.05 0 0.05
0.1 0.0025 0.10000625
0.15 0.007500313 0.150056255
0.2 0.015003125 0.200225094
0.25 0.02501438 0.250625719
0.3 0.0375 0.301409677
 d) y (0.5)

x y f(x,y)
0 0 0
0.05 0 0.05
0.1 0.0025 0.10000625
0.15 0.007500313 0.150056255
0.2 0.015003125 0.200225094
0.25 0.02501438 0.250625719
0.3 0.037545666 0.301409677
0.35 0.05261615 0.352768459
0.4 0.070254573 0.404935705
0.45 0.090501358 0.458190496
0.5 0.1134 0.512862028

Imagen 2. Primeras cuatro iteraciones con h = 0.05. Método de Euler.

1.2. Método de Euler mejorado.

1.2.1. Con h = 0.1.

a) y (0.1)

x y y* f(Xn,Yn)
0 0
0.1 0.005 0 0

b) y (0.2)

x y y* f(Xn,Yn)
0 0
0.1 0.005 0 0
0.2 0.0200 0.0150025 0.100025
 c) y (0.3)

x y y* f(Xn,Yn)
0 0
0.1 0.005 0 0
0.2 0.020012504 0.0150025 0.100025
0.3 0.0451 0.040052554 0.2004005

d) y (0.5)

x y y* f(Xn,Yn)
0 0
0.1 0.005 0 0
0.2 0.020012504 0.0150025 0.100025
0.3 0.045112739 0.040052554 0.2004005
0.4 0.080498124 0.075316255 0.302035159
0.5 0.1266 0.121146119 0.406479948

Imagen 3. Primeras cuatro iteraciones con h = 0.1. Método de Euler mejorado.

1.2.2. Con h = 0.05.

a) y (0.1)

x y y* f(Xn,Yn)
0 0
0.05 0.00125 0 0
0.1 0.0050 0.003750078 0.050001563

b) y (0.2)

x y y* f(Xn,Yn)
0 0
0.05 0.00125 0 0
0.1 0.005000391 0.003750078 0.050001563
0.15 0.011253517 0.010001641 0.100025004
0.2 0.0200 0.018759849 0.150126642
 c) y (0.3)

x y y* f(Xn,Yn)
0 0
0.05 0.00125 0 0
0.1 0.005000391 0.003750078 0.050001563
0.15 0.011253517 0.010001641 0.100025004
0.2 0.020015481 0.018759849 0.150126642
0.25 0.03129805 0.030035512 0.200400619
0.3 0.0451 0.043847028 0.250979568

d) y (0.5)

x y y* f(Xn,Yn)
0 0
0.05 0.00125 0 0
0.1 0.005000391 0.003750078 0.050001563
0.15 0.011253517 0.010001641 0.100025004
0.2 0.020015481 0.018759849 0.150126642
0.25 0.03129805 0.030035512 0.200400619
0.3 0.045120603 0.043847028 0.250979568
0.35 0.061512168 0.060222396 0.302035869
0.4 0.080513583 0.079201355 0.353783747
0.45 0.10217985 0.100837705 0.406482437
0.5 0.1266 0.125201886 0.460440722

Imagen 4. Primeras cuatro iteraciones con h = 0.05. Método de Euler mejorado.

 1.3. Método de Runge-Kutta de cuarto orden
1.3.1. Con h = 0.1.

a) y (0.1)

x y k1 k2 k3 k4
0 0 0 0.05 0.05000625 0.10002501
0.1 0.0050 0.100025006 0.15010004 0.15015639 0.20040065

b) y (0.2)

x y k1 k2 k3 k4
0 0 0 0.05 0.05000625 0.10002501
0.1 0.005000625 0.100025006 0.15010004 0.15015639 0.20040065
0.2 0.0200 0.200400651 0.25090218 0.25106024 0.30203602

c) y (0.3)

x y k1 k2 k3 k4
0 0 0 0.05 0.05000625 0.10002501
0.1 0.005000625 0.100025006 0.15010004 0.15015639 0.20040065
0.2 0.020016267 0.200400651 0.25090218 0.25106024 0.30203602
0.3 0.0451 0.302036021 0.35362694 0.3539443 0.40648294

d) y (0.5)

x y k1 k2 k3 k4
0 0 0 0.05 0.05000625 0.10002501
0.1 0.005000625 0.100025006 0.15010004 0.15015639 0.20040065
0.2 0.020016267 0.200400651 0.25090218 0.25106024 0.30203602
0.3 0.045122292 0.302036021 0.35362694 0.3539443 0.40648294
0.4 0.08051665 0.406482931 0.46016887 0.46071744 0.51602462
0.5 0.1266 0.516024518 0.57322247 0.57410229 0.63385534

Imagen 5. Primeras cuatro iteraciones con h = 0.1. Runge-Kutta de cuarto orden.

 1.3.2. Con h = 0.05.
a) y (0.1)

x y k1 k2 k3 k4
0 0 0 0.025 0.02500039 0.05000156
0.05 0.00125002 0.050001563 0.07500625 0.07500977 0.10002501
0.1 0.0050 0.100025005 0.12505627 0.12506605 0.15012665

b) y (0.2)

x y k1 k2 k3 k4
0 0 0 0.025 0.02500039 0.05000156
0.05 0.00125002 0.050001563 0.07500625 0.07500977 0.10002501
0.1 0.005000508 0.100025005 0.12505627 0.12506605 0.15012665
0.15 0.01125381 0.150126648 0.17522521 0.17524444 0.20040064
0.2 0.0200 0.200400642 0.2256263 0.22565827 0.25097962

c) y (0.3)

x y k1 k2 k3 k4
0 0 0 0.025 0.02500039 0.05000156
0.05 0.00125002 0.050001563 0.07500625 0.07500977 0.10002501
0.1 0.005000508 0.100025005 0.12505627 0.12506605 0.15012665
0.15 0.01125381 0.150126648 0.17522521 0.17524444 0.20040064
0.2 0.020016032 0.200400642 0.2256263 0.22565827 0.25097962
0.25 0.031298943 0.250979624 0.27641176 0.27645995 0.30203599
0.3 0.0451 0.302035989 0.32777443 0.32784263 0.35378398

d) y (0.5)

x y k1 k2 k3 k4
0 0 0 0.025 0.02500039 0.05000156
0.05 0.00125002 0.050001563 0.07500625 0.07500977 0.10002501
0.1 0.005000508 0.100025005 0.12505627 0.12506605 0.15012665
0.15 0.01125381 0.150126648 0.17522521 0.17524444 0.20040064
0.2 0.020016032 0.200400642 0.2256263 0.22565827 0.25097962
0.25 0.031298943 0.250979624 0.27641176 0.27645995 0.30203599
0.3 0.045121935 0.302035989 0.32777443 0.32784263 0.35378398
0.35 0.061514053 0.353783979 0.37995034 0.38004282 0.40648286
0.4 0.080516162 0.406482852 0.43322254 0.43334422 0.46044144
0.45 0.102183311 0.460441429 0.4879264 0.48808312 0.51602439
0.5 0.1266 0.516024358 0.54445689 0.5446557 0.57366063
 Imagen 6. Primeras cuatro iteraciones con h = 0.05. Runge-Kutta de cuarto orden.

2. 𝑦 ′ = 𝑥𝑦 + √𝑦 𝑐𝑜𝑛 𝑦(0) = 1 𝑦(0.6), 𝑦(0.7), 𝑦(0.8), 𝑦(1.0)

2.1. Método de Euler.

2.1.1. Con h = 0.1.

a) y (0.6)

x y f(x,y)
0 1 1
0.1 1.1 1.158808848
0.2 1.21588088 1.345846066
0.3 1.35046549 1.56723495
0.4 1.50718899 1.830551849
0.5 1.69024417 2.145215994
0.6 1.9048 2.522991982

b) y (0.7)

x y f(x,y)
0 1 1
0.1 1.1 1.158808848
0.2 1.21588088 1.345846066
0.3 1.35046549 1.56723495
0.4 1.50718899 1.830551849
0.5 1.69024417 2.145215994
0.6 1.90476577 2.522991982
0.7 2.1571 2.978640467
 c) y (0.8)

x y f(x,y)
0 1 1
0.1 1.1 1.158808848
0.2 1.21588088 1.345846066
0.3 1.35046549 1.56723495
0.4 1.50718899 1.830551849
0.5 1.69024417 2.145215994
0.6 1.90476577 2.522991982
0.7 2.15706497 2.978640467
0.8 2.4549 3.530764522

d) y (1.0)

x y f(x,y)
0 1 1
0.1 1.1 1.158808848
0.2 1.21588088 1.345846066
0.3 1.35046549 1.56723495
0.4 1.50718899 1.830551849
0.5 1.69024417 2.145215994
0.6 1.90476577 2.522991982
0.7 2.15706497 2.978640467
0.8 2.45492902 3.530764522
0.9 2.80800547 4.202915358
1 3.2283 5.02504323

Imagen 7. Primeras cuatro iteraciones con h = 0.1. Método de Euler.

2.1.2. Con h = 0.05
a) y (0.6)

x y f(x,y)
0 1 1
0.05 1.05 1.077195077
0.1 1.10385975 1.161033278
0.15 1.16191142 1.252206662
0.2 1.22452175 1.351485459
0.25 1.29209602 1.459728024
0.3 1.36508242 1.577892146
0.35 1.44397703 1.707047915
0.4 1.52932943 1.848392367
0.45 1.62174905 2.003266182
0.5 1.72191236 2.173172759
0.55 1.83057099 2.359800001
0.6 1.9486 2.565045258

b) y (0.7)

x y f(x,y)
0 1 1
0.05 1.05 1.077195077
0.1 1.10385975 1.161033278
0.15 1.16191142 1.252206662
0.2 1.22452175 1.351485459
0.25 1.29209602 1.459728024
0.3 1.36508242 1.577892146
0.35 1.44397703 1.707047915
0.4 1.52932943 1.848392367
0.45 1.62174905 2.003266182
0.5 1.72191236 2.173172759
0.55 1.83057099 2.359800001
0.6 1.94856099 2.565045258
0.65 2.07681326 2.791043899
0.7 2.2164 3.040202084
c) y (0.8)

x y f(x,y)
0 1 1
0.05 1.05 1.077195077
0.1 1.10385975 1.161033278
0.15 1.16191142 1.252206662
0.2 1.22452175 1.351485459
0.25 1.29209602 1.459728024
0.3 1.36508242 1.577892146
0.35 1.44397703 1.707047915
0.4 1.52932943 1.848392367
0.45 1.62174905 2.003266182
0.5 1.72191236 2.173172759
0.55 1.83057099 2.359800001
0.6 1.94856099 2.565045258
0.65 2.07681326 2.791043899
0.7 2.21636545 3.040202084
0.75 2.36837556 3.315234413
0.8 2.5341 3.619207203

d) y (1.0)

x y f(x,y)
0 1 1
0.05 1.05 1.077195077
0.1 1.10385975 1.161033278
0.15 1.16191142 1.252206662
0.2 1.22452175 1.351485459
0.25 1.29209602 1.459728024
0.3 1.36508242 1.577892146
0.35 1.44397703 1.707047915
0.4 1.52932943 1.848392367
0.45 1.62174905 2.003266182
0.5 1.72191236 2.173172759
0.55 1.83057099 2.359800001
0.6 1.94856099 2.565045258
0.65 2.07681326 2.791043899
0.7 2.21636545 3.040202084
0.75 2.36837556 3.315234413
0.8 2.53413728 3.619207203
0.85 2.71509764 3.955588324
0.9 2.91287705 4.328304629
0.95 3.12929228 4.741808247
1 3.3664 5.201153171
 Imagen 8. Primeras cuatro iteraciones con h = 0.05. Método de Euler.

2.2. Método de Euler mejorado.

2.2.1. Con h = 0.1.

a) y (0.6)

x y y* f(Xn,Yn)
0 1 1.1 1
0.1 1.10794044 1.224278597 1.16338154
0.2 1.23367587 1.369420347 1.3574448
0.3 1.38060053 1.539517504 1.58916975
0.4 1.55288801 1.739618463 1.86730448
0.5 1.755691 1.975978039 2.2028704
0.6 1.9954 2.256381098 2.60982495

b) y (0.7)

x y y* f(Xn,Yn)
0 1 1.1 1
0.1 1.10794044 1.224278597 1.16338154
0.2 1.23367587 1.369420347 1.3574448
0.3 1.38060053 1.539517504 1.58916975
0.4 1.55288801 1.739618463 1.86730448
0.5 1.755691 1.975978039 2.2028704
0.6 1.9953986 2.256381098 2.60982495
0.7 2.2799 2.590563005 3.1059354
 c) y (0.8)

x y y* f(Xn,Yn)
0 1 1.1 1
0.1 1.10794044 1.224278597 1.16338154
0.2 1.23367587 1.369420347 1.3574448
0.3 1.38060053 1.539517504 1.58916975
0.4 1.55288801 1.739618463 1.86730448
0.5 1.755691 1.975978039 2.2028704
0.6 1.9953986 2.256381098 2.60982495
0.7 2.27996946 2.590563005 3.1059354
0.8 2.6194 2.990758597 3.71393711

d) y (1.0)

x y y* f(Xn,Yn)
0 1 1.1 1
0.1 1.10794044 1.224278597 1.16338154
0.2 1.23367587 1.369420347 1.3574448
0.3 1.38060053 1.539517504 1.58916975
0.4 1.55288801 1.739618463 1.86730448
0.5 1.755691 1.975978039 2.2028704
0.6 1.9953986 2.256381098 2.60982495
0.7 2.27996946 2.590563005 3.1059354
0.8 2.61936489 2.990758597 3.71393711
0.9 3.02611493 3.472422591 4.46307664
1 3.5161 4.055179938 5.39117863

Imagen 9. Primeras cuatro iteraciones con h = 0.1. Método de Euler mejorado.

2.2.2. Con h = 0.05.
a) y (0.6)

x y y* f(Xn,Yn)
0 1 1.05 1
0.05 1.05192988 1.10584152 1.07823282
0.1 1.10794005 1.16610912 1.16338132
0.15 1.16839413 1.23120322 1.25618193
0.2 1.23369459 1.30156744 1.35745697
0.25 1.30428738 1.37769367 1.46812586
0.3 1.38066702 1.46012792 1.58921799
0.35 1.46338253 1.5494769 1.72188737
0.4 1.55304398 1.64641545 1.86742944
0.45 1.65033007 1.75169508 2.02730026
0.5 1.75599667 1.8661536 2.20313858
0.55 1.87088656 1.99072611 2.39679116
0.6 1.9959 2.12645759 2.61034188

b) y (0.7)

x y y* f(Xn,Yn)
0 1 1.05 1
0.05 1.05192988 1.10584152 1.07823282
0.1 1.10794005 1.16610912 1.16338132
0.15 1.16839413 1.23120322 1.25618193
0.2 1.23369459 1.30156744 1.35745697
0.25 1.30428738 1.37769367 1.46812586
0.3 1.38066702 1.46012792 1.58921799
0.35 1.46338253 1.5494769 1.72188737
0.4 1.55304398 1.64641545 1.86742944
0.45 1.65033007 1.75169508 2.02730026
0.5 1.75599667 1.8661536 2.20313858
0.55 1.87088656 1.99072611 2.39679116
0.6 1.9959405 2.12645759 2.61034188
0.65 2.13220993 2.27451719 2.84614532
0.7 2.2809 2.43621464 3.10686536
c) y (0.8)

x y y* f(Xn,Yn)
0 1 1.05 1
0.05 1.05192988 1.10584152 1.07823282
0.1 1.10794005 1.16610912 1.16338132
0.15 1.16839413 1.23120322 1.25618193
0.2 1.23369459 1.30156744 1.35745697
0.25 1.30428738 1.37769367 1.46812586
0.3 1.38066702 1.46012792 1.58921799
0.35 1.46338253 1.5494769 1.72188737
0.4 1.55304398 1.64641545 1.86742944
0.45 1.65033007 1.75169508 2.02730026
0.5 1.75599667 1.8661536 2.20313858
0.55 1.87088656 1.99072611 2.39679116
0.6 1.9959405 2.12645759 2.61034188
0.65 2.13220993 2.27451719 2.84614532
0.7 2.28087137 2.43621464 3.10686536
0.75 2.44324297 2.61301897 3.39551987
0.8 2.6208 2.80658005 3.71553232

d) y (1.0)

x y y* f(Xn,Yn)
0 1 1.05 1
0.05 1.05192988 1.10584152 1.07823282
0.1 1.10794005 1.16610912 1.16338132
0.15 1.16839413 1.23120322 1.25618193
0.2 1.23369459 1.30156744 1.35745697
0.25 1.30428738 1.37769367 1.46812586
0.3 1.38066702 1.46012792 1.58921799
0.35 1.46338253 1.5494769 1.72188737
0.4 1.55304398 1.64641545 1.86742944
0.45 1.65033007 1.75169508 2.02730026
0.5 1.75599667 1.8661536 2.20313858
0.55 1.87088656 1.99072611 2.39679116
0.6 1.9959405 2.12645759 2.61034188
0.65 2.13220993 2.27451719 2.84614532
0.7 2.28087137 2.43621464 3.10686536
0.75 2.44324297 2.61301897 3.39551987
0.8 2.62080344 2.80658005 3.71553232
0.85 2.8152137 3.01875327 4.0707915
0.9 3.02834183 3.25162787 4.46572081
0.95 3.26229169 3.50755962 4.90535863
1 3.5194 3.78920844 5.3954518
Imagen 10. Primeras cuatro iteraciones con h = 0.05. Método de Euler mejorado.

2.3. Método de Runge-Kutta de cuarto orden.

2.3.1. Con h = 0.1.

a) y (0.6)

x y k1 k2 k3 k4
0 1 1 1.07719508 1.0792697 1.1633738
0.1 1.10793839 1.163380362 1.25478064 1.25758009 1.35745815
0.2 1.23369772 1.357459006 1.46625667 1.46999826 1.58924014
0.3 1.38068454 1.589230697 1.71941617 1.72438484 1.86749277
0.4 1.55308996 1.867466287 2.02405433 2.0306248 2.20327524
0.5 1.75609163 2.203221887 2.3925476 2.40121437 2.61060188
0.6 1.9961 2.61050747 2.84061585 2.85203385 3.10732531

b) y (0.7)

x y k1 k2 k3 k4
0 1 1 1.07719508 1.0792697 1.1633738
0.1 1.10793839 1.163380362 1.25478064 1.25758009 1.35745815
0.2 1.23369772 1.357459006 1.46625667 1.46999826 1.58924014
0.3 1.38068454 1.589230697 1.71941617 1.72438484 1.86749277
0.4 1.55308996 1.867466287 2.02405433 2.0306248 2.20327524
0.5 1.75609163 2.203221887 2.3925476 2.40121437 2.61060188
0.6 1.99611409 2.61050747 2.84061585 2.85203385 3.10732531
0.7 2.2812 3.10716944 3.38833072 3.40337087 3.71630844
 c) y (0.8)

x y k1 k2 k3 k4
0 1 1 1.07719508 1.0792697 1.1633738
0.1 1.10793839 1.163380362 1.25478064 1.25758009 1.35745815
0.2 1.23369772 1.357459006 1.46625667 1.46999826 1.58924014
0.3 1.38068454 1.589230697 1.71941617 1.72438484 1.86749277
0.4 1.55308996 1.867466287 2.02405433 2.0306248 2.20327524
0.5 1.75609163 2.203221887 2.3925476 2.40121437 2.61060188
0.6 1.99611409 2.61050747 2.84061585 2.85203385 3.10732531
0.7 2.28116629 3.10716944 3.38833072 3.40337087 3.71630844
0.8 2.6213 3.716061833 4.06145694 4.08128214 4.46698808

d) y (1.0)

x y k1 k2 k3 k4
0 1 1 1.07719508 1.0792697 1.1633738
0.1 1.10793839 1.163380362 1.25478064 1.25758009 1.35745815
0.2 1.23369772 1.357459006 1.46625667 1.46999826 1.58924014
0.3 1.38068454 1.589230697 1.71941617 1.72438484 1.86749277
0.4 1.55308996 1.867466287 2.02405433 2.0306248 2.20327524
0.5 1.75609163 2.203221887 2.3925476 2.40121437 2.61060188
0.6 1.99611409 2.61050747 2.84061585 2.85203385 3.10732531
0.7 2.28116629 3.10716944 3.38833072 3.40337087 3.71630844
0.8 2.62128098 3.716061833 4.06145694 4.08128214 4.46698808
0.9 3.02908978 4.466608844 4.89324597 4.91941576 5.39747249
1 3.5206 5.396900687 5.92684749 5.96146285 6.55737045

Imagen 11. Primeras cuatro iteraciones con h = 0.1. Runge-Kutta de cuarto orden.
2.3.2. Con h = 0.05.
a) y (0.6)

x y k1 k2 k3 k4
0 1 1 1.03804784 1.03854127 1.07823131
0.05 1.05192841 1.078232035 1.11960983 1.12018525 1.16337995
0.1 1.10793843 1.163380386 1.20844065 1.20910955 1.2561818
0.15 1.16839395 1.256181825 1.30531788 1.30609344 1.35745959
0.2 1.23369782 1.357459069 1.41111038 1.41200767 1.46813331
0.25 1.30429639 1.46813206 1.52679033 1.5278266 1.589233
0.3 1.38068471 1.589230824 1.65344688 1.65464196 1.72191354
0.35 1.4634124 1.721910173 1.79230194 1.79367855 1.86747138
0.4 1.55309025 1.867466518 1.94472807 1.9463123 2.02736387
0.45 1.65039784 2.02735714 2.1122691 2.11409097 2.20323133
0.5 1.75609208 2.203222288 2.29666378 2.29875782 2.39692246
0.55 1.87101698 2.396910569 2.49987293 2.50227891 2.61052353
0.6 1.9961 2.610508147 2.72411074 2.7268745 2.84639202

b) y (0.7)

x y k1 k2 k3 k4
0 1 1 1.03804784 1.03854127 1.07823131
0.05 1.05192841 1.078232035 1.11960983 1.12018525 1.16337995
0.1 1.10793843 1.163380386 1.20844065 1.20910955 1.2561818
0.15 1.16839395 1.256181825 1.30531788 1.30609344 1.35745959
0.2 1.23369782 1.357459069 1.41111038 1.41200767 1.46813331
0.25 1.30429639 1.46813206 1.52679033 1.5278266 1.589233
0.3 1.38068471 1.589230824 1.65344688 1.65464196 1.72191354
0.35 1.4634124 1.721910173 1.79230194 1.79367855 1.86747138
0.4 1.55309025 1.867466518 1.94472807 1.9463123 2.02736387
0.45 1.65039784 2.02735714 2.1122691 2.11409097 2.20323133
0.5 1.75609208 2.203222288 2.29666378 2.29875782 2.39692246
0.55 1.87101698 2.396910569 2.49987293 2.50227891 2.61052353
0.6 1.9961148 2.610508147 2.72411074 2.7268745 2.84639202
0.65 2.13243872 2.846372373 2.97188077 2.9750552 3.10719541
0.7 2.2812 3.107170566 3.24601754 3.24966377 3.39595595
 c) y (0.8)

x y k1 k2 k3 k4
0 1 1 1.03804784 1.03854127 1.07823131
0.05 1.05192841 1.078232035 1.11960983 1.12018525 1.16337995
0.1 1.10793843 1.163380386 1.20844065 1.20910955 1.2561818
0.15 1.16839395 1.256181825 1.30531788 1.30609344 1.35745959
0.2 1.23369782 1.357459069 1.41111038 1.41200767 1.46813331
0.25 1.30429639 1.46813206 1.52679033 1.5278266 1.589233
0.3 1.38068471 1.589230824 1.65344688 1.65464196 1.72191354
0.35 1.4634124 1.721910173 1.79230194 1.79367855 1.86747138
0.4 1.55309025 1.867466518 1.94472807 1.9463123 2.02736387
0.45 1.65039784 2.02735714 2.1122691 2.11409097 2.20323133
0.5 1.75609208 2.203222288 2.29666378 2.29875782 2.39692246
0.55 1.87101698 2.396910569 2.49987293 2.50227891 2.61052353
0.6 1.9961148 2.610508147 2.72411074 2.7268745 2.84639202
0.65 2.13243872 2.846372373 2.97188077 2.9750552 3.10719541
0.7 2.28116738 3.107170566 3.24601754 3.24966377 3.39595595
0.75 2.44362146 3.395924802 3.54973461 3.55392332 3.71610247
0.8 2.6213 3.716063693 3.88668015 3.89149313 4.07153035

d) y (1.0)

x y k1 k2 k3 k4
0 1 1 1.03804784 1.03854127 1.07823131
0.05 1.05192841 1.078232035 1.11960983 1.12018525 1.16337995
0.1 1.10793843 1.163380386 1.20844065 1.20910955 1.2561818
0.15 1.16839395 1.256181825 1.30531788 1.30609344 1.35745959
0.2 1.23369782 1.357459069 1.41111038 1.41200767 1.46813331
0.25 1.30429639 1.46813206 1.52679033 1.5278266 1.589233
0.3 1.38068471 1.589230824 1.65344688 1.65464196 1.72191354
0.35 1.4634124 1.721910173 1.79230194 1.79367855 1.86747138
0.4 1.55309025 1.867466518 1.94472807 1.9463123 2.02736387
0.45 1.65039784 2.02735714 2.1122691 2.11409097 2.20323133
0.5 1.75609208 2.203222288 2.29666378 2.29875782 2.39692246
0.55 1.87101698 2.396910569 2.49987293 2.50227891 2.61052353
0.6 1.9961148 2.610508147 2.72411074 2.7268745 2.84639202
0.65 2.13243872 2.846372373 2.97188077 2.9750552 3.10719541
0.7 2.28116738 3.107170566 3.24601754 3.24966377 3.39595595
0.75 2.44362146 3.395924802 3.54973461 3.55392332 3.71610247
0.8 2.62128265 3.716063693 3.88668015 3.89149313 4.07153035
0.85 2.81581549 4.071482348 4.26100134 4.26653342 4.46667104
0.9 3.02909235 4.466611897 4.67741911 4.68378029 4.90657276
0.95 3.26322221 4.906500198 5.14131497 5.14863303 5.39699441
1 3.5206 5.39690568 5.65883198 5.66725551 5.94451479
Imagen 12. Primeras cuatro iteraciones con h = 0.05. Runge-Kutta de cuarto orden.

3. 𝑦 ′ = 𝑦 − 𝑦 2 𝑐𝑜𝑛 𝑦(0) = 0.5 𝑦(0.1), 𝑦(0.2), 𝑦(0.3), 𝑦(0.5)

3.1. Método de Euler.

3.1.1. Con h = 0.1.

a) y (0.1)

x y f(x,y)
0 0.5 0.25
0.1 0.525 0.249375

b) y (0.2)

x y f(x,y)
0 0.5 0.25
0.1 0.525 0.249375
0.2 0.5499 0.24750625

c) y (0.3)

x y f(x,y)
0 0.5 0.25
0.1 0.525 0.249375
0.2 0.5499375 0.24750625
0.3 0.5747 0.24442168

d) y (0.5)

x y f(x,y)
0 0.5 0.25
0.1 0.525 0.249375
0.2 0.5499375 0.24750625
0.3 0.57468812 0.24442168
0.4 0.59913029 0.24017319
0.5 0.6231 0.23483467
 Imagen 13. Primeras cuatro iteraciones con h = 0.1. Método de Euler.

3.1.2. Con h = 0.05.

a) y (0.1)

x y f(x,y)
0 0.5 0.25
0.05 0.5125 0.24984375
0.1 0.5249 0.24937539

b) y (0.2)

x y f(x,y)
0 0.5 0.25
0.05 0.5125 0.24984375
0.1 0.52499219 0.24937539
0.15 0.53746096 0.24859668
0.2 0.5499 0.24751091

c) y (0.3)

x y f(x,y)
0 0.5 0.25
0.05 0.5125 0.24984375
0.1 0.52499219 0.24937539
0.15 0.53746096 0.24859668
0.2 0.54989079 0.24751091
0.25 0.56226634 0.2461229
0.3 0.5746 0.24443895
 d) y (0.5)

x y f(x,y)
0 0.5 0.25
0.05 0.5125 0.24984375
0.1 0.52499219 0.24937539
0.15 0.53746096 0.24859668
0.2 0.54989079 0.24751091
0.25 0.56226634 0.2461229
0.3 0.57457248 0.24443895
0.35 0.58679443 0.24246673
0.4 0.59891777 0.24021528
0.45 0.61092853 0.23769486
0.5 0.6228 0.2349169

Imagen 14. Primeras cuatro iteraciones con h = 0.05. Método de Euler.

3.2. Método de Euler mejorado.

3.2.1. Con h = 0.1.

a) y (0.1)

x y y* f(Xn,Yn)
0 0.5 0.525 0.25
0.1 0.5249 0.54990641 0.24937656

b) y (0.2)

x y y* f(Xn,Yn)
0 0.5 0.525 0.25
0.1 0.52496875 0.54990641 0.24937656
0.2 0.5498 0.57456491 0.24751866
 c) y (0.3)

x y y* f(Xn,Yn)
0 0.5 0.525 0.25
0.1 0.52496875 0.54990641 0.24937656
0.2 0.54981305 0.57456491 0.24751866
0.3 0.5744 0.59885728 0.24446301

d) y (0.5)

x y y* f(Xn,Yn)
0 0.5 0.525 0.25
0.1 0.52496875 0.54990641 0.24937656
0.2 0.54981305 0.57456491 0.24751866
0.3 0.57441098 0.59885728 0.24446301
0.4 0.59864549 0.6226724 0.24026907
0.5 0.6224 0.64590819 0.23501664

Imagen 15. Primeras cuatro iteraciones con h = 0.1. Método de Euler mejorado.

3.2.2. Con h = 0.05.

a) y (0.1)

x y y* f(Xn,Yn)
0 0.5 0.5125 0.25
0.05 0.51249609 0.52498829 0.24984385
0.1 0.5249 0.53744539 0.24937617

b) y (0.2)

x y y* f(Xn,Yn)
0 0.5 0.5125 0.25
0.05 0.51249609 0.52498829 0.24984385
0.1 0.52497658 0.53744539 0.24937617
0.15 0.53742593 0.54985589 0.2485993
0.2 0.5498 0.56220463 0.24751709
 c) y (0.3)

x y y* f(Xn,Yn)
0 0.5 0.5125 0.25
0.05 0.51249609 0.52498829 0.24984385
0.1 0.52497658 0.53744539 0.24937617
0.15 0.53742593 0.54985589 0.2485993
0.2 0.54982877 0.56220463 0.24751709
0.25 0.56216996 0.57447671 0.2461349
0.3 0.5744 0.58665764 0.24445948

d) y (0.5)

x y y* f(Xn,Yn)
0 0.5 0.5125 0.25
0.05 0.51249609 0.52498829 0.24984385
0.1 0.52497658 0.53744539 0.24937617
0.15 0.53742593 0.54985589 0.2485993
0.2 0.54982877 0.56220463 0.24751709
0.25 0.56216996 0.57447671 0.2461349
0.3 0.57443467 0.58665764 0.24445948
0.35 0.58660842 0.59873336 0.24249898
0.4 0.59867718 0.61069032 0.24026281
0.45 0.61062744 0.62251552 0.23776157
0.5 0.6224 0.63419658 0.23500692

Imagen 16. Primeras cuatro iteraciones con h = 0.05. Método de Euler mejorado.

3.3. Método de Runge-Kutta de cuarto orden.

3.3.1. Con h = 0.1.

a) y (0.1)

x y k1 k2 k3 k4
0 0.5 0.25 0.24984375 0.24984395 0.24937578
0.1 0.5249 0.24937604 0.24859765 0.24860056 0.24751605
 b) y (0.2)

x y k1 k2 k3 k4
0 0.5 0.25 0.24984375 0.24984395 0.24937578
0.1 0.52497919 0.24937604 0.24859765 0.24860056 0.24751605
0.2 0.5498 0.24751657 0.24612994 0.24613856 0.24445752

c) y (0.3)

x y k1 k2 k3 k4
0 0.5 0.25 0.24984375 0.24984395 0.24937578
0.1 0.52497919 0.24937604 0.24859765 0.24860056 0.24751605
0.2 0.54983399 0.24751657 0.24612994 0.24613856 0.24445752
0.3 0.5744 0.24445831 0.2424891 0.24250616 0.24025967

d) y (0.5)

x y k1 k2 k3 k4
0 0.5 0.25 0.24984375 0.24984395 0.24937578
0.1 0.52497919 0.24937604 0.24859765 0.24860056 0.24751605
0.2 0.54983399 0.24751657 0.24612994 0.24613856 0.24445752
0.3 0.57444251 0.24445831 0.2424891 0.24250616 0.24025967
0.4 0.59868765 0.24026075 0.23774536 0.23777319 0.23500233
0.5 0.6225 0.23500371 0.23198781 0.23202826 0.22878254

Imagen 17. Primeras cuatro iteraciones con h = 0.1. Runge-Kutta de cuarto orden.

3.3.2. Con h = 0.05.

a) y (0.1)

x y k1 k2 k3 k4
0 0.5 0.25 0.24996094 0.24996095 0.2498438
0.05 0.5124974 0.24984382 0.24964868 0.24964886 0.24937601
0.1 0.5249 0.24937604 0.24902571 0.24902626 0.24859896
 b) y (0.2)

x y k1 k2 k3 k4
0 0.5 0.25 0.24996094 0.24996095 0.2498438
0.05 0.5124974 0.24984382 0.24964868 0.24964886 0.24937601
0.1 0.52497919 0.24937604 0.24902571 0.24902626 0.24859896
0.15 0.53742985 0.24859901 0.24809513 0.24809623 0.24751651
0.2 0.5498 0.24751657 0.24686155 0.24686338 0.246134

c) y (0.3)

x y k1 k2 k3 k4
0 0.5 0.25 0.24996094 0.24996095 0.2498438
0.05 0.5124974 0.24984382 0.24964868 0.24964886 0.24937601
0.1 0.52497919 0.24937604 0.24902571 0.24902626 0.24859896
0.15 0.53742985 0.24859901 0.24809513 0.24809623 0.24751651
0.2 0.549834 0.24751657 0.24686155 0.24686338 0.246134
0.25 0.5621765 0.24613408 0.24533103 0.24533377 0.24445821
0.3 0.5744 0.24445831 0.24351106 0.24351487 0.24249728

d) y (0.5)

x y k1 k2 k3 k4
0 0.5 0.25 0.24996094 0.24996095 0.2498438
0.05 0.5124974 0.24984382 0.24964868 0.24964886 0.24937601
0.1 0.52497919 0.24937604 0.24902571 0.24902626 0.24859896
0.15 0.53742985 0.24859901 0.24809513 0.24809623 0.24751651
0.2 0.549834 0.24751657 0.24686155 0.24686338 0.246134
0.25 0.5621765 0.24613408 0.24533103 0.24533377 0.24445821
0.3 0.57444252 0.24445831 0.24351106 0.24351487 0.24249728
0.35 0.58661758 0.2424974 0.24141042 0.24141545 0.24026061
0.4 0.59868766 0.24026075 0.23903913 0.23904552 0.2377588
0.45 0.61063923 0.23775896 0.23640836 0.23641623 0.23500354
0.5 0.6225 0.23500371 0.23353028 0.23353973 0.23200745

Imagen 18. Primeras cuatro iteraciones con h = 0.05. Runge-Kutta de cuarto orden.
Conclusiones

• La solución de la ecuación diferencial 𝑦 ′ = 𝑥 + 𝑦 2 𝑐𝑜𝑛 𝑦(0) = 0, en el punto y

𝑦(0.1) es de 0.005.

• La solución de la ecuación diferencial 𝑦 ′ = 𝑥𝑦 + √𝑦 𝑐𝑜𝑛 𝑦(0) = 1, en el punto y

𝑦(0.6) es de 1.9961.

• La solución de la ecuación diferencial 𝑦 ′ = 𝑦 − 𝑦 2 𝑐𝑜𝑛 𝑦(0) = 0.5, en el punto y

𝑦(0.1) es de 0.5249.

• Al utilizar valor de h muy pequeños, se obtienen valores más cercanos al valor

teórico.

• El método de Runge-Kutta de cuarto orden resulta más eficiente que los otros
métodos, ya que, se utiliza el método de integración de Simpson. Al combinar
el método con valor de h muy pequeños, el valor que se obtiene posee un error
despreciable.

Comentario

Los métodos numéricos para resolver ecuaciones diferenciales ordinarias pueden ser
de gran ayuda en diferentes áreas. La primera opción seria para comparar resultados
obtenidos, con los métodos numéricos se pueden obtener resultados muy cercanos al
valor real y así poder comprobar si el resultado al que se llego es correcto o no.

Otra opción seria cuando no se pueden resolver ecuaciones diferenciales, los métodos
numéricos nos podrían ayudar a visualizar el resultado de la ecuación diferencial y
analizar hacia que valor tiende y llegar a un posible resultado con errores
despreciables.

También hay que recordar que el método de integración de Simpson es uno de los
más exactos, y es curioso ver que al usar e este método en otros métodos (Runge-
Kutta de cuarto orden) el margen de error se sigue manteniendo, ya que el método de
Runge-Kutta de cuarto orden es uno de los más exactos para resolver ecuaciones
diferenciales.
Bibliografía
Desconocido. (s.f.). Capmus USAL. Obtenido de Capmus USAL:
http://campus.usal.es/~mpg/Personales/PersonalMAGL/Docencia/MetNumTema4Teo(09-
10).pdf

Mecánicos, I. (12 de Mayo de 2012). Google Sites. Obtenido de Google Sites:

https://sites.google.com/site/metodosnumericosmecanica/home/unidad-vi/62-mtodos-de-
un-paso-mtodo-de-euler-mtodo-de-euler-mejorado-y-mtodo-de-runge-kutta

Varios. (s.f.). Wikipedia. Obtenido de Wikipedia :

https://es.wikipedia.org/wiki/M%C3%A9todo_de_Euler#:~:text=La%20soluci%C3%B3n%20de
%20las%20ecuaciones,un%20error%20propio%20del%20m%C3%A9todo.