相关代码 |
|
void cholsl(a,n,p,b,x)
float **a,b[],p[],x[];
int n;
{
int i,k;
float sum;
for (i=1;i for (sum=b[i],k=i-1;k>=1;k--) sum -= a[i][k]*x[k];
x[i]=sum/p[i];
}
for (i=n;i>=1;i--) {
for (sum=x[i],k=i+1;k x[i]=sum/p[i];
}
}
/* (C) Copr. 1986-92 Numerical Recipes Software . */
相关资源 |
|