Runge-Kutta method: Difference between revisions
→{{header|Perl 6}}: golfing a bit
(→{{header|Perl 6}}: golfing a bit) |
|||
Line 500:
<lang perl6>sub runge-kutta(&yp, \δt) {
return -> \t, \y {
map {
state $yp = 0;
}
}
my &RK = runge-kutta { $^t * sqrt($^y) }, .1;
loop ( my ($t, $y) = (0, 1); $t <= 10; ($t, $y) = RK($t, $y)) {
printf "y(%2.0f) =
if $t == $t.Int;
}</lang>▼
▲</lang>
{{out}}
<pre>y( 0) = 1.000000 ± 0.000000e+00
| |||