Euler method: Difference between revisions
→{{header|Oforth}}
(Added Erlang implementation) |
|||
Line 1,648:
100.000 20.000 20.07295
- : unit = ()
</pre>
=={{header|Oforth}}==
<lang oforth>: euler(f, y, a, b, h)
| t |
a b h step: t [
System.Out t <<wjp(6, JUSTIFY_RIGHT, 3) " : " << y << cr
t y f perform h * y + ->y
] ;</lang>
Usage :
<lang oforth>: newtonCoolingLaw(t, y)
y 20 - -0.07 * ;
: test
euler(#newtonCoolingLaw, 100.0, 0.0, 100.0, 2)
euler(#newtonCoolingLaw, 100.0, 0.0, 100.0, 5)
euler(#newtonCoolingLaw, 100.0, 0.0, 100.0, 10) ;</lang>
{{out}}
<pre>
....
0 : 100
10 : 44
20 : 27.2
30 : 22.16
40 : 20.648
50 : 20.1944
60 : 20.05832
70 : 20.017496
80 : 20.0052488
90 : 20.00157464
100 : 20.000472392
</pre>
| |||