The documentation for elliptic123 leads me to believe that this call should be in the range that Passes. [K,E,Pi]=elliptic123(0.358, 0.99958) However, when I compare it with the built-in ellipticPi(0.9958, 0.358), I get very different results.