elle.c

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

#include <math.h>
#include "nrutil.h"

float elle(phi,ak)
float ak,phi;
{
	float rd(),rf();
	float cc,q,s;

	s=sin(phi);
	cc=SQR(cos(phi));
	q=(1.0-s*ak)*(1.0+s*ak);
	return s*(rf(cc,q,1.0)-(SQR(s*ak))*rd(cc,q,1.0)/3.0);
}
/* (C) Copr. 1986-92 Numerical Recipes Software . */