XLPack 7.0
XLPack 数値計算ライブラリ (C API) リファレンスマニュアル
読み取り中…
検索中…
一致する文字列を見つけられません

◆ cpzero()

void cpzero ( int  n,
doublecomplex  a[],
doublecomplex  r[],
int  iflag,
int  maxiter,
int *  iter,
double  s[],
doublecomplex  work[],
int *  info 
)

高次代数方程式 (複素係数, 複素解) (連立法)

目的
cpzeroは複素係数の多項式p(z)のすべてのゼロ点を連立法により求める.
p(z) = a0*z^n + a1*z^(n-1) + ... + an
引数
[in]n多項式の次数. (n >= 1)
[in]a[]配列 a[la] (la >= n + 1)
p(z)の複素係数ベクトル (a0 〜 an).
[in,out]r[]配列 r[lr] (lr >= n)
[in] 方程式の解の初期推定値. 不明であれば, iflag = 0 として推定値を設定しなくてもよい.
  注 - 初期推定値はすべて異なった値であること.
[out] 求められた解.
[in]iflag方程式の解の初期推定値の入力フラグ.
= 0: 初期推定値は与えられていない.
!= 0: r[]に初期推定値が与えられている.
[in]maxiter最大反復回数. (maxiter >= 1)
[out]iter収束に要した反復回数.
[out]s[]配列 s[ls] (ls >= n)
求められた解の誤差限界.
[out]work[]配列 work[lwork] (lwork >= 2*(n + 1))
作業領域.
[out]info= 0: 正常終了
= -1: 入力パラメータ n の誤り (n < 1)
= -2: 入力パラメータ a の誤り (a[0] = 0)
= -5: 入力パラメータ maxiter の誤り (maxiter <= 0)
= 1: maxiter回の反復で収束しなかった。ゼロ点の最終推定値がr[]に入る. 誤差限界s[]は計算されない.
出典
SLATEC