尝试在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])