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
I was able implement the test only on the ufuncs function (ellip_harm) and not on ellip_normal. I get the following error,   
File "_ellip_harm_2.pyx", line 126, in scipy.special._ellip_harm_2._ellipsoid_norm (scipy/special/_ellip_harm_2.c:2416)
TypeError: only length-1 arrays can be converted to Python scalars
How can we go around this?