Steenstrup
Steenstrup

Reputation: 644

EllipticE in maple compared to matematica and python

Trying to implement Euler Elastica in maple where Jacobi elliptic functions and Elliptic Integration is vital. The Jacobi elliptic functions: JacobiSN, JacobiCN, JacobiDN, and JacobiAM seems to work as I would expect, but Elliptic Integration: EllipticE behaves inacurate to what I would think it should do.

Have made a comparison between maple, matematica and python buildin function:

Maple: 
plot([t, EllipticE(t), t = 0 .. 1])

enter image description here

Matematica: 
plot[EllipticE[t], {t, 0, 1}]

enter image description here

Python(scipy.special):     
t = arange(0, 1, 0.001)
plot(t, ellipe(t))

enter image description here

The only thing that I could find was that: Maple calculate the Incomplete and complete elliptic integrals of the second kind.

Where python and matematica calculate the complete elliptic integrals of the second kind

Upvotes: 2

Views: 422

Answers (1)

Steenstrup
Steenstrup

Reputation: 644

Love and hade when I realise the answer 4 min after posting the question. maple use m where python and matematica uses k. where k^2 = m

plot([t, EllipticE(sqrt(t)), t = 0 .. 1])

enter image description here

Upvotes: 2

Related Questions