对阵正定矩阵的楚列斯基（Cholesky）分解（C语言）

//定理2.2.3：对阵正定矩阵的楚列斯基（Cholesky）分解
//设A为n阶对阵正定矩阵，则存在一个可逆的下三角矩阵G，使得
//A=GG’，当限定G的对角元为正时，这种分解是唯一的
//
//--------A=GG’的分解算法-------
//参考教材：《数值分析》李乃成，梅立泉，科学出版社
//    《计算方法教程》第二版 凌永祥，陈明逵
#include<stdio.h>
#include<math.h>

int min(int a,int b);
int main(void)
{
 int size=10;
 int i,j,p=0;
 double A[10][10] = {0.0};
 double sum,up=0.0;
 printf("矩阵A[10][10]:\n");
 for(i=0;i<size;i++)
 {
  for(j=0;j<size;j++)
  {
   if(i==j)
    A[i][i]=i+1;
   else
    A[i][j]=min(i,j)-1;
   printf("%5.4f ",A[i][j]);
  }
  printf("\n");
 }
 A[0][0]=sqrt(A[0][0]);
 for(i=1;i<size;i++)
 {