尝试在枫树中实现Euler Elastica,Jacobi椭圆函数和椭圆积分至关重要。 雅可比椭圆函数:JacobiSN,JacobiCN,JacobiDN和JacobiAM似乎按照我的预期工作,但椭圆积分:EllipticE表现为我认为应该做的事情。
对maple,matematica和python buildin函数进行了比较:
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计算第二种完整的椭圆积分
答案 0 :(得分:1)
在发布问题4分钟后我意识到答案时爱和恨。 枫使用m,其中python和matematica使用k。其中k ^ 2 = m
plot([t, EllipticE(sqrt(t)), t = 0 .. 1])