#CH0402P3533. 方程求解

方程求解

Description

编写一个程序,用Jacobi迭代法求解n阶线性方程组Ax=b。输入数据确保Jacobi迭代法收敛。
Jacobi迭代法的迭代公式:

Input Format

第1行有一个正整数n和一个正的双精度数 ε,n是方程组的阶数、 ε是结果的精度要求。接下来是n行数据,每行包含n+1个双精度数。第i行的前n个数据是第i个方程的系数向量,最后一个数据是第i个方程的右式值。

Output Format

一个整数,表示所计算出的解向量使得方程组中多少个方程的左式值与右式值相差超过右式值的万分之一。

5 0.000294118
-0.329221 -0.054962 0.0968168 -0.0723648 0.0520982 0.0529794 
-0.00123477 0.161426 0.0619855 0.050214 0.0195094 -0.0284821 
0.0251155 -0.0314213 0.104726 -0.034825 0.0124324 -0.000931532 
-0.0265115 -0.00221929 0.0267514 0.132498 0.0205734 -0.0564419 
0.0711287 -0.173507 0.161147 0.122201 0.54361 -0.015625

4
因为计算出的解向量为[-0.0890094 -0.00092776 -0.146684 -0.433396 0.123449],
由此计算出的5个方程的右式值分别是0.0529473、-0.0284863、-0.000940252、-0.0564465、-0.0156611