asolve.c

来自「Numerical Recipes Software 提供的算法子程序集」· C语言 代码 · 共 14 行

extern unsigned long ija[];
extern double sa[];

void asolve(n,b,x,itrnsp)
double b[],x[];
int itrnsp;
unsigned long n;
{
	unsigned long i;

	for(i=1;i<=n;i++) x[i]=(sa[i] != 0.0 ? b[i]/sa[i] : b[i]);
}
/* (C) Copr. 1986-92 Numerical Recipes Software . */