用C語言實(shí)現(xiàn)雅可比迭代法

1、用 C 語言實(shí)現(xiàn)雅可比迭代法用 C 語言實(shí)現(xiàn)雅可比迭代法， 源程序代碼如下 : #include "stdlib.h"#include "stdio.h"#include "conio.h"#include "string.h"#include "math.h"#define N 100float Table(int n,float aNN,float bN) int i,j;float cNN;printf("Please input the matrix A by row!n&qu

2、ot;); for(i=0;i<n;i+) printf("Row %d:",i+1);for(j=0;j<n;j+) scanf("%f",&aij);printf("Please input the vector b:"); for(i=0;i<n;i+) scanf("%f",&bi);for(i=0;i<n;i+) for(j=0;j<n;j+)if(i=j)cij=0; continue;cij=-aij/aii;printf("nThe matri

3、x A and vector b:n"); for(i=0;i<n;i+)for(j=0;j<n;j+) printf("%10.5f",aij);printf("%10.5f",bi); printf("n");printf("nThe Jacobi iterative scheme:n"); for(i=0;i<n;i+)for(j=0;j<n;j+) printf("%10.5f",cij);printf("%10.5f",bi/aii

4、); printf("n");float init_vec(int n,float xN)int i;printf("Please input the initial iteration vector x:"); for(i=0;i<n;i+)scanf("%f",&xi);printf("nThe initial iteration vector x:n"); for(i=0;i<n;i+)printf("%10.5f",xi);printf("n");

5、float jacobi(int n,float aNN,float bN,float xN)int i,j,k;float tmp,x2N;for(k=0;k+)for(i=0;i<n;i+)x2i=xi;for(i=0;i<n;i+)tmp=0.0; for(j=0;j<n;j+) if(j=i) continue; tmp+=aij*x2j; xi=(bi-tmp)/aii; for(i=0,j=0;i<n;i+) if(fabs(x2i-xi)<0.00001) j+;if(j=n) printf("nThis Jacobi iterative

6、scheme is convergent!n"); printf("Number of iterations: %d",k+1);break;if(k=499)printf("nThis Jacobi iterative scheme may be not convergent!"); break; printf("nThe results:n");for(i=0;i<n;i+) printf("%12.7f",xi); int main()int n;float xN,aNN,bN; printf("Input