尝试在 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])
我是一名优秀的程序员,十分优秀!