inputandoutput.h

Jacobi GaussSeidel法求算式的根

#define	MAX_n	100
#define PRECISION	0.005
#define MAX_Number	1000
void Input(float x[],int n)
{
	int i;
	for(i=1;i<=n;++i)
	{
		printf("x[%d]=",i);
		scanf("%f",&x[i]);
	}
	//return x;
}
void Matrix(float A[][MAX_n],int m,int n)
{
	//int i,j;
	//printf("\n===请输入矩阵元素===\n");
	//for(i=1;i<=m;++i)
	//{
	//	printf("请输入第 %d 行: ",i);
	//	for(j=1;j<=n;++j)
	//		scanf("%f",&A[i][j]);
	//}
	A[1][1]=-1;
	A[1][2]=0;
	A[1][3]=.1;
	A[1][4]=-0.9;
	A[2][1]=0.2;
	A[2][2]=-1;
	A[2][3]=.1;
	A[2][4]=-0.7;
	A[3][1]=0;
	A[3][2]=.2;
	A[3][3]=-1;
	A[3][4]=-0.8;
	                  
	

	//return A;
}
void Output(float x[],int n)
{
	int i;
	for(i=1;i<=n;++i)
		printf("\nx[%d]=%f",i,x[i]);
    printf("\n\n");
}
int Is(float x1[],float x2[],int n)
{
	int i;
	for(i=1;i<=n;++i)
		if(fabs(x1[i]-x2[i])>PRECISION) return 1;
	return 0;
}