Pythonを使用したODES、ISOCLINESのプロット

Question

以下に示すJavaアプレットに似たものをプロットできるようにするPythonパッケージを探しています。

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

誰もがこのためにパッケージをプロットすることを知っていますか？ numpy、matplotlibを使用してゼロから何かをコーディングできますが、最初に尋ねたかったのです。

ありがとう、

Answer 1

私はこのようなものを書いた、それは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 2

セージはこれを行います：

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 3

これらの回答には、ドラッグツールを使用してパラメーターを変更するオプションがありません。このオプションが必要な場合は、これらの2つの例の動的システムがその方法を示します。それらはPython Sageで書かれています。数学のための多くの既製の機能を備えたPythonと考えてください。

セージの例1-フェーズプロット.
SAGE例2--トラムプロット.