2015-02-04 1 views
1

Jacobi 타원 함수와 Elliptic Integration이 필수적인 단풍 나무에 Euler Elastica를 구현하려고합니다. Jacobi 타원 함수 : JacobiSN, JacobiCN, JacobiDN 및 JacobiAM이 예상대로 작동하는 것처럼 보입니다. 그러나 Elliptic Integration : EllipticE는 내가해야한다고 생각하는 것만 큼 부적절한 행동을합니다.maple의 EllipticE와 matematica 및 python

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

:

유무 기능 buildin 메이플 matematica 파이썬를 비교했다

찾을 수있는 유일한 것은 : Maple은 불완전하고 두 번째 종류의 타원 적분을 계산합니다.

어디 파이썬 내가이 질문을 게시 한 후 답 4 분을 실현할 때 두 번째 종류

답변

1

사랑과 HADE의 전체 타원 적분을 계산 matematica. maple use m 파이썬과 matematica는 k를 사용합니다. 여기서 K^2 = m

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

enter image description here