第1種ベッセル関数実装した
もちろんJで。
J =: 1 : '(i.0)H.(1+m)@(_0.25&*)@*: * ^&m@-:%(!m)"_' BesselJ =: 1 : 'if. (0>m)*.4=3!:0 m do. ((|m) J)*(_1^|m)"_ else. m J end.' NB. 第1種ベッセル関数 PBesselJ=: 1 : '((%~&m) * (m BesselJ)) - ((>:m) BesselJ)' NB. BesselJの微分
青が次数0、赤が次数1、緑が次数2の第1種ベッセル関数
青が次数0の第1種ベッセル関数で、赤がその微分