zep1b

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

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

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

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