cholesky_resolution.h

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

void Cholesky_resolution(double A[][n],double bb[n])
{
	int i,j;
	double LL[n][n],UU[n][n];
	int Cholesky(double A[][n]);
	Cholesky(A);
	for(i=0;i<n;i++)
	{
		for(j=0;j<=i;j++)
		{
			LL[i][j]=A[i][j];
			UU[j][i]=LL[i][j];
		}
	}
	down_triangle(LL,bb);
	up_triangle(UU,bb);
}