scipy.special.ellipeinc

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

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