simp3.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 20 行

void simp3(a,i1,k1,ip,kp)
float **a;
int i1,ip,k1,kp;
{
	int kk,ii;
	float piv;

	piv=1.0/a[ip+1][kp+1];
	for (ii=1;ii<=i1+1;ii++)
		if (ii-1 != ip) {
			a[ii][kp+1] *= piv;
			for (kk=1;kk<=k1+1;kk++)
				if (kk-1 != kp)
					a[ii][kk] -= a[ip+1][kk]*a[ii][kp+1];
		}
	for (kk=1;kk<=k1+1;kk++)
		if (kk-1 != kp) a[ip+1][kk] *= -piv;
	a[ip+1][kp+1]=piv;
}