dawlg

    \(W(x)=\displaystyle \log \frac{1}{x}\) を重み関数とし，\([0,1]\)を積分区間とする定積分 \(I=\displaystyle \int_0^1 f(\xi)\log \frac{1}{\xi}\,d\xi\) に対する Gauss 積分の積分点 \(\xi_i\) と重み \(w_i\) \((i=1,2,\cdots,n)\)を自動的に計算する倍精度サブルーチン．
     
  • サブルーチン名： dawlg
  • 呼び出し方：call dawlg(nval, lmax, epsln, gzi, weit)
  • 引　数( dawlgルーチン)：
    • 引数
    意　味
    入出力の区分
    備　考
    nval 積分点の数 入力 2*nval-1 が積分公式の次数になる
    lmax ２分法の最大反復回数 入力 \(1000\) 程度で十分である
    epsln ２分法の収束判定子 入力 \(10^{-14}\) 程度
    gzi(*) 積分点 \(\xi_i\) \((i=1,2,\cdots,n)\)を格納する配列の配列名 出力
    weit(*) 重み \(w_i\) \((i=1,2,\cdots,n)\)を格納する配列の配列名 出力