scipy.special.ellipkinc

scipy.special.ellipkinc(x1, x2[, out])

y=ellipkinc(phi,m) returns the incomplete elliptic integral of the first kind: integral(1/sqrt(1-m*sin(t)**2),t=0..phi)