Bessel関数Jν・Yνの正の零点を次数と順位から計算

使い方

1. 通常Bessel関数の実数次数νを0以上10以下で指定し、JνとYνそれぞれについて求めたい正側の順位sを1以上10以下から入力します。
2. 計算後は円筒BesselのJν零点とYν零点を別々の結果欄で読み、同じsでも二つの値が一致しないことを確認します。
3. 境界値問題へ転記するときは、通常型であること、次数ν、各関数内の順位s、原点を除く正の実零点という条件を残します。

計算式・考え方

Jν(x)=0の正側s番目をj(ν,s)、Yν(x)=0の正側s番目をy(ν,s)として返します。NIST DLMF §10.21(i) https://dlmf.nist.gov/10.21.i の単純零点と交錯、§10.74(vi) https://dlmf.nist.gov/10.74.vi の零点計算指針を照合しました。

注意事項

• 通常Bessel関数は円筒座標の波動方程式や境界値問題で現れます。半径方向の境界でJνまたはYνを0にする条件を置くと、許される波数が正の零点で決まります。どちらの関数を使うかは原点を含む領域か外側領域かで変わります。
• JνとYνの零点は正の実軸上で交互に現れますが、同じ順位sを一つの共通列として数えません。Jの符号変化とYの符号変化を別々に数え、それぞれs番目の区間だけを有界反復で狭めます。
• ν>0ではJν(0)=0と書ける場合がありますが、このページが求めるのは厳密に正のxです。原点を第1零点へ含めると文献のj(ν,s)と順位が一つずれるため、探索は小さな正数から開始します。Yνは原点で特異です。
• 公開域0≤ν≤10、s≤10では対象根がx<50へ収まることを全格子oracleで確認しています。canonicalの大次数や高順位を無断で近似せず、遷移領域と一様漸近展開の追加監査が終わるまで範囲を広げません。
• ν=1、s=1ではJ零点が約3.831705970207512、Y零点が約2.197141326031017です。ν=0ならJの第1零点は約2.404825557695773、Yは約0.893576966279168となり、よく使われる基本検算になります。
• 探索はx=0を避け、有限な走査点の符号変化から順位を確定した後、bracketを保持したまま最大64回で狭めます。関数値の残差が小さいだけで合格せず、根を挟む区間幅と有限値を同時に確認します。
• DLMF §10.6 https://dlmf.nist.gov/10.6 の導関数漸化式はbracket内の補助検算に使えますが、Newton反復だけを無制限に行いません。初期値が悪いと別順位へ飛ぶため、符号変化の計数を順位契約の正本にします。
• 円板や円柱の内部で原点を含む問題では、原点で有限なJνを採用し、外周x=kRでJν(kR)=0を置くDirichlet型境界が代表例です。一方のYνは原点で特異なので、内半径を持つ環状領域や外部解で現れます。この計算器は零点だけを返し、半径R、波数k、境界条件の選択を自動で補いません。円筒殻の内外二境界を同時に課す場合はJνとYνの線形結合が対象となり、どちらか一方の零点表だけでは固有条件を満たしません。
• 単純零点の近くではJνまたはYνの符号が根をまたいで変わります。表示値を関数へ再代入する検算では、値が小さいことに加えて左右の点が反対符号であることを確かめると、停留点や別順位を根と誤認しにくくなります。微小な再代入残差をそのまま根の絶対誤差とはみなしません。さらに同じ関数のs−1番目とs+1番目が表示根の両側に並ぶことを確認すると、根そのものだけでなく順位の数え方も検査できます。
• 次数νは整数表の行番号ではなく連続な実数入力です。ν=0とν=1の零点表からν=0.5の値を線形補間しても、関数そのものの次数依存を再現できません。本ページは指定されたνでJνとYνを直接評価してbracketを作るため、既成表の丸め値を補間する方法とは結果の由来が異なります。νを少し変えて連続性を見る用途でも、毎回その次数の関数を再評価し、前回の零点を無条件に次の正解として流用しません。
• 一次資料は2026年7月16日に確認しました。結果は数学的に厳密な0ではなく、監査済みbinary64近似です。関数値を再代入するときは、評価器側の丸め誤差も含めて絶対残差とbracket幅を確認してください。
• 通常Bessel Jν・Yνの正の零点について: 表示値はIEEE 754 binary64を基礎とする有限精度の近似値です。最大有効桁を表示しても全桁が数学的に保証される意味ではなく、零点付近、極付近、大きな相殺がある入力では相対誤差が増える場合があります。
• 通常Bessel Jν・Yνの正の零点について: 入力欄の上下限は単なる画面上の目安ではありません。独立した高精度oracleと境界探索で監査した公開域であり、範囲外を最後の有限値へ丸めたり、計算途中の値を成功結果として表示したりしません。
• 通常Bessel Jν・Yνの正の零点について: この計算器は外部API、外部CDN、fetch、XMLHttpRequest、WebSocket、eval、new Functionを使用しません。数値計算と結果表の組み立ては利用者のブラウザ内で完結し、共有サーバーへ数値計算の負荷を移しません。
• 通常Bessel Jν・Yνの正の零点について: 端末内計算でも、古いブラウザ、改変されたJavaScript、拡張機能の干渉、極端な省電力設定まで保証するものではありません。重要な研究、設計、試験、規格適合では別実装と一次資料を使って再計算してください。
• 通常Bessel Jν・Yνの正の零点について: 結果のコピー機能は表示値の転記を助けますが、単位、関数種、入力条件までは自動で保存しません。再現性が必要なら、入力一式、確認日、ページ名、表示された値を同じ記録へ残してください。
• 通常Bessel Jν・Yνの正の零点について: 出典URLは定義と数値上の注意を確認するための一次資料です。参照ページの説明文や実装を複製せず、入力意図だけを同定し、式、公開域、例、注意、FAQは監査結果に基づいて独自に作成します。
• 0≤ν≤10、1≤s≤10の正の実零点をbinary64近似で返します。負・複素零点、任意精度、高順位は扱いません。

よくある質問