xsphfpt.c

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

/* Auxiliary routine for sphfpt */

extern int m,n;
extern float c2,dx,gmma;

void derivs(float x,float y[],float dydx[])
{
	dydx[1]=y[2];
	dydx[2]=(2.0*x*(m+1.0)*y[2]-(y[3]-c2*x*x)*y[1])/(1.0-x*x);
	dydx[3]=0.0;
}
/* (C) Copr. 1986-92 Numerical Recipes Software . */