第2種Bessel関数Y_ν(x)を実数次数でグラフ・表計算

第2種Bessel関数Y_ν(x)を実数次数でグラフ・表計算

実数次数νと正のx格子から第2種Bessel関数Y_ν(x)を計算し、原点近くの次数依存挙動と振動を確認します。

入力

−20以上20以下の実数です。整数次数だけでなく0.5なども指定できます。

0.0001以上50以下の第0点です。

0以外の有限値です。負なら大きいxから小さいxへ並べます。

初期値を含む1〜201点です。

使い方

  1. 実数次数νと正のx格子を入力します。既定はν=0.5、x=1〜3です。
  2. Y_(1/2)(x)=−sqrt(2/(πx))cos xを表の代表点と照合します。
  3. 原点近くの急増を調べるときはxの下限0.0001を守り、対数軸ではないことに注意します。

計算式・考え方

第2種Bessel関数Y_ν(x)はJ_νと独立なx²y''+xy'+(x²−ν²)y=0の標準解です。非整数νではY_ν(x)={cos(πν)J_ν(x)−J_(−ν)(x)}/sin(πν)で定義でき、整数次数はその極限です。本ページは実数次数νへ拡張し、次数で異なる原点極限と直接積分の共通安全域を踏まえて0.0001≤x≤50だけを評価します。

注意事項

  • NIST DLMF §10.2、§10.9および実数積分表示10.9.E6・E7を2026年7月15日に確認しました。Neumann関数N_νと表記する資料もありますが、本ページではDLMFのY_νへ統一します。
  • Y_νのx→0+の極限は次数に依存し、負の半整数などでは有限値や0へ近づく例もあります。本ツールは全次数に共通する直接積分の安全域として入力下限を0.0001に固定します。xが小さくν>0なら主要部は−Γ(ν)(2/x)^ν/πとなり、グラフの縦軸が非常に大きな負値へ引かれる場合があります。これは次数に応じた原点近傍の挙動で、描画エラーとは限りません。
  • νが整数に近いと接続公式の分母sin(πν)が小さくなります。実装は整数近傍をそのまま不安定な差で計算せず、級数・漸近式・対称関係を切り替える有界evaluatorを使います。それでも厳しい桁が必要なら任意精度計算で確認してください。
  • J_νとY_νのWronskianはJ_νY_ν'−J_ν'Y_ν=2/(πx)です。導関数は本ページの出力に含めないため、この恒等式を使う検証では別の導関数実装と同じ次数・xを指定してください。
  • このwaveでは、次数νを整数に丸めず−20から20までの実数として扱います。古典的な球ベッセル関数は非負整数nで紹介されることが多い一方、本ページ群の球関数は円柱ベッセル関数との関係式を定義として実数次数νへ拡張しています。したがって、整数次数だけを載せた数表と比較するときも、入力欄のνが階数nなのか円柱側のν+1/2なのかを確認してください。
  • 横軸はx_i=xの初期値+i×増分(i=0,1,…,点数−1)で作ります。点数は1〜201、xは関数ごとの許可範囲内です。増分を負にすれば右から左へ並べられますが、初期点だけでなく最終点も範囲内でなければ計算しません。増分0は同じ点を重複させるだけなので入力エラーにします。
  • グラフの線は指定した標本点を見やすく結んだものです。点間を追加評価した補間曲線、零点探索、極値探索ではありません。振動が速い区間やx=0近くの急変を調べる場合は表示範囲を狭め、増分を小さくしてください。見かけ上の交点をそのまま根とせず、表の符号と専用の高精度計算で再確認してください。
  • 同じBesselという名称でも、第1種と第2種、通常型と変形型、円柱型と球型では微分方程式、原点での挙動、無限遠での増減が異なります。さらに資料によって球型の正規化やmodified spherical kの係数が違います。比較時は関数記号だけでなく、次数、引数、平方根係数、位相、正規化を式単位で照合してください。
  • 計算はschema v4のsampled-series-v1と、固定名で許可した有界なBessel evaluatorを使います。任意の式やコードを実行するeval、new Function、外部API、CDN、fetch、XMLHttpRequest、WebSocketは使いません。反復回数、入力範囲、有限性を検査し、倍精度で安全に返せない標本は未定義セルとして線を分断します。
  • 入力した次数、x格子、生成した系列は利用者のブラウザ内だけで計算・表示され、サーバーや第三者へ送信・保存されません。WordPressの本文、追加CSS、追加JavaScriptで完結できる構成を前提とし、このbuilderはローカル成果物だけを生成します。下書き作成、記事更新、本公開は実行しません。
  • 結果の表とTSVはグラフと同じ内部数値を使います。表示は最大15桁の有効数字で、binary64の丸め、漸近展開や級数の打切り、近接する大きな項の相殺による誤差を含みます。極端に小さい値や大きい値では相対誤差と絶対誤差のどちらが用途に重要かを決め、NIST DLMFの表や任意精度ライブラリとも照合してください。
  • NIST Digital Library of Mathematical Functions(DLMF)第10章の定義、接続公式、実数引数の積分表示、球関数との関係を2026年7月15日に一次資料として確認しました。ページごとの本文に参照節を記載しています。本ツールは公式定義の学習と代表値の検算を補助しますが、解析接続、複素引数、導関数、零点、証明、任意精度計算を提供するものではありません。
  • Bessel関数は波動方程式、円筒・球対称の境界値問題、熱伝導、電磁気、拡散モデルなどに現れます。ただし、このグラフだけから物理量を決めることはできません。モデル側の境界条件、単位、無次元化、線形結合係数を確認し、採用した関数規約と入力値を結果と一緒に記録してください。
  • 有界な倍精度による実数次数の標本グラフです。任意精度、複素数、零点探索、物理モデルの妥当性を保証しません。

よくある質問

Neumann関数と同じですか?

一般にN_νとも書かれる第2種Bessel関数Y_νです。

なぜx=0を入力できませんか?

原点極限は次数で異なりますが、本ツールでは全次数に共通する直接積分の安全域としてx=0を一律に除外しています。

半整数次数の検算式は?

Y_(1/2)(x)=−sqrt(2/(πx))cos xです。

整数に近いνも使えますか?

使えますが、高精度用途では整数近傍を任意精度実装とも比較してください。

入力した実数次数とx格子、計算結果は端末内だけで処理され、外部へ送信または保存されません。個人情報の入力欄もありません。