尝试在 maple 中实现 Euler Elastica，其中 Jacobi 椭圆函数和椭圆积分至关重要。 Jacobi 椭圆函数:JacobiSN、JacobiCN、JacobiDN 和 JacobiAM 似乎按我的预期工作，但椭圆积分:EllipticE 的行为与我认为它应该做的不准确。
比较了maple、matematica和python的内置函数:
Maple:
plot([t, EllipticE(t), t = 0 .. 1])
Matematica:
plot[EllipticE[t], {t, 0, 1}]
Python(scipy.special):
t = arange(0, 1, 0.001)
plot(t, ellipe(t))
我唯一能找到的是: Maple计算第二类不完全和完全椭圆积分。
其中python和matematica计算第二类完全椭圆积分
请您参考如下方法:
当我在发布问题 4 分钟后意识到答案时，爱和哈迪。 maple 使用 m，而 python 和 matematica 使用 k。其中 k^2 = m
plot([t, EllipticE(sqrt(t)), t = 0 .. 1])
