雅各比和高斯塞代尔迭代法.cpp

数值分析的实验报告

//雅各比和高斯塞代尔迭代法
#include<iostream.h>
#include <stdio.h>
#include <math.h>
#define NP 5
#define ERR 1e-15
int k;
double X[NP];
void jacobi(double A[NP][NP],double f[NP])
{ 
	  int i,j;
	  double sum=0,temp=0,X0,err_max,err_x;
	  do
	  {
		 err_max=0;
		 for(i=0; i<NP; i++)
		 { 
			 sum=0.0;
			 for(j=0; j<NP;j++)
			 {   
				 if(j!=i) 
				 sum=sum+A[i][j]*X[j];
			 }
			 X0=X[i];
			 X[i]=(f[i]-sum)/A[i][ i];
			err_x=fabs(X[i]-X0);
			 if (err_x>err_max)
				 err_max=err_x;
		 }
	  }
	  while (err_max>ERR);
}
int main()
{ 
	int i,numb,k,n;
	cout<<"请输入矩阵的行列个数"<<endl;
	cin>>numb;
	static double a1[NP][NP];
	for(k=1;k<numb+1;k++)
	{
		for(n=1;n<numb+1;n++)
		{
			if(k==n)
				a1[k][n]=2;
			else if (k==(n-1))
				a1[k][n]=1;
			else if (k==(n+1))
				a1[k][n]=1;
			else
				a1[k][n]=0;
		}
	}
	static double f[NP]={-7,-5,0,-5};
	jacobi(a1,f);
	for (i=0;i<NP;i++)
	   cout<<X[i]<<endl;
	return 0;
}