Jacobi点迭代.docVIP

Jacobi点迭代法 该线性方程组的增广矩阵为： #include #include #include #include #define n 3 #define precision 1e-16 static double aa[n][n+1] 2,2,5,5 , 3,4,7,6 , 1,3,3,5 ; void main int i,j,det;double a[n+1][n+2],x[n+1]; int GaussElimination_ColumnSelect ; for i 1;i n;i++ for j 1;j n+1;j++ 【消元过程】 a[i][j] aa[i-1][j-1]; det GaussElimination_ColumnSelect a,x ; if det! 0 for i 1;i n;i++ printf \nx[%d] %f\n,i,x[i] ;printf \n ; getch ; int GaussElimination_ColumnSelect double a[][n+2],double x[n+1] int i,j,k,r;double c; for k 1;k n-1;k++ r k; for i k;i n;i++ if fabs a[i][k] fabs a[r][k] r i; if fabs a[r][k] precision printf \n det A 0.Elimination faild! ;exit 0 ; if r! k for j k;j n+1;j++ c a[k][j];a[k][j] a[r][j];a[r][j] c; for i k+1;i n;i++ c a[i][k]/a[k][k]; for j k+1;j n+1;j++ a[i][j] a[i][j]-c*a[k][j]; if fabs a[n][n] precision printf \n det A 0.Algorithm faild! ;exit 0 ; for k n;k 1;k-- 【回代过程】 x[k] a[k][n+1]; for j k+1;j n;j++ x[k] x[k]-a[k][j]*x[j]; x[k] x[k]/a[k][k]; return 1 ;