equation_selfconsistent プログラム仕様

自己無撞着法（不動点反復法）を用いて方程式の解を求めるスクリプト。

このモジュールは、x = f(x) の形式で表現された方程式を、反復式 x_{n+1} = f(x_n) を用いて解きます。 初期値から開始し、指定された収束条件または最大反復回数に達するまで計算を繰り返します。 混合比を用いて反復を安定化させることも可能です。

equation-selfconsistent_usage

cms.equation.equation_selfconsistent.func(x)

自己無撞着法で解く方程式の右辺を定義します。

この関数は、x = f(x) 形式の方程式における f(x) に相当します。 具体的には、f(x) = 0.25 * (-x^3 + x^2 - 2.0) を計算します。

パラメータ:

x -- float 独立変数。

戻り値:

float 計算された f(x) の値。

cms.equation.equation_selfconsistent.main()

自己無撞着法を実行し、方程式の解を求めます。

グローバル変数 x0 (初期値)、kmix (混合比)、eps (収束判定条件)、 nmaxiter (最大反復回数) を使用して、反復計算を行います。 計算の過程と結果は標準出力に表示されます。

戻り値:

int 収束した場合は1、収束しなかった場合は0。