Question

尝试在枫树中实现Euler Elastica，Jacobi椭圆函数和椭圆积分至关重要。雅可比椭圆函数：JacobiSN，JacobiCN，JacobiDN和JacobiAM似乎按照我的预期工作，但椭圆积分：EllipticE表现为我认为应该做的事情。

对maple，matematica和python buildin函数进行了比较：

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

我唯一能找到的是： Maple计算第二类不完全和完整的椭圆积分。

其中python和matematica计算第二种完整的椭圆积分

Answer 1

在发布问题4分钟后我意识到答案时爱和恨。枫使用m，其中python和matematica使用k。其中k ^ 2 = m

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

enter image description here