ENH: scipy.special, Implemented Ellipsoidal harmonic function, the first kind

Here is an attempt to implement Ellipsoidal Harmonic Function (Lame's Function) the first kind.
The following DE is taken as the Lame's Equation
(s^2 - h^2)*(s^2 - k^2) d2E(s)/ds2 + s*(2*s^2 - h^2 - k^2) dE(s)/ds + (p - q*s) E(s) = 0
Other assumptions and references have been included in the Docs.















 
Is there a strong reason not to use E''(s) for the second derivative etc