zep1s

    １次元割線法を用いて，超越方程式 \(f(x)=v\) を解く倍精度サブルーチン．
     
     
  • サブルーチン名： zep1s

    •  
  • 呼び出し方：call zep1s(y0, y1, xd, value, epsl, lmax, func, icode)

    •  
  • 引　数：
    • 引数
    意　味
    入出力の区分
    備　考
    y0 解が存在する区間の下限
    入力
    		 解が１つだけ存在する区間 (y0,y1) を予め指定しなければならない．
    y1 解が存在する区間の上限
    入力
    xd 解
    出力
    value 超越方程式の右辺の値 \(v\)
    入力
    epsl 収束判定子
    入力
    lmax 最大反復回数
    入力
    func 超越方程式の左辺の関数 \(f(x)\) を示す関数名
    入力
    icode エラー・コード
    出力
    		 解が収束した場合には，  icode =   0.
    反復が不十分な場合には，icode = 100.

  • 注　意

    • 本ルーチンを参照するプログラム単位中では，関数 func を external 文で宣言しておかねばならない．