高斯消去法

# include<stdio.h>
# include<math.h>
# define delta 1e-6
#define N 100
int main()
{
	int i,j,t,r,n,u,c=0;
	float p,L,max,s; 
	float X[N];
	float a[N][N+1];
	printf("请输入系数矩阵的阶数：\n"); 
	scanf("%d",&n);
	printf("输入的增广矩阵系数,中间用空格隔开：\n"); 
	for(i=0;i<n;i++)
		for(j=0;j<n+1;j++)
       	scanf("%f",&a[i][j]);		
    printf("你输入的增广矩阵为：\n");
	for(i=0;i<n;i++)
	{  
		for(j=0;j<n+1;j++)
			printf("%f\t",a[i][j]); 
		printf("\n");
	}  
	for(j=0;j<n-1;j++)  
	{    
		{      
			max=fabs(a[j][j]); 
	        r=j; 
		}    
		for(i=j+1;i<n;i++)
		if(fabs(a[i][j])>max)  
		{
			max=a[i][j]; 
		    r=i;    
		}  
		if(fabs(a[i][j])<delta)
			printf("矩阵奇异");
		for(t=j;t<n+1;t++)
		{
			 p=a[j][t]; 
			 a[j][t]=a[r][t];
			 a[r][t]=p; 
		}        
		for(i=j+1;i<N;i++)  
		{       
			L=a[i][j]/a[j][j]; 
			for(t=j;t<n+1;t++)   
				a[i][t]=a[i][t]-L*a[j][t]; 
		}   
	}
	printf("输出原方程的解为：\n");
	X[n-1]=a[n-1][n]/a[n-1][n-1]; 
	for(i=n-2;i>=0;i--) 
	{
	    s=a[i][n];
	    for(j=i+1;j<n;j++)
	    s=s-a[i][j]*X[j];
	    X[i]=s/a[i][i]; 
	} 
	for(u=0;u<n;u++)
		for(j=0;j<n+1;j++)
		{
			printf("%12f",a[u][j]);
			c++;
			if(c%(n+1)==0)
			printf("\n");
		}
	for(i=0;i<n;i++)
    {  
		printf("x(%d)=%.4f\n",i+1,X[i]);
		if(i==n-1)
		printf("\n"); 
	}
	system("pause");
	return 0;
}