もちろん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の微分…
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。