SymPy の初心者として、次の SymPy 式を検討しています (これは 1 if に等しくなければなりません0<x<1
):
f = Sum((2/(lambda_m*besselj(1, lambda_m)))*besselj(0,x*lambda_m), (m, 1, oo))
どこlambda_m
は のm
- 番目のゼロですbesselj(0,x)
。
今、これを正確に計算mpmath
する関数があります。besseljzero(0,m,0)
残念ながら、上の式で (手動で、つまり...) を置き換えるとlambda_m
、besseljzero(0,m,0)
SymPy はm
整数ではないため、エラーが発生します...
次の関数を作成することで、この問題を解決できると思います。
- 整数に評価されない
j0(m)
場合に返すm
besseljzero(0,m,0)
あれば返す
しかし、私は進む方法がわかりません。
これは良い考えですか、誰かが私を助けてくれますか?