高斯消元法（不列主元）.txt

来自「高斯消元法（不列主元）大学计算方法里面的一个算法」· 文本 代码 · 共 43 行

//高斯消元法（不列主元）
#include <iostream.h>
#include <math.h>

int main()
{
    int  i,j,k,n;
    double matrix[50][50];
    cout<<"请输入方程阶数：";
    cin>>n;
    for (i=1;i<=n;i++)                              //输入矩阵
       {
        for (j=1;j<=n+1;j++)
           {
            cout<<"a["<<i<<"]["<<j<<"]";
            cin>>matrix[i][j];
           }
	}
    for(k=1;k<=n-1;k++)                           //迭代
	{
		if (matrix[k][k]==0) {cout<<"对角元为零，消元失败。";return 0;}
        for(i=k+1;i<=n;i++)
		{
			matrix[i][k]=matrix[i][k]/matrix[k][k];    //l[i][k]=matrix[i][k]
			for (j=k+1;j<=n+1;j++)
			{
				matrix[i][j]=matrix[i][j]-matrix[i][k]*matrix[k][j];
			}
		}
	}
 	for(k=n;k>=1;k--)                          //回代
	{
		double temp=0;
		for(i=k+1;i<=n;i++)temp=temp+matrix[k][i]*matrix[i][n+1];
		matrix[k][n+1]=(matrix[k][n+1]-temp)/matrix[k][k];
	}
	cout<<"方程组的解为："<<endl;
	for(k=1;k<=n;k++)
	{
		cout<<"a["<<k<<"]="<<matrix[k][n+1]<<"  ";
	}
	return 0;
}