Rungekutte

Download as xlsx, pdf, or txt
Download as xlsx, pdf, or txt
You are on page 1of 5

Runge Kutta

1. input initial condition dan time increment


Δt 0.1 h
Time (x) y
0 1
0.1

2. calculate the four intermediate δs

Time y k1 k2 k3 k4
0 1 -0.1 -0.094429 -0.095489 -0.089546 1.2
0.1 0.905103 -0.08976 -0.084234 -0.085186 -0.079572
1
0.2 0.820408 -0.079733 -0.074612 -0.075413 -0.070346
0.3 0.745386 -0.070463 -0.06587 -0.066523 -0.062051
0.8
0.4 0.679169 -0.062135 -0.058083 -0.058609 -0.054706
0.5 0.620798 -0.054765 -0.051223 -0.051644 -0.048257 0.6
0.6 0.569339 -0.048299 -0.045219 -0.045555 -0.042624
Y

0.7 0.523927 -0.042654 -0.039981 -0.04025 -0.037717 0.4


0.8 0.483788 -0.037738 -0.035422 -0.035637 -0.033447
0.9 0.448237 -0.033463 -0.031455 -0.031628 -0.029734 0.2
1 0.416676 -0.029745 -0.028004 -0.028144 -0.026503
1.1 0.388586 -0.026511 -0.025 -0.025114 -0.02369 0
0 1 2 3 4
1.2 0.363514 -0.023696 -0.022382 -0.022475 -0.021237
1.3 0.341073 -0.021242 -0.020097 -0.020173 -0.019096 time Title
1.4 0.320927 -0.019099 -0.0181 -0.018163 -0.017222
1.5 0.302785 -0.017225 -0.016351 -0.016403 -0.01558
1.6 0.2864 -0.015582 -0.014816 -0.014859 -0.014137
1.7 0.271555 -0.014139 -0.013466 -0.013501 -0.012867
1.8 0.258065 -0.012868 -0.012275 -0.012305 -0.011745
1.9 0.24577 -0.011746 -0.011222 -0.011247 -0.010752
2 0.23453 -0.010753 -0.010289 -0.01031 -0.009871
2.1 0.224226 -0.009872 -0.00946 -0.009478 -0.009087
2.2 0.214754 -0.009088 -0.008721 -0.008736 -0.008388
2.3 0.206023 -0.008388 -0.00806 -0.008073 -0.007762
2.4 0.197953 -0.007762 -0.007469 -0.00748 -0.007201
2.5 0.190476 -0.007201 -0.006937 -0.006947 -0.006696
2.6 0.183532 -0.006696 -0.006458 -0.006467 -0.00624
2.7 0.177068 -0.00624 -0.006025 -0.006033 -0.005828
2.8 0.171037 -0.005828 -0.005634 -0.00564 -0.005455
2.9 0.165399 -0.005455 -0.005278 -0.005284 -0.005115
3 0.160117 -0.005115 -0.004954 -0.004959 -0.004806
3.1 0.155158 -0.004806 -0.004659 -0.004663 -0.004523
3.2 0.150496 -0.004523 -0.004389 -0.004393 -0.004264
3.3 0.146104 -0.004264 -0.004141 -0.004145 -0.004027
3.4 0.141961 -0.004027 -0.003914 -0.003917 -0.003808
3.5 0.138044 -0.003808 -0.003704 -0.003707 -0.003607
3.6 0.134338 -0.003607 -0.003511 -0.003514 -0.003421
3.7 0.130825 -0.003422 -0.003333 -0.003335 -0.00325
3.8 0.12749 -0.00325 -0.003167 -0.00317 -0.00309
3.9 0.124321 -0.00309 -0.003014 -0.003016 -0.002942
4 0.121306 -0.002942 -0.002872 -0.002873 -0.002805
4.1 0.118433 -0.002805 -0.002739 -0.00274 -0.002677
4.2 0.115693 -0.002677 -0.002615 -0.002617 -0.002557
4.3 0.113077 -0.002557 -0.0025 -0.002501 -0.002445
4.4 0.110576 -0.002445 -0.002391 -0.002393 -0.002341
4.5 0.108184 -0.002341 -0.00229 -0.002291 -0.002243
4.6 0.105893 -0.002243 -0.002195 -0.002196 -0.002151
4.7 0.103697 -0.002151 -0.002106 -0.002107 -0.002064
4.8 0.10159 -0.002064 -0.002022 -0.002023 -0.001983
4.9 0.099567 -0.001983 -0.001943 -0.001944 -0.001906
5 0.097623 -0.001906 -0.001869 -0.00187 -0.001834
y

2 3 4 5 6
time Title
Heun's Method

x1 0.1 0.2
y1

f (x, y) =1+ y + x2

yo = y(xo) = 1, h=0,1

step 1

You might also like