cholsl.c

统计模式识别算法包

/* This was taken from the new version of the book */

#include <math.h>

void cholsl(a,n,p,b,x)
float **a,*p,*b,*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];
	}
}