plot (2 (EllipticE[-a^2] + Sqrt[1 + a^2] EllipticE[a^2/(1 + a^2)]))/a for a=1 to 100