down_triangle.h

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

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