18 lines
271 B
C
18 lines
271 B
C
|
|
void cholsl(a,n,p,b,x)
|
|
float **a,b[],p[],x[];
|
|
int n;
|
|
{
|
|
int i,k;
|
|
float sum;
|
|
|
|
for (i=1;i<=n;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<=n;k++) sum -= a[k][i]*x[k];
|
|
x[i]=sum/p[i];
|
|
}
|
|
}
|