步骤

Gauss-Seidel迭代法

步骤1

输入误差精度 $\epsilon \left(\epsilon ={10}^{-5}\right)$ ， ${\mu }_2$ ， ${\mu }_b$ ，c， $\beta$ ， $\eta$ ， $\alpha$ ， $p_1$ ， $p_2$ ， $p_3$ ，K，M， $\theta$ 和 $R=0$

步骤2

输入 $A$ ， $B$ ， $C$

步骤3

定义 $n=1$

$R_{n-1}=R$

$R_n=-\left(R_{n-1}^{2}B+C\right)A^{-1}$

步骤4

当 ${\Vert R_n-R_{n-1}\Vert }_{\infty }>\epsilon$ 时

$n=n+1$

回到步骤3

否则

回到步骤5

步骤5

$R=R_n$