Subst: y=b²/a²
∫ sqrt(1+y*z²/(1-z²))dz = (sqrt(1-z^2)*sqrt((z^2-y*z^2-1)/(z^2-1))*EllipticE2(asin(z), 1-y))/sqrt(1-z^2+y*z^2)
EllipticE(asin(z), 1-y),z=1 ergibt EllipticE(1-y) und bei z=0 ergibt das 0
sqrt(1-z^2)*sqrt((z^2-y*z^2-1)/(z^2-1))/sqrt(1-z^2+y*z^2),z=1 ergibt +/-1 bei z=0 ergibt das 1
Differenz=EllipticE(1-y)
also Endergebnis
L=a*EllipticE(1-b²/a²)
Weg2 (was laut Dir einfacher sein soll):
∫ sqrt(1-y*cos(x)²)dx = ((1-y) sqrt((-2+y (1+cos(2 x)))/(-1+y)) EllipticE2(x, y/(-1+y)))/sqrt(2-y (1+cos(2 x)))
ergibt auch was mit EllipticE2... was soll da einfacher sein?
Differenz: sqrt(1/(1-y)) (1-y)*EllipticE(y/(-1+y)) - 0
also Endergebnis mit Rücksubst:
L=a*b/a*EllipticE(1-a²/b²)
Probe mit numerischer Integration und a=7, b=5
L(7,5)=7*1.3557631762415649910305354634865588
=9.490342233690954937213748244405...
symbolisch:
7*EllipticE(1-5²/7²)=7*EllipticE(24/49)
z.B. http://www.lamprechts.de/gerd/php/RechnerMitUmkehrfunktion.php
=7*1.3557631762415649910305354634865588
=9.490342233690954937213748244405... stimmt
Probe 2. Integral erst numerisch:
integrate sqrt(1-(1-5²/7²)cos(x)²)dx,x=0...Pi/2
=7*1.3557631762415649910305354634865588
=9.490342233690954937213748244405...
symbolisch:
7*5/7*EllipticE(1-7²/5²)=5*EllipticE(-24/25)
=5*1.898068446738190987442749648881182...
=9.490342233690954937213748244405...
Alle 4 stimmen, weil
EllipticE[1 - 1/z] = EllipticE[1 - z]/Sqrt[z] für Abs[Arg[z]] < Pi