#include <iostream>  
#include <math.h>  
using namespace std; 
  
void main()  
{   
    int i=0,j=0,k=0,n=1,t=0,p=0;
    double s=0.0;
	double q=0.0;
    double a[100][100];
    double c[100];
       
    cout<<"請輸入你想處理方程的元數N：\n";
    cin>>n;  
    for(i=0;i<n;i++)
     {
       cout<<"系數"<<i+1<<"  ";
     } 
    cout<<"常數C"<<endl;
    for(i=0;i<n;i++)
     {
      cout<<"請依次輸入第"<<i+1<<"個方程的系數A和常數C："<<endl;
      for(j=0;j<n;j++)
        {
          cin>>a[i][j];
        }  
      cin>>c[i];   
     }
 
  for(i=0;i<n;i++)
  { s=a[i][i];
	for(j=0;j<n;j++)
	{ 	   
	   a[i][j]/=s;
	}	
	c[i]=c[i]/s; 
	for(p=0;p<n;p++)
	{ q=a[p][i];
      for(t=0;t<n;t++)
	  {
		if(p!=i) 		  
		  a[p][t]-=(a[i][t]*q);
	  }
      if(p!=i) 
		c[p]-=c[i]*q;
	   
	}

  }    
    cout<<"未知數為："<<endl;
    for(i=0;i<n;i++)
	{		
		cout<<"x["<<i+1<<"]="<<c[i]/a[i][i]<<endl;
	}   
}