zep1n

    1次元 Newton 法を用いて,超越方程式 f(x) = v を解く倍精度サブルーチン.
     
     
  • サブルーチン名: zep1n

  •  
  • 呼び出し方:call zep1n(x0,value,epsl,det0,lmax,func,dfun,icode)

  •  
  • 引 数:
  • 引数
    意 味
    入出力の区分
    備 考
    x0 入力時は,近似解の初期推定値
    出力時は,近似解
    入力&
    出力 
    value 超越方程式の右辺の値 v
    入力
    epsl 収束判定子
    入力
    det0 特異性判定子
    入力
    | f'(x0) | ≦ det0 のとき,Newton法を実行しない.
    lmax 最大反復回数
    入力
    func 超越方程式の左辺の関数 f(x) を示す関数名
    入力
    2≦n≦ndim を満足するように選ぶ.
    dfun 導関数 f(x) を示す関数名
    入力
    icode エラー・コード
    出力
    解が収束した場合には,  icode =   0.
    反復が不十分な場合には,icode = 100.
  • 注 意

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