Symmetrised elliptic integral of 3rd kind RJ(x,y,z,r)
double elliptic_integral_rj( double x, double y, double z, double r )
Returns a value of the symmetrised elliptic integral of 3rd kind
EX1
void elliptic_integral_rj_ex1() { double dx[6] = {0.5,0.5,0.5,1.0,1.0,1.5}; double dy[6] = {0.5,1.0,1.5,0.5,1.0,1.5}; double dz[6] = {0.5,0.5,1.0,1.0,1.5,1.5}; double dr = 2.0; double dResults[6]; for(int ii = 0; ii<6 ; ii++) { dResults[ii] = elliptic_integral_rj(dx[ii],dy[ii],dz[ii],dr); printf("%lf %lf %lf %lf: %lf\n",dx[ii],dy[ii],dz[ii],dr,dResults[ii]); } }
elliptic_integral_rc, elliptic_integral_rf, elliptic_integral_rd
origin.h