高斯消元模板

link

#include
using namespace std;
const int N=110;
double a[N][N];
int n;
void gauss()
{
	for(int i=1;i<=n;i++){
		int id=i;
		for(int j=i+1;j<=n;j++)
			if(fabs(a[j][i])>fabs(a[id][i]))
				id=j;
		if(fabs(a[id][i])<1e-6){
			puts("No Solution");
			return ;
		}
		for(int j=1;j<=n+1;j++)
			swap(a[i][j],a[id][j]);
		for(int j=1;j<=n;j++){
			if(j!=i)
			{
				double temp=a[j][i]/a[i][i];
				for(int k=i+1;k<=n+1;k++)
					a[j][k]-=temp*a[i][k];
			}
		}
	}
	for(int i=1;i<=n;i++)
	{
		printf("%.2lfn",a[i][n+1]/a[i][i]);
	}
}
int main(){
	cin>>n;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n+1;j++)
			cin>>a[i][j];
		gauss();
		
}