//使用Gauss?Jordan消去法求解線性方程組

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

double?**A?*b?*x;
unsigned?int?RANK?=?4;
int?N=4;?????????????????????//保留N位小數(shù)
unsigned?int?sign?=?0;???????????//保留小數(shù)標(biāo)志位

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?	{
		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（“請輸入矩陣行列數(shù)，用空格隔開：“）;
	scanf_s（“%d?%d“?&r?&c）;

	A?=?（double**）malloc（sizeof（double*）*r）;//創(chuàng)建一個指針數(shù)組，把指針數(shù)組的地址賦值給a?*r是乘以r的意思
	for?（int?i?=?0;?i?	A[i]?=?（double*）malloc（sizeof（double）*c）;//給第二維分配空間
	for?（int?i?=?0;?i?		???for?（int?j?=?0;?j?			A[i