SciPy

scipy.special.pro_rad1

scipy.special.pro_rad1(x1, x2, x3, x4[, out1, out2]) = <ufunc 'pro_rad1'>

(s,sp)=pro_rad1(m,n,c,x) computes the prolate sheroidal radial function of the first kind and its derivative (with respect to x) for mode parameters m>=0 and n>=m, spheroidal parameter c and |x| < 1.0.