cholesky.h

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

int Cholesky(double A[][n])
{
	int i,j,k;
	for(k=0;k<n;k++)
	{
		A[k][k]=sqrt(A[k][k]);
		for(i=k+1;i<n;i++)
		{
			A[i][k]=A[i][k]/A[k][k];
		}
		for(j=k+1;j<n;j++)
		{
			for(i=j;i<n;i++)
			{
				A[i][j]=A[i][j]-A[i][k]*A[j][k];
			}
		}
	}
	return(0);
}