up_triangle.h

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

int up_triangle(double UU[][n],double bb[])
{
	int i,j;
	for(j=n-1;j>=1;j=j-1)
	{		
	if(fabs(UU[j][j])<eps)
		 return(-1);
		else 
		{
			bb[j]=bb[j]/UU[j][j];
		for(i=0;i<=j-1;i++)
			bb[i]=bb[i]-bb[j]*UU[i][j];		
		}
	}	
	bb[0]=bb[0]/UU[0][0];
return(0);
}