枫树中的椭圆形与mathematica和python相比

时间:2015-02-04 08:48:06

标签: python math wolfram-mathematica integration maple

尝试在枫树中实现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计算第二种完整的椭圆积分

1 个答案:

答案 0 :(得分:1)

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

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

enter image description here