RosettaCodeData/Task/Euler-method/Forth/euler-method.fth
Ingy döt Net 6f050a029e update
2013-06-05 21:47:54 +00:00

14 lines
329 B
Forth

: newton-cooling-law ( f: temp -- f: temp' )
20e f- -0.07e f* ;
: euler ( f: y0 xt step end -- )
1+ 0 do
cr i . fdup f.
fdup over execute
dup s>f f* f+
dup +loop
2drop fdrop ;
100e ' newton-cooling-law 2 100 euler cr
100e ' newton-cooling-law 5 100 euler cr
100e ' newton-cooling-law 10 100 euler cr