improved_cholesky.h

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

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