After some more digging I found out that these two functions take slightly different inputs. ellipticE() assumes the input is square, so the following will be zero:
mfun('EllipticE',0.5)-ellipticE(0.5^2)
Case closed!
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!