ENH: scipy.special, Implemented Ellipsoidal harmonic function: the first kind, the second kind and normalization constant 

This is an attempt to implement ellipsoidal harmonic functionality. 
The functions implemented are
ellip_harm, ellip_harm_2, ellip_normal
The sources referred to are:
http://www.geology.osu.edu/~jekeli.1/OSUReports/reports/report_499.pdf
http://arxiv.org/abs/1204.0267
for the absence of any of the sources on the list, the build fails at present,
The errors are:
undefined symbol: wrap_PyUFunc_getfperr
undefined symbol: sf_error_check_fpe (sf_error) 
undefined symbol: expitf (logit)
undefined symbol: cdft1_wrap (cdfwrappers)
undefined symbol: oblate_radial2_nocv_wrap (specfun_wrappers),
undefined symbol: cbesy_wrap_e_real (amos_wrappers.c)