c++线性方程组求解源代码 原创

#include 
#define M 3
#define N （2*M）
int main（）
{
	int ijk;
	double a[M][M]={123221343};;
	double result1[M]={-214};
	double result[M][M];
	double x[M]={0};
	double b[M][N];
	cout<<“请输入矩阵的值（默认大小为3*3的矩阵）：“<	for（i=0;i	{
		for（j=0;j		{
			cin>>a[i][j];
			b[i][j]=a[i][j];
		}
	}
/*****************扩展矩阵***********************/
	for（i=0;i	{
		for（j=M;j		{
			if（i==（j-M））
			{
				b[i][j]=1;
			}
			else
			{
				b[i][j]=0;
			}
		}
	}
/*****************扩展矩阵***********************/

/*****************求逆模块***********************/
	for（i=0;i	{
		if（b[i][i]==0）
		{
			for（k=i;k			{
				if（b[k][k]!=0）
				{
					for（int j=0;j					{
						double temp;