double ellipsoid_area ( double a, double b, double c );
double elliptic_inc_em ( double phi, double m );
double elliptic_inc_fm ( double phi, double m );
double rd ( double x, double y, double z, double errtol, int *ierr );
double rf ( double x, double y, double z, double errtol, int *ierr );