Scribd
Upload
25 views13 pages

Euler Sistem As

This document contains the results of numerically solving several ordinary differential equations using the Euler method and Runge-Kutta fourth-order method. The solutions are presented in tables showing the iteration number, independent variable value, and the dependent variable values. Graphs are also shown comparing the solutions from each method.

Uploaded by

Diego Montenegro
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as XLSX, PDF, TXT or read online on Scribd
25 views13 pages

Euler Sistem As

This document contains the results of numerically solving several ordinary differential equations using the Euler method and Runge-Kutta fourth-order method. The solutions are presented in tables showing the iteration number, independent variable value, and the dependent variable values. Graphs are also shown comparing the solutions from each method.

Uploaded by

Diego Montenegro
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as XLSX, PDF, TXT or read online on Scribd
You are on page 1/ 13

y1'= -2(y1-y2) y1(0)=293 h= 0.

2
y2'= -(y2-y1)/2 y2(0)=400 x en [0, 1]

Iteracin x y1 y2
0 0 293 400
1 0.2 335.8 389.3
2 0.4 357.2 383.95
3 0.6 367.9 381.275
4 0.8 373.25 379.9375
5 1 375.925 379.26875

y''= (-y'/r)+5y y(0)=10 h= 0.1


y'(0)=0 x en [0 ; 0,5]

Iteracin x y1 y2
0 1.000E-009 10 0
1 0.1 10 5
2 0.2 10.5 5
3 0.3 11 7.75
4 0.4 11.775 10.666667
5 0.5 12.8416667 13.8875

y''=(-3y'/r)-2y y(0)=5 h= 0.25


y'(0)=0 [0; 1,5]

Iteracin x y1 y2
0 1.000E-009 5 0
1 0.25 5 -2.5
2 0.5 4.375 2.5
3 0.75 5 -3.4375
4 1 4.140625 -2.5
5 1.25 3.515625 -2.695313
6 1.5 2.84179688 -2.835938

y''=(2/x)+(1/2y) y(1)=1/2 h= 1
y'(1)=2 x en [1, 5]

Iteracin x y1 y2
0 1 0.5 2
1 2 2.5 5
2 3 7.5 6.2
3 4 13.7 6.9333333
4 5 20.6333333 7.4698297
u=y'
u'=y''
y1 = y y1'=y'=y2 y1'=y2
y2=y' y2'=y'' y2'=(-y2/x)+5y1

y1'=y2
u=y' y2'=(-y2/x)+5y1
u'=y''

y'=u
u'=(-u/x)+5y

y1 = y y1'=y'=y2
y2=y' y2'=y''

y1'=y2
y2'= (-3*y2/x)-2*y1

y1'=y2
y2'=(2/x)+(1/2y1)
y'=u
u'=(2/x)+(1/2y)
x y1 y2
0 293 400
390
0.2 335.8 389.3
0.4 357.2 383.95
0.6 367.9 381.275 370
0.8 373.25 379.9375
1 375.925 379.26875 350

x y1 y2
330
0 293 400
rk 0.2 326.660417 391.584896
0.4 347.084576 386.478856 310
0.6 359.47736 383.38066
0.8 366.99694 381.500765 290
1 371.559601 380.3601

270
0 0.2 0

405

400

395

390

385

380

375

370

365
0 0.2 0.4 0.6 0.8 1
390

370

350

330 Euler
RK4

310

290

270
0 0.2 0.4 0.6 0.8 1 1.2

Euler
RK4

6 0.8 1 1.2
euler
Iteracin x y1 y2
0 1.0000E-009 10 0
1 0.1 10 5
2 0.2 10.5 5.00000005
3 0.3 11 7.75000004
4 0.4 11.775 10.6666667
5 0.5 12.8416667 13.8875

rk
x y u
1.0000E-009 10 0
0.1 10.1677083 2.51041669
0.2 10.5491178 5.13933794
0.3 11.2019028 7.95698358
0.4 12.1504937 11.0746063
0.5 13.4308074 14.6128965

16

14

12

10

0
0 0.1 0.2 0.3 0.4
14

13

12

11 Euler
RK 4

10

8
0 0.1 0.2 0.3 0.4 0.5 0.6

Euler
RK4

2 0.3 0.4 0.5 0.6


Euler
RK 4

0.6
euler
x y1 y2
6
1.0000E-009 5 0
0.25 5 -2.5
0.5 4.375 2.5 5
0.75 5 -3.4375
1 4.140625 -2.5 4
1.25 3.515625 -2.6953125
1.5 2.84179688 -2.8359375 3

rk4
2
x y u
1.0000E-009 5 0
0.25 4.37825522 3.68489573 1
0.5 4.35701371 0.16112854
0.75 4.18486148 -1.2271881 0
0 0.2 0.4 0.6
1 3.78885868 -1.87489093
1.25 3.26832889 -2.25430782
1.5 2.67673854 -2.45008493 5

0
0 0.2 0.4 0.6
-1

-2

-3

-4
Euler
RK4

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6

Euler
RK4
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
euler
Iteracin x y1 y2
25
0 1 0.5 2
1 2 2.5 5
2 3 7.5 6.2 20
3 4 13.7 6.93333333
4 5 20.6333333 7.46982968
15
rk4
x y u
1 0.5 2 10
2 3.53703704 3.76574074
3 7.78965317 4.67112997
5
4 12.7893831 5.29637918
5 18.334509 5.77523793
0
0.5 1 1.5 2

0
0.5 1 1.5 2 2.5 3 3.5 4
25

20

15

Euler
RK4
10

0
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5

Euler
RK4

2.5 3 3.5 4 4.5 5 5.5

You might also like