高斯选主元消去法.cpp

计算方法中

#include<iostream.h>
#include<iomanip.h>
#define N 20

void main( )
{
	int n,i,j,k;
	cout<<"输入方程的阶数:"<<endl;
	cin>>n;
	float A[N][N],X[N],b[N],a,l;
	for(i=1;i<=n;i++)
	{
		cout<<"输入第"<<i<<"个方程的系数和b的值:"<<endl;
		for(j=1;j<=n;j++)cin>>A[i][j];
		cin>>b[i];
	}
	for(k=1;k<=n;k++)
	{
		if(A[k][k]==0)
		{
			cout<<"error!";
			break;
		}
		for(i=k+1;i<=n;i++)
		{
			l=A[i][k]/A[k][k];
			for(j=k+1;j<=n;j++)
		        A[i][j]=A[i][j]-l*A[k][j];
			b[i]=b[i]-l*b[k];
		}
	}
    X[n]=b[n]/A[n][n];
	for(i=n-1;i;i--)
	{
		a=0;
		for(j=i+1;j<=n;j++)
			  a+=A[i][j]*X[j];
		X[i]=(b[i]-a)/A[i][i];
	}
	cout<<endl;
	for(i=1;i<=n;i++)
	   cout<<"X["<<i<<"]="<<X[i]<<endl;
}