Euler method: Difference between revisions
Content added Content deleted
m (→version 2: added a missed cut n' paste output line.) |
m (→version 2: updated the output.) |
||
Line 2,698: | Line 2,698: | ||
$=Ti; inc= word(ss, sSize) /*the 1st value; obtain the increment.*/ |
$=Ti; inc= word(ss, sSize) /*the 1st value; obtain the increment.*/ |
||
do t=0 to Ti by inc /*step through calculations by the inc.*/ |
do t=0 to Ti by inc /*step through calculations by the inc.*/ |
||
a= format(Tr + (Ti-Tr)/exp(cc*t), |
a= format(Tr + (Ti-Tr)/exp(cc*t),6,10) /*calculate the analytic (exact) value.*/ |
||
say center(t,11) format($,6,3) 'ºC ' a "ºC" format(abs(a-$)/a*100,6,2) '%' |
say center(t,11) format($,6,3) 'ºC ' a "ºC" format(abs(a-$)/a*100,6,2) '%' |
||
$= $ + inc * cc * (Tr-$) /*calc. next value via Euler's method. */ |
$= $ + inc * cc * (Tr-$) /*calc. next value via Euler's method. */ |
||
Line 2,712: | Line 2,712: | ||
time in |
time in |
||
══seconds══ ══Euler method══ ═════analytic═════ ══difference══ |
══seconds══ ══Euler method══ ═════analytic═════ ══difference══ |
||
0 100.000 ºC 100. |
0 100.000 ºC 100.0000000000 ºC 0.00 % |
||
2 88.800 ºC 89. |
2 88.800 ºC 89.5486588319 ºC 0.84 % |
||
4 79.168 ºC 80. |
4 79.168 ºC 80.4626993165 ºC 1.61 % |
||
6 70.884 ºC 72. |
6 70.884 ºC 72.5637455852 ºC 2.31 % |
||
8 63.761 ºC 65. |
8 63.761 ºC 65.6967251079 ºC 2.95 % |
||
10 57.634 ºC 59. |
10 57.634 ºC 59.7268243033 ºC 3.50 % |
||
12 52.365 ºC 54. |
12 52.365 ºC 54.5368418743 ºC 3.98 % |
||
14 47.834 ºC 50. |
14 47.834 ºC 50.0248879081 ºC 4.38 % |
||
16 43.937 ºC 46. |
16 43.937 ºC 46.1023835698 ºC 4.70 % |
||
18 40.586 ºC 42. |
18 40.586 ºC 42.6923221200 ºC 4.93 % |
||
20 37.704 ºC 39. |
20 37.704 ºC 39.7277571153 ºC 5.09 % |
||
22 35.226 ºC 37. |
22 35.226 ºC 37.1504881142 ºC 5.18 % |
||
24 33.094 ºC 34. |
24 33.094 ºC 34.9099180832 ºC 5.20 % |
||
26 31.261 ºC 32. |
26 31.261 ºC 32.9620600747 ºC 5.16 % |
||
28 29.684 ºC 31. |
28 29.684 ºC 31.2686736737 ºC 5.07 % |
||
30 28.328 ºC 29. |
30 28.328 ºC 29.7965142602 ºC 4.93 % |
||
32 27.163 ºC 28. |
32 27.163 ºC 28.5166803503 ºC 4.75 % |
||
34 26.160 ºC 27. |
34 26.160 ºC 27.4040462008 ºC 4.54 % |
||
36 25.297 ºC 26. |
36 25.297 ºC 26.4367685400 ºC 4.31 % |
||
38 24.556 ºC 25. |
38 24.556 ºC 25.5958577396 ºC 4.06 % |
||
40 23.918 ºC 24. |
40 23.918 ºC 24.8648050100 ºC 3.81 % |
||
42 23.369 ºC 24. |
42 23.369 ºC 24.2292582991 ºC 3.55 % |
||
44 22.898 ºC 23. |
44 22.898 ºC 23.6767405319 ºC 3.29 % |
||
46 22.492 ºC 23. |
46 22.492 ºC 23.1964046609 ºC 3.04 % |
||
48 22.143 ºC 22. |
48 22.143 ºC 22.7788207156 ºC 2.79 % |
||
50 21.843 ºC 22. |
50 21.843 ºC 22.4157906738 ºC 2.55 % |
||
52 21.585 ºC 22. |
52 21.585 ºC 22.1001875173 ºC 2.33 % |
||
54 21.363 ºC 21. |
54 21.363 ºC 21.8258153140 ºC 2.12 % |
||
56 21.172 ºC 21. |
56 21.172 ºC 21.5872875795 ºC 1.92 % |
||
58 21.008 ºC 21. |
58 21.008 ºC 21.3799215292 ºC 1.74 % |
||
60 20.867 ºC 21. |
60 20.867 ºC 21.1996461456 ºC 1.57 % |
||
62 20.746 ºC 21. |
62 20.746 ºC 21.0429222563 ºC 1.41 % |
||
64 20.641 ºC 20. |
64 20.641 ºC 20.9066730524 ºC 1.27 % |
||
66 20.551 ºC 20. |
66 20.551 ºC 20.7882236849 ºC 1.14 % |
||
68 20.474 ºC 20. |
68 20.474 ºC 20.6852487518 ºC 1.02 % |
||
70 20.408 ºC 20. |
70 20.408 ºC 20.5957266457 ºC 0.91 % |
||
72 20.351 ºC 20. |
72 20.351 ºC 20.5178998655 ºC 0.81 % |
||
74 20.302 ºC 20. |
74 20.302 ºC 20.4502405132 ºC 0.73 % |
||
76 20.259 ºC 20. |
76 20.259 ºC 20.3914202980 ºC 0.65 % |
||
78 20.223 ºC 20. |
78 20.223 ºC 20.3402844596 ºC 0.58 % |
||
80 20.192 ºC 20. |
80 20.192 ºC 20.2958290973 ºC 0.51 % |
||
82 20.165 ºC 20. |
82 20.165 ºC 20.2571814620 ºC 0.45 % |
||
84 20.142 ºC 20. |
84 20.142 ºC 20.2235828220 ºC 0.40 % |
||
86 20.122 ºC 20. |
86 20.122 ºC 20.1943735676 ºC 0.36 % |
||
88 20.105 ºC 20. |
88 20.105 ºC 20.1689802617 ºC 0.32 % |
||
90 20.090 ºC 20. |
90 20.090 ºC 20.1469043822 ºC 0.28 % |
||
92 20.078 ºC 20. |
92 20.078 ºC 20.1277125344 ºC 0.25 % |
||
94 20.067 ºC 20. |
94 20.067 ºC 20.1110279436 ºC 0.22 % |
||
96 20.057 ºC 20. |
96 20.057 ºC 20.0965230571 ºC 0.19 % |
||
98 20.049 ºC 20. |
98 20.049 ºC 20.0839131146 ºC 0.17 % |
||
100 20.042 ºC 20. |
100 20.042 ºC 20.0729505572 ºC 0.15 % |
||
time in |
time in |
||
══seconds══ ══Euler method══ ═════analytic═════ ══difference══ |
══seconds══ ══Euler method══ ═════analytic═════ ══difference══ |
||
0 100.000 ºC 100. |
0 100.000 ºC 100.0000000000 ºC 0.00 % |
||
5 72.000 ºC 76. |
5 72.000 ºC 76.3750471775 ºC 5.73 % |
||
10 53.800 ºC 59. |
10 53.800 ºC 59.7268243033 ºC 9.92 % |
||
15 41.970 ºC 47. |
15 41.970 ºC 47.9950199289 ºC 12.55 % |
||
20 34.281 ºC 39. |
20 34.281 ºC 39.7277571153 ºC 13.71 % |
||
25 29.282 ºC 33. |
25 29.282 ºC 33.9019154760 ºC 13.63 % |
||
30 26.034 ºC 29. |
30 26.034 ºC 29.7965142602 ºC 12.63 % |
||
35 23.922 ºC 26. |
35 23.922 ºC 26.9034869199 ºC 11.08 % |
||
40 22.549 ºC 24. |
40 22.549 ºC 24.8648050100 ºC 9.31 % |
||
45 21.657 ºC 23. |
45 21.657 ºC 23.4281701494 ºC 7.56 % |
||
50 21.077 ºC 22. |
50 21.077 ºC 22.4157906738 ºC 5.97 % |
||
55 20.700 ºC 21. |
55 20.700 ºC 21.7023789151 ºC 4.62 % |
||
60 20.455 ºC 21. |
60 20.455 ºC 21.1996461456 ºC 3.51 % |
||
65 20.296 ºC 20. |
65 20.296 ºC 20.8453763507 ºC 2.64 % |
||
70 20.192 ºC 20. |
70 20.192 ºC 20.5957266457 ºC 1.96 % |
||
75 20.125 ºC 20. |
75 20.125 ºC 20.4198014719 ºC 1.44 % |
||
80 20.081 ºC 20. |
80 20.081 ºC 20.2958290973 ºC 1.06 % |
||
85 20.053 ºC 20. |
85 20.053 ºC 20.2084672415 ºC 0.77 % |
||
90 20.034 ºC 20. |
90 20.034 ºC 20.1469043822 ºC 0.56 % |
||
95 20.022 ºC 20. |
95 20.022 ºC 20.1035217684 ºC 0.40 % |
||
100 20.014 ºC 20. |
100 20.014 ºC 20.0729505572 ºC 0.29 % |
||
time in |
time in |
||
══seconds══ ══Euler method══ ═════analytic═════ ══difference══ |
══seconds══ ══Euler method══ ═════analytic═════ ══difference══ |
||
0 100.000 ºC 100. |
0 100.000 ºC 100.0000000000 ºC 0.00 % |
||
10 44.000 ºC 59. |
10 44.000 ºC 59.7268243033 ºC 26.33 % |
||
20 27.200 ºC 39. |
20 27.200 ºC 39.7277571153 ºC 31.53 % |
||
30 22.160 ºC 29. |
30 22.160 ºC 29.7965142602 ºC 25.63 % |
||
40 20.648 ºC 24. |
40 20.648 ºC 24.8648050100 ºC 16.96 % |
||
50 20.194 ºC 22. |
50 20.194 ºC 22.4157906738 ºC 9.91 % |
||
60 20.058 ºC 21. |
60 20.058 ºC 21.1996461456 ºC 5.38 % |
||
70 20.017 ºC 20. |
70 20.017 ºC 20.5957266457 ºC 2.81 % |
||
80 20.005 ºC 20. |
80 20.005 ºC 20.2958290973 ºC 1.43 % |
||
90 20.002 ºC 20. |
90 20.002 ºC 20.1469043822 ºC 0.72 % |
||
100 20.000 ºC 20. |
100 20.000 ºC 20.0729505572 ºC 0.36 % |
||
</pre> |
</pre> |
||