cholesky法求解方程组

for(i=2;i<(n+1);i++) { for(j=1;j<i;j++) { sigma1=0.0; for(k=1;k<j;k++) sigma1+=t[i-1][k-1]*L[j-1][k-1]; t[i-1][j-1]=a[i-1][j-1]-sigma1; L[i-1][j-1]=t[i-1][j-1]/d[j-1]; } sigma2=0.0; for(k=1;k<i;k++) sigma2+=t[i-1][k-1]*L[i-1][k-1]; d[i-1]=a[i-1][i-1]-sigma2; } y[0]=b[0]; for(i=2;i<(n+1);i++) { sigma3=0.0; for(k=1;k<i;k++) sigma3+=L[i-1][k-1]*y[k-1]; y[i-1]=b[i-1]-sigma3; } x[n-1]=y[n-1]/d[n-1]; for(i=(n-1);i>0;i--) { sigma4=0.0; for(k=(i+1);k<(n+1);k++) sigma4+=L[k-1][i-1]*x[k-1]; x[i-1]=y[i-1]/d[i-1]-sigma4; } for(i=0;i<n;i++) *(xx+i)=x[i];