高斯 Jordan消元法 c语言代码

//使用Gauss Jordan消去法求解线性方程组

#include “stdio.h“
#include “stdlib.h“
#include “math.h“

double **A *b *x;
unsigned int RANK = 4;
int N=4;                     //保留N位小数
unsigned int sign = 0;           //保留小数标志位

void baoliu（double &sp）
{

	int x = pow（10.0N）;
	if（sign）
	{
		if（sp>0）
		{
			sp=（ （double）（ （int）（ （sp+0.5/x）*x ） ） ）/x;
		}
		else
		{
			sp=（ （double）（ （int）（ （sp-0.5/x）*x ） ） ）/x;
		}
	}

}

void xianshi（）
{
	for （int x = 0; x < RANK; x++）
	{
		for （int y = 0; y		{
			printf（“%f\t“ A[x][y]）;
		}
		printf（“    %f“ b[x]）;
		printf（“\n“）;
	}


	printf（“\n“）;
	printf（“\n“）;
}

unsigned int makematrix（）
{
	unsigned int r c;

	printf（“请输入矩阵行列数，用空格隔开：“）;
	scanf_s（“%d %d“ &r &c）;

	A = （double**）malloc（sizeof（double*）*r）;//创建一个指针数组，把指针数组的地址赋值给a *r是乘以r的意思
	for （int i = 0; i < r; i++）
	A[i] = （double*）malloc（sizeof（double）*c）;//给第二维分配空间
	for （int i = 0; i < r; i++） {
		   for （int j = 0; j < c; j++）
			A[i