xsphfpt.cpp

来自「C++数值算法原书代码」· C++ 代码 · 共 14 行

// Auxiliary routine for sphfpt

#include "nrtypes.h"

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

void derivs(DP x,DP y[],DP dydx[])
{
        dydx[0]=y[1];
        dydx[1]=(2.0*x*(m+1.0)*y[1]-(y[2]-c2*x*x)*y[0])/(1.0-x*x);
        dydx[2]=0.0;
}