Jacobi迭代和Gauss迭代 c语言实现

/*
*TestMain.cpp
*功能：利用Jacobi迭代和Gauss迭代求解方程组
*时间：2013.4.1
*/
#include "Jacobi.h"
#include "Gauss.h"
int main ()
{
	printf("Jacobi迭代的结果：\n");
	Jacobi(3 , 0.0001);
	printf("*****************************************************\n");
	printf("Gauss迭代的结果：\n");
    Gauss(3,0.0001);
	return 0;
}
/*
*Gauss.cpp
*功能：高斯迭代
*时间：2013.4.1
*/
#include "Gauss.h"

void Gauss(int n , double E)
{
	int i; // 用来循环
	int k=0; //记录迭代次数
	double maxE;

	//设置初始值
	double *a = (double *)malloc(sizeof(double)*(n+1));
	for(i=1;i<=n;i++)
	{
		a[i]=0;
		printf("%d :X0=%f ,",k,a[i]);
	}

	printf("\n");
	//存放新的解
	double *temp = (double *)malloc(sizeof(double)*(n+1));
    do
	{
		CalculateG(a,temp);
		k++;
		printf("%d :",k);
		for(i=1;i<=n;i++)
		{
			printf("X%d=%f ,",i,temp[i]);
		}
		printf