guass_c_resolution.h

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

void Guass_C_resolution(double A[][n],double bb[n])
 {
	int i,j,p;
	int u[n];
	double LL[n][n],UU[n][n],temp;
	int Gauss_C(double A[][n],int u[n]);
	for(i=0;i<n;i++)
	  u[i]=i;
	 Gauss_C(A,u);
	for(i=0;i<n;i++)
	  LL[i][i]=1;
	for(i=0;i<n;i++)
	{		
		for(j=0;j<i;j++)
			LL[i][j]=A[i][j];
		for(j=i;j<n;j++)
			UU[i][j]=A[i][j];
	}
    for(i=0;i<n-1;i++)
	{
	p=u[i];
	if(i!=p)
	{
	  temp=bb[i];
	  bb[i]=bb[p];
	  bb[p]=temp;
	}
	}
	down_triangle(LL,bb);    
	up_triangle(UU,bb);
}