Trazar ODES, isoclinas usando Python

Question

Estoy buscando un paquete de Python que me permita trazar algo similar al applet Java que se ve a continuación:

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

¿Alguien sabe algún paquete de trazado de oda para esto? Puedo codificar algo desde cero usando numpy, matplotlib, pero quería preguntar primero.

Gracias,

Answer 1

Escribí algo como esto, parece funcionar para 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

Sage hará esto:

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)

Aunque es simplemente envolver la funcionalidad matplotlib para los gráficos. (Para ser perfectamente honesto, el envoltorio de matplotlib no es tan bueno como podría ser, lo que a menudo me causa dolores de cabeza).

Answer 3

Estas respuestas no tienen la opción de cambiar los parámetros utilizando la herramienta de arrastre. Si desea esta opción, estos dos sistemas dinámicos de ejemplo le muestran cómo. Están escritos en Python Sage. Piense en ello como Python con muchas funciones prefabricadas para las matemáticas.

SAGE Ejemplo 1-Parcela de fase.
SAGE Ejemplo 2-Gráfico de trayectoria.