Построение ODE, Isoclines с использованием Python

Question

Я ищу пакет Python, который позволил бы мне построить что-то похожее на Java-апплет, показанный ниже:

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

Кто-нибудь знает какие-либо пакеты для печати ODE для этого? Я могу написать что-то с нуля, используя 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

В этих ответах нет возможности изменить параметры с помощью инструмента перетаскивания. Если вы хотите эту опцию, то эти два примера динамических систем покажут вам, как. Они написаны на Python Sage. Просто представьте, что это Python с множеством готовых функций для математики.

Sage Пример 1 - фазовый график .
Sage Example 2 - траектория построения .