Python绘制ODE，使用等值线

Question

我正在寻找一个Python包，它允许我绘制类似于下面所示的Java applet的内容：

http://math.mit.edu/mathlets/mathlets/isoclines/

有人知道有什么颂歌绘图包吗？我可以用Numpy，Matplotlib从头开始编写代码，但我想先问问周围的人。在

谢谢

Answer 1

Sage会这样做：

x,y = var("x y")
eq = y^3-3*y-x
p = implicit_plot(eq==0,(x,-4,4),(y,-4,4))
p += plot_slope_field(eq, (x,-4,4),(y,-4,4), headlength=1e-8)
p.show(aspect_ratio=1)

尽管它只是包装了图形的matplotlib功能。（老实说，matplotlib包装还没有达到预期的效果，这常常让我头疼。）

Answer 2

我写了这样的东西，似乎对y'=y^2-x有效

from pylab import *
xmax = 4.0
xmin = -xmax
D = 20
ymax = 4.0
ymin = -ymax
x = linspace(xmin, xmax, D)
y = linspace(ymin, ymax, D)
X, Y = meshgrid(x, y)
deg = arctan(Y**2 - X)
QP = quiver(X,Y,cos(deg),sin(deg))
show()

Answer 3

这些答案不能使用拖动工具更改参数。如果你想要这个选项，这两个例子动力系统会告诉你如何。它们是用Python Sage编写的。你可以把它想象成Python，里面有很多数学函数。在

Sage Example 1 phase plot。
Sage Example 2 trajectory plot。在