The J contribution on Rosetta Code
http://rosettacode.org/wiki/Runge-Kutta_method#J for subject is a loopy
one. My goal is to change that into a more functional and J-ish one.
Here's my attempt for that particular task:
p4=: 1 :(':';'({:y)+6%~+/1 2 2 1*(x*[: u y+(*x&,@{:))/\.1 0.5 0.5,x*u y')
rk4M=: 1 :0
'Y0 a b h'=. 4{. y
(,.[:h&(u p4)@(,{:)/\.Y0,~}.)&.|. a+i.@>:&.(%&h) b-a
)
(fyp rk4 -: fyp rk4M) 1 0 10 0.001
1
Any comments/ideas before entering this as an alternative?
--
Met vriendelijke groet,
@@i = Arie Groeneveld
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
