Simple Runge-Kutta solution of an ODE

J. M. Cimbala, September 2001

The equation to be solved is dy/dt = B - Ay, where in general A and B can be functions of time.
The data are plotted on a separate sheet called "plot of y vs t".

tstart = 0
deltat = 0.01 This case is for A and B constant
y(tstart) = 7.6

ti yi A(t) B(t) Di ya Da yb Db yc Dc yi+1

0 7.6 3.82 10 -19.032 7.50484 -18.66849 7.506658 -18.67543 7.413246 -18.3186 7.413269
0.01 7.413269 3.82 10 -18.31869 7.321676 -17.9688 7.323425 -17.97548 7.233514 -17.63203 7.233537
0.02 7.233537 3.82 10 -17.63211 7.145377 -17.29534 7.14706 -17.30177 7.060519 -16.97118 7.060541
0.03 7.060541 3.82 10 -16.97127 6.975685 -16.64712 6.977306 -16.65331 6.894008 -16.33511 6.894029
0.04 6.894029 3.82 10 -16.33519 6.812353 -16.02319 6.813913 -16.02915 6.733738 -15.72288 6.733758
0.05 6.733758 3.82 10 -15.72296 6.655143 -15.42265 6.656645 -15.42838 6.579474 -15.13359 6.579494
0.06 6.579494 3.82 10 -15.13367 6.503825 -14.84461 6.505271 -14.85013 6.430992 -14.56639 6.431011
0.07 6.431011 3.82 10 -14.56646 6.358179 -14.28824 6.35957 -14.29356 6.288075 -14.02045 6.288094
0.08 6.288094 3.82 10 -14.02052 6.217991 -13.75273 6.21933 -13.75784 6.150515 -13.49497 6.150533
0.09 6.150533 3.82 10 -13.49503 6.083057 -13.23728 6.084346 -13.2422 6.01811 -12.98918 6.018127
0.1 6.018127 3.82 10 -12.98925 5.953181 -12.74115 5.954421 -12.74589 5.890668 -12.50235 5.890684
0.11 5.890684 3.82 10 -12.50241 5.828172 -12.26362 5.829366 -12.26818 5.768003 -12.03377 5.768018
0.12 5.768018 3.82 10 -12.03383 5.707849 -11.80398 5.708998 -11.80837 5.649934 -11.58275 5.649949
0.13 5.649949 3.82 10 -11.58281 5.592035 -11.36157 5.593141 -11.3658 5.536291 -11.14863 5.536306
0.14 5.536306 3.82 10 -11.14869 5.480562 -10.93575 5.481627 -10.93981 5.426907 -10.73079 5.426921
0.15 5.426921 3.82 10 -10.73084 5.373267 -10.52588 5.374292 -10.5298 5.321623 -10.3286 5.321637
0.16 5.321637 3.82 10 -10.32865 5.269993 -10.13137 5.27098 -10.13514 5.220285 -9.94149 5.220298
0.17 5.220298 3.82 10 -9.941538 5.17059 -9.751655 5.17154 -9.755282 5.122745 -9.568887 5.122758
0.18 5.122758 3.82 10 -9.568934 5.074913 -9.386167 5.075827 -9.389658 5.028861 -9.210249 5.028873
0.19 5.028873 3.82 10 -9.210294 4.982821 -9.034378 4.983701 -9.037738 4.938495 -8.865053 4.938507
0.2 4.938507 3.82 10 -8.865096 4.894181 -8.695773 4.895028 -8.699007 4.851517 -8.532794 4.851528
0.21 4.851528 3.82 10 -8.532836 4.808864 -8.369859 4.809678 -8.372972 4.767798 -8.212989 4.767809
0.22 4.767809 3.82 10 -8.213029 4.726743 -8.05616 4.727528 -8.059156 4.687217 -7.905169 4.687227
0.23 4.687227 3.82 10 -7.905208 4.647701 -7.754219 4.648456 -7.757102 4.609656 -7.608887 4.609666
0.24 4.609666 3.82 10 -7.608924 4.571621 -7.463594 4.572348 -7.46637 4.535002 -7.323709 4.535012
0.25 4.535012 3.82 10 -7.323745 4.498393 -7.183861 4.499092 -7.186533 4.463146 -7.049219 4.463155
0.26 4.463155 3.82 10 -7.049254 4.427909 -6.914613 4.428582 -6.917185 4.393984 -6.785018 4.393992
0.27 4.393992 3.82 10 -6.785051 4.360067 -6.655456 4.360715 -6.657932 4.327413 -6.530718 4.327421
0.28 4.327421 3.82 10 -6.53075 4.294768 -6.406013 4.295391 -6.408395 4.263338 -6.285949 4.263346
0.29 4.263346 3.82 10 -6.28598 4.231916 -6.165918 4.232516 -6.168211 4.201663 -6.050355 4.201671
0.3 4.201671 3.82 10 -6.050384 4.171419 -5.934822 4.171997 -5.937029 4.142301 -5.82359 4.142308
0.31 4.142308 3.82 10 -5.823618 4.11319 -5.712387 4.113747 -5.714512 4.085163 -5.605324 4.085171
0.32 4.085171 3.82 10 -5.605352 4.057144 -5.498289 4.057679 -5.500334 4.030167 -5.395239 4.030174
0.33 4.030174 3.82 10 -5.395265 4.003198 -5.292216 4.003713 -5.294184 3.977232 -5.193028 3.977239
0.34 3.977239 3.82 10 -5.193053 3.951274 -5.093866 3.95177 -5.09576 3.926281 -4.998395 3.926288
0.35 3.926288 3.82 10 -4.99842 3.901296 -4.90295 3.901773 -4.904773 3.87724 -4.811057 3.877246
0.36 3.877246 3.82 10 -4.811081 3.853191 -4.719189 3.85365 -4.720944 3.830037 -4.630741 3.830043
0.37 3.830043 3.82 10 -4.630764 3.806889 -4.542316 3.807331 -4.544005 3.784603 -4.457183 3.784609
0.38 3.784609 3.82 10 -4.457205 3.762323 -4.372072 3.762748 -4.373698 3.740872 -4.290129 3.740877
0.39 3.740877 3.82 10 -4.29015 3.719426 -4.208209 3.719836 -4.209774 3.698779 -4.129337 3.698785
0.4 3.698785 3.82 10 -4.129357 3.678138 -4.050487 3.678532 -4.051993 3.658265 -3.974571 3.65827
0.41 3.65827 3.82 10 -3.974591 3.638397 -3.898676 3.638776 -3.900126 3.619269 -3.825606 3.619274
0.42 3.619274 3.82 10 -3.825625 3.600145 -3.752555 3.600511 -3.753951 3.581734 -3.682224 3.581739
0.43 3.581739 3.82 10 -3.682242 3.563328 -3.611911 3.563679 -3.613254 3.545606 -3.544216 3.545611
0.44 3.545611 3.82 10 -3.544233 3.52789 -3.476538 3.528228 -3.477831 3.510832 -3.41138 3.510837
0.45 3.510837 3.82 10 -3.411397 3.49378 -3.346239 3.494106 -3.347484 3.477362 -3.283523 3.477366
0.46 3.477366 3.82 10 -3.283539 3.460949 -3.220823 3.461262 -3.222021 3.445146 -3.160458 3.44515
0.47 3.44515 3.82 10 -3.160473 3.429348 -3.100108 3.42965 -3.101261 3.414137 -3.042005 3.414141
0.48 3.414141 3.82 10 -3.04202 3.398931 -2.983918 3.399222 -2.985027 3.384291 -2.927992 3.384295
0.49 3.384295 3.82 10 -2.928006 3.369655 -2.872082 3.369934 -2.87315 3.355563 -2.818252 3.355567
0.5 3.355567 3.82 10 -2.818266 3.341476 -2.764437 3.341745 -2.765465 3.327912 -2.712625 3.327916
0.51 3.327916 3.82 10 -2.712639 3.314353 -2.660827 3.314612 -2.661817 3.301298 -2.610957 3.301301
0.52 3.301301 3.82 10 -2.61097 3.288246 -2.5611 3.288496 -2.562053 3.275681 -2.5131 3.275684
0.53 3.275684 3.82 10 -2.513112 3.263118 -2.465112 3.263358 -2.466028 3.251023 -2.41891 3.251027
0.54 3.251027 3.82 10 -2.418922 3.238932 -2.37272 3.239163 -2.373603 3.227291 -2.32825 3.227294
0.55 3.227294 3.82 10 -2.328261 3.215652 -2.283792 3.215875 -2.284641 3.204447 -2.240988 3.20445
0.56 3.20445 3.82 10 -2.240999 3.193245 -2.198196 3.193459 -2.199014 3.18246 -2.156997 3.182463
0.57 3.182463 3.82 10 -2.157007 3.171678 -2.115809 3.171884 -2.116595 3.161297 -2.076153 3.161299
0.58 3.161299 3.82 10 -2.076164 3.150919 -2.036509 3.151117 -2.037266 3.140927 -1.99834 3.140929
0.59 3.140929 3.82 10 -1.99835 3.130938 -1.960181 3.131128 -1.96091 3.12132 -1.923443 3.121323
0.6 3.121323 3.82 10 -1.923453 3.111705 -1.886715 3.111889 -1.887416 3.102449 -1.851353 3.102451
0.61 3.102451 3.82 10 -1.851362 3.093194 -1.816001 3.093371 -1.816677 3.084284 -1.781965 3.084286
0.62 3.084286 3.82 10 -1.781974 3.075377 -1.747938 3.075547 -1.748589 3.066801 -1.715178 3.066803
0.63 3.066803 3.82 10 -1.715186 3.058227 -1.682426 3.058391 -1.683052 3.049972 -1.650894 3.049974
0.64 3.049974 3.82 10 -1.650902 3.04172 -1.61937 3.041878 -1.619972 3.033775 -1.589019 3.033777
0.65 3.033777 3.82 10 -1.589027 3.025832 -1.558676 3.025983 -1.559256 3.018184 -1.529463 3.018186
0.66 3.018186 3.82 10 -1.529471 3.010539 -1.500258 3.010685 -1.500816 3.003178 -1.47214 3.00318
0.67 3.00318 3.82 10 -1.472147 2.995819 -1.444029 2.99596 -1.444566 2.988734 -1.416964 2.988736
0.68 2.988736 3.82 10 -1.416971 2.981651 -1.389907 2.981786 -1.390424 2.974832 -1.363857 2.974833
0.69 2.974833 3.82 10 -1.363864 2.968014 -1.337814 2.968144 -1.338312 2.96145 -1.31274 2.961452
0.7 2.961452 3.82 10 -1.312747 2.954888 -1.287673 2.955014 -1.288152 2.948571 -1.263539 2.948572
0.71 2.948572 3.82 10 -1.263546 2.942254 -1.239412 2.942375 -1.239873 2.936173 -1.216183 2.936175
0.72 2.936175 3.82 10 -1.216189 2.930094 -1.192959 2.93021 -1.193403 2.924241 -1.170601 2.924242
0.73 2.924242 3.82 10 -1.170606 2.918389 -1.148248 2.918501 -1.148675 2.912756 -1.126727 2.912757
0.74 2.912757 3.82 10 -1.126732 2.907124 -1.105212 2.907231 -1.105623 2.901701 -1.084498 2.901702
0.75 2.901702 3.82 10 -1.084503 2.89628 -1.063789 2.896383 -1.064185 2.891061 -1.043851 2.891062
0.76 2.891062 3.82 10 -1.043856 2.885843 -1.023919 2.885942 -1.024299 2.880819 -1.004728 2.88082
0.77 2.88082 3.82 10 -1.004733 2.875796 -0.985543 2.875892 -0.985909 2.870961 -0.967071 2.870962
0.78 2.870962 3.82 10 -0.967076 2.866127 -0.948605 2.866219 -0.948958 2.861473 -0.930826 2.861474
0.79 2.861474 3.82 10 -0.93083 2.85682 -0.913052 2.856909 -0.913391 2.85234 -0.895939 2.852341
0.8 2.852341 3.82 10 -0.895943 2.847861 -0.878831 2.847947 -0.879158 2.84355 -0.862359 2.843551
0.81 2.843551 3.82 10 -0.862364 2.839239 -0.845893 2.839321 -0.846207 2.835089 -0.830039 2.83509
0.82 2.83509 3.82 10 -0.830043 2.830939 -0.814189 2.831019 -0.814492 2.826945 -0.798929 2.826946
0.83 2.826946 3.82 10 -0.798933 2.822951 -0.783673 2.823027 -0.783965 2.819106 -0.768986 2.819107
0.84 2.819107 3.82 10 -0.768989 2.815262 -0.754302 2.815336 -0.754582 2.811561 -0.740164 2.811562
0.85 2.811562 3.82 10 -0.740168 2.807861 -0.726031 2.807932 -0.726301 2.804299 -0.712423 2.8043
0.86 2.8043 3.82 10 -0.712427 2.800738 -0.698819 2.800806 -0.699079 2.797309 -0.685722 2.79731
0.87 2.79731 3.82 10 -0.685725 2.793882 -0.672628 2.793947 -0.672878 2.790582 -0.660021 2.790582
0.88 2.790582 3.82 10 -0.660025 2.787282 -0.647418 2.787345 -0.647659 2.784106 -0.635284 2.784107
0.89 2.784107 3.82 10 -0.635287 2.78093 -0.623153 2.780991 -0.623385 2.777873 -0.611474 2.777874
0.9 2.777874 3.82 10 -0.611477 2.774816 -0.599798 2.774875 -0.600021 2.771873 -0.588556 2.771874
0.91 2.771874 3.82 10 -0.588559 2.768931 -0.577317 2.768987 -0.577532 2.766099 -0.566497 2.766099
0.92 2.766099 3.82 10 -0.5665 2.763267 -0.55568 2.763321 -0.555887 2.760541 -0.545265 2.760541
0.93 2.760541 3.82 10 -0.545268 2.757815 -0.534853 2.757867 -0.535052 2.755191 -0.524829 2.755191
0.94 2.755191 3.82 10 -0.524831 2.752567 -0.514807 2.752617 -0.514999 2.750041 -0.505158 2.750042
0.95 2.750042 3.82 10 -0.505161 2.747516 -0.495512 2.747565 -0.495697 2.745085 -0.486225 2.745086
0.96 2.745086 3.82 10 -0.486228 2.742655 -0.476941 2.742701 -0.477118 2.740315 -0.468002 2.740315
0.97 2.740315 3.82 10 -0.468004 2.737975 -0.459065 2.73802 -0.459236 2.735723 -0.450461 2.735723
0.98 2.735723 3.82 10 -0.450464 2.733471 -0.44186 2.733514 -0.442024 2.731303 -0.433578 2.731304
0.99 2.731304 3.82 10 -0.43358 2.729136 -0.425299 2.729177 -0.425457 2.727049 -0.417328 2.72705
1 2.72705 3.82 10 -0.41733 2.724963 -0.409359 2.725003 -0.409511 2.722955 -0.401687 2.722955
8

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Exact solution:
yss = 2.617801

t y
0 7.6
0.01 7.413269
0.02 7.233537
0.03 7.060541
0.04 6.894029
0.05 6.733758
0.06 6.579494
0.07 6.431011
0.08 6.288094
0.09 6.150532
0.1 6.018127
0.11 5.890684
0.12 5.768018
0.13 5.649949
0.14 5.536306
0.15 5.426921
0.16 5.321637
0.17 5.220298
0.18 5.122758
0.19 5.028873
0.2 4.938507
0.21 4.851528
0.22 4.767809
0.23 4.687227
0.24 4.609666
0.25 4.535012
0.26 4.463155
0.27 4.393992
0.28 4.327421
0.29 4.263346
0.3 4.201671
0.31 4.142308
0.32 4.085171
0.33 4.030174
0.34 3.977239
0.35 3.926288
0.36 3.877246
0.37 3.830043
0.38 3.784609
0.39 3.740877
0.4 3.698785
0.41 3.65827
0.42 3.619273
0.43 3.581739
0.44 3.545611
0.45 3.510837
0.46 3.477366
0.47 3.44515
0.48 3.414141
0.49 3.384295
0.5 3.355567
0.51 3.327916
0.52 3.301301
0.53 3.275684
0.54 3.251027
0.55 3.227294
0.56 3.20445
0.57 3.182463
0.58 3.161299
0.59 3.140929
0.6 3.121323
0.61 3.102451
0.62 3.084286
0.63 3.066803
0.64 3.049974
0.65 3.033777
0.66 3.018186
0.67 3.00318
0.68 2.988736
0.69 2.974833
0.7 2.961452
0.71 2.948572
0.72 2.936175
0.73 2.924242
0.74 2.912757
0.75 2.901702
0.76 2.891062
0.77 2.88082
0.78 2.870962
0.79 2.861474
0.8 2.852341
0.81 2.843551
0.82 2.83509
0.83 2.826946
0.84 2.819107
0.85 2.811562
0.86 2.8043
0.87 2.79731
0.88 2.790582
0.89 2.784107
0.9 2.777874
0.91 2.771874
0.92 2.766099
0.93 2.760541
0.94 2.755191
0.95 2.750042
0.96 2.745086
0.97 2.740315
0.98 2.735723
0.99 2.731304
1 2.72705
Simple Runge-Kutta solution of an ODE
J. M. Cimbala, September 2001

The equation to be solved is dy/dt = B - Ay, where in general A and B can be functions of time.
The data are plotted on a separate sheet called "plot of y vs t".

tstart = 0
deltat = 0.01 This case is for B = constant, but A = 3.82 + 2.0sin(10pt)
y(tstart) = 7.6

ti yi A(t) B(t) Di ya Da yb Db yc
0 7.6 3.82 10 -19.032 7.50484 -18.66849 7.506658 -18.67543 7.413246
0.01 7.413269 4.438034 10 -22.90034 7.298768 -22.39218 7.301308 -22.40345 7.189235
0.02 7.189273 4.995571 10 -25.91452 7.059701 -25.26723 7.062937 -25.2834 6.936439
0.03 6.936494 5.438034 10 -27.72089 6.79789 -26.96716 6.801659 -26.98765 6.666618
0.04 6.666688 5.722113 10 -28.14754 6.52595 -27.34223 6.529977 -27.36527 6.393035
0.05 6.393114 5.82 10 -27.20793 6.257075 -26.41618 6.261034 -26.43922 6.128722
0.06 6.128801 5.722113 10 -25.06969 6.003453 -24.35244 6.007039 -24.37296 5.885072
0.07 5.885142 5.438034 10 -22.0036 5.775124 -21.40532 5.778115 -21.42159 5.670926
0.08 5.670982 4.995571 10 -18.32979 5.579333 -17.87195 5.581622 -17.88339 5.492148
0.09 5.492187 4.438034 10 -14.37451 5.420315 -14.05554 5.421909 -14.06262 5.351561
0.1 5.351585 3.82 10 -10.44305 5.29937 -10.24359 5.300367 -10.2474 5.249111
0.11 5.249124 3.201966 10 -6.807516 5.215086 -6.698529 5.215631 -6.700274 5.182121
0.12 5.182127 2.644429 10 -3.70377 5.163608 -3.654798 5.163853 -3.655445 5.145573
0.13 5.145575 2.201966 10 -1.330381 5.138923 -1.315734 5.138996 -1.315895 5.132416
0.14 5.132416 1.917887 10 0.156605 5.133199 0.155104 5.133192 0.155118 5.133968
0.15 5.133968 1.82 10 0.656179 5.137248 0.650208 5.137219 0.650262 5.14047
0.16 5.14047 1.917887 10 0.14116 5.141176 0.139806 5.141169 0.139819 5.141868
0.17 5.141868 2.201966 10 -1.322219 5.135257 -1.307661 5.13533 -1.307822 5.12879
0.18 5.12879 2.644429 10 -3.562725 5.110977 -3.515618 5.111212 -3.516241 5.093628
0.19 5.09363 3.201966 10 -6.309631 5.062082 -6.208615 5.062587 -6.210232 5.031528
0.2 5.031533 3.82 10 -9.220457 4.985431 -9.044347 4.986312 -9.04771 4.941056
0.21 4.941068 4.438034 10 -11.92863 4.881425 -11.66393 4.882748 -11.6698 4.82437
0.22 4.82439 4.995571 10 -14.10058 4.753887 -13.74838 4.755648 -13.75717 4.686818
0.23 4.686848 5.438034 10 -15.48724 4.609412 -15.06614 4.611517 -15.07759 4.536072
0.24 4.536111 5.722113 10 -15.95614 4.456331 -15.49963 4.458613 -15.51269 4.380984
0.25 4.381029 5.82 10 -15.49759 4.303541 -15.04661 4.305796 -15.05973 4.230432
0.26 4.230477 5.722113 10 -14.20727 4.159441 -13.80079 4.161473 -13.81242 4.092353
0.27 4.092393 5.438034 10 -12.25457 4.03112 -11.92137 4.032786 -11.93043 3.973088
0.28 3.973119 4.995571 10 -9.847998 3.923879 -9.602016 3.925109 -9.60816 3.877038
0.29 3.877059 4.438034 10 -7.206518 3.841026 -7.046605 3.841826 -7.050153 3.806557
0.3 3.806569 3.82 10 -4.541095 3.783864 -4.45436 3.784297 -4.456016 3.762009
0.31 3.762015 3.201966 10 -2.045843 3.751786 -2.01309 3.751949 -2.013614 3.741879
0.32 3.74188 2.644429 10 0.104861 3.742405 0.103475 3.742398 0.103493 3.742915
0.33 3.742915 2.201966 10 1.758228 3.751706 1.73887 3.75161 1.739083 3.760306
0.34 3.760305 1.917887 10 2.788159 3.774246 2.761422 3.774112 2.761679 3.787922
0.35 3.787921 1.82 10 3.105983 3.803451 3.077719 3.80331 3.077976 3.818701
0.36 3.8187 1.917887 10 2.676165 3.832081 2.650502 3.831953 2.650748 3.845208
0.37 3.845207 2.201966 10 1.532985 3.852872 1.516107 3.852787 1.516293 3.86037
0.38 3.860369 2.644429 10 -0.208474 3.859327 -0.205718 3.859341 -0.205754 3.858312
0.39 3.858312 3.201966 10 -2.354183 3.846541 -2.316493 3.846729 -2.317096 3.835141
0.4 3.835143 3.82 10 -4.650245 3.811892 -4.561426 3.812336 -4.563122 3.789512
0.41 3.789517 4.438034 10 -6.818007 3.755427 -6.666714 3.756184 -6.670071 3.722817
0.42 3.722828 4.995571 10 -8.59765 3.67984 -8.382899 3.680914 -8.388263 3.638945
0.43 3.638964 5.438034 10 -9.788809 3.59002 -9.522649 3.591351 -9.529886 3.543665
0.44 3.54369 5.722113 10 -10.27739 3.492303 -9.983351 3.493773 -9.991764 3.443772
0.45 3.443801 5.82 10 -10.04292 3.393586 -9.750672 3.395048 -9.759177 3.346209
0.46 3.346238 5.722113 10 -9.147554 3.300501 -8.885837 3.301809 -8.893325 3.257305
0.47 3.257331 5.438034 10 -7.713475 3.218763 -7.503744 3.219812 -7.509447 3.182236
0.48 3.182256 4.995571 10 -5.897183 3.15277 -5.749884 3.153506 -5.753563 3.12472
0.49 3.124733 4.438034 10 -3.86767 3.105394 -3.781846 3.105823 -3.78375 3.086895
0.5 3.086902 3.82 10 -1.791964 3.077942 -1.757738 3.078113 -1.758392 3.069318
0.51 3.06932 3.201966 10 0.172142 3.070181 0.169386 3.070167 0.16943 3.071014
0.52 3.071014 2.644429 10 1.87892 3.080409 1.854076 3.080285 1.854405 3.089558
0.53 3.089557 2.201966 10 3.1969 3.105542 3.161703 3.105366 3.16209 3.121178
0.54 3.121177 1.917887 10 4.013936 3.141246 3.975445 3.141054 3.975814 3.160935
0.55 3.160934 1.82 10 4.247101 3.182169 4.208452 3.181976 4.208804 3.203022
0.56 3.20302 1.917887 10 3.856969 3.222305 3.819983 3.22212 3.820337 3.241224
0.57 3.241223 2.201966 10 2.862938 3.255537 2.831417 3.25538 2.831765 3.26954
0.58 3.269539 2.644429 10 1.353934 3.276309 1.336032 3.276219 1.336269 3.282902
0.59 3.282901 3.201966 10 -0.511737 3.280342 -0.503545 3.280383 -0.503676 3.277864
0.6 3.277865 3.82 10 -2.521443 3.265257 -2.473283 3.265498 -2.474203 3.253123
0.61 3.253126 4.438034 10 -4.437483 3.230938 -4.339014 3.231431 -4.341199 3.209714
0.62 3.209721 4.995571 10 -6.034389 3.179549 -5.883663 3.180303 -5.887427 3.150847
0.63 3.15086 5.438034 10 -7.134483 3.115187 -6.940495 3.116157 -6.945769 3.081402
0.64 3.08142 5.722113 10 -7.632235 3.043259 -7.413872 3.044351 -7.420119 3.007219
0.65 3.00724 5.82 10 -7.502139 2.96973 -7.283827 2.970821 -7.29018 2.934339
0.66 2.93436 5.722113 10 -6.790742 2.900407 -6.596455 2.901378 -6.602014 2.86834
0.67 2.868359 5.438034 10 -5.598235 2.840368 -5.446018 2.841129 -5.450157 2.813858
0.68 2.813872 4.995571 10 -4.056896 2.793587 -3.955563 2.794094 -3.958094 2.774291
0.69 2.7743 4.438034 10 -2.312436 2.762737 -2.261123 2.762994 -2.262261 2.751677
0.7 2.751681 3.82 10 -0.511421 2.749124 -0.501653 2.749173 -0.501839 2.746662
0.71 2.746663 3.201966 10 1.205278 2.75269 1.185982 2.752593 1.186291 2.758526
0.72 2.758525 2.644429 10 2.705275 2.772051 2.669506 2.771873 2.669979 2.785225
0.73 2.785223 2.201966 10 3.867033 2.804558 3.824458 2.804345 3.824927 2.823472
0.74 2.823471 1.917887 10 4.584902 2.846395 4.540935 2.846176 4.541357 2.868884
0.75 2.868883 1.82 10 4.778633 2.892776 4.735147 2.892559 4.735543 2.916238
0.76 2.916237 1.917887 10 4.406987 2.938272 4.364726 2.938061 4.365132 2.959888
0.77 2.959887 2.201966 10 3.482429 2.977299 3.444088 2.977108 3.44451 2.994332
0.78 2.994331 2.644429 10 2.081703 3.004739 2.054179 3.004602 2.054543 3.014876
0.79 3.014875 3.201966 10 0.346473 3.016607 0.340926 3.01658 0.341015 3.018285
0.8 3.018285 3.82 10 -1.529848 3.010636 -1.500628 3.010782 -1.501186 3.003273
0.81 3.003275 4.438034 10 -3.328636 2.986632 -3.254773 2.987001 -3.256412 2.970711
0.82 2.970716 4.995571 10 -4.840423 2.946514 -4.719519 2.947119 -4.722539 2.923491
0.83 2.923501 5.438034 10 -5.898099 2.894011 -5.737729 2.894813 -5.742089 2.86608
0.84 2.866095 5.722113 10 -6.400121 2.834095 -6.21701 2.83501 -6.222249 2.803873
0.85 2.803891 5.82 10 -6.318644 2.772298 -6.134772 2.773217 -6.140122 2.74249
0.86 2.742508 5.722113 10 -5.69294 2.714043 -5.530062 2.714858 -5.534722 2.687161
0.87 2.687177 5.438034 10 -4.612958 2.664112 -4.487531 2.664739 -4.490941 2.642267
0.88 2.642279 4.995571 10 -3.19969 2.62628 -3.119769 2.62668 -3.121765 2.611061
0.89 2.611068 4.438034 10 -1.588009 2.603128 -1.552771 2.603304 -1.553553 2.595533
0.9 2.595535 3.82 10 0.085056 2.59596 0.083431 2.595952 0.083462 2.59637
0.91 2.59637 3.201966 10 1.686512 2.604802 1.659512 2.604667 1.659944 2.612969
0.92 2.612968 2.644429 10 3.090191 2.628419 3.049332 2.628214 3.049872 2.643466
0.93 2.643465 2.201966 10 4.179181 2.66436 4.133169 2.66413 4.133675 2.684801
0.94 2.6848 1.917887 10 4.850858 2.709054 4.804341 2.708821 4.804787 2.732848
0.95 2.732846 1.82 10 5.02622 2.757977 4.980482 2.757748 4.980898 2.782655
0.96 2.782654 1.917887 10 4.663185 2.80597 4.618468 2.805746 4.618896 2.828843
0.97 2.828841 2.201966 10 3.770988 2.847696 3.72947 2.847488 3.729927 2.86614
0.98 2.866139 2.644429 10 2.420698 2.878242 2.388691 2.878082 2.389114 2.89003
0.99 2.890029 3.201966 10 0.746227 2.89376 0.73428 2.8937 0.734471 2.897373
1 2.897373 3.82 10 -1.067963 2.892033 -1.047565 2.892135 -1.047955 2.886893
 Dc yi+1
-18.3186 7.413269
-21.90607 7.189273
-24.65147 6.936494
-26.25329 6.666688
-26.58167 6.393114
-25.66916 6.128801
-23.67505 5.885142
-20.83869 5.670982
-17.43641 5.492187
-13.75041 5.351585
-10.0516 5.249124
-6.592976 5.182127
-3.607104 5.145575
-1.301405 5.132416
0.15363 5.133968
0.644344 5.14047
0.138478 5.141868
-1.293421 5.12879
-3.46974 5.09363
-6.110781 5.031533
-8.874835 4.941068
-11.41072 4.82439
-13.41333 4.686848
-14.66731 4.536111
-15.06849 4.381029
-14.62111 4.230477
-13.4169 4.092393
-11.60579 3.973119
-9.368015 3.877059
-6.89363 3.806569
-4.370875 3.762015
-1.981368 3.74188
0.102124 3.742915
1.719934 3.760305
2.735193 3.787921
3.049964 3.8187
2.625326 3.845207
1.499597 3.860369
-0.203033 3.858312
-2.27999 3.835143
-4.475934 3.789517
-6.521987 3.722828
-8.178608 3.638964
-9.27057 3.54369
-9.705653 3.443801
-9.474937 3.346238
-8.638668 3.257331
-7.305109 3.182256
-5.60976 3.124733
-3.699746 3.086902
-1.724794 3.06932
0.166717 3.071014
1.829881 3.089557
3.127272 3.121177
3.937684 3.160934
4.170501 3.20302
3.783699 3.241223
2.800583 3.269539
1.318598 3.282901
-0.49561 3.277865
-2.426928 3.253126
-4.244819 3.209721
-5.740278 3.15086
-6.756769 3.08142
-7.207647 3.00724
-7.077851 2.93436
-6.412967 2.868359
-5.301854 2.813872
-3.859166 2.7743
-2.212036 2.751681
-0.492251 2.746663
1.167293 2.758525
2.634669 2.785223
3.78281 2.823471
4.497804 2.868883
4.692446 2.916237
4.323269 2.959887
3.406582 2.994331
2.027373 3.014875
0.335554 3.018285
-1.472503 3.003275
-3.184115 2.970716
-4.604505 2.923501
-5.585842 2.866095
-6.044077 2.803891
-5.961289 2.742508
-5.376237 2.687177
-4.368739 2.642279
-3.043741 2.611068
-1.519062 2.595535
0.081867 2.59637
1.633362 2.612968
3.009539 2.643465
4.088159 2.6848
4.758707 2.732846
4.935568 2.782654
4.5746 2.828841
3.688856 2.866139
2.357519 2.890029
0.722709 2.897373
-1.027931 2.886894