gauss.h

利用gauss分解求解线性方程组

int Gauss(double A[][n])
{
	int i,j,k;
	double max;
	for(k=0;k<n-1;k++)
	{
		max=0;
		if(fabs(A[k][k])<eps)
		  return(-1);
		else
		{
			for(i=k+1;i<=n-1;i++)
			{
				A[i][k]=A[i][k]/A[k][k];
				for(j=k+1;j<=n-1;j++)
					A[i][j]=A[i][j]-A[i][k]*A[k][j];
			}
		}
	}
	return(0);
}