Trama a odi, isocline usando Python

Question

Sto cercando un pacchetto Python che mi permetta di tracciare qualcosa di simile all'applet Java visto di seguito:

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

Qualcuno conosce qualche pacchetti che tracciano ODE per questo? Posso codificare qualcosa da zero usando Numpy, Matplotlib, ma volevo prima chiedere in giro.

Grazie,

Answer 1

Ho scritto qualcosa del genere, sembra funzionare per 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 lo farà:

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)

Anche se sta semplicemente avvolgendo la funzionalità di matplotlib per la grafica. (Ad essere perfettamente onesti, la confezione di matplotlib non è così buona come potrebbe ancora essere, il che spesso mi provoca mal di testa.)

Answer 3

Queste risposte non hanno la possibilità di modificare i parametri utilizzando lo strumento di trascinamento. Se vuoi questa opzione, questi due sistemi dinamici di esempio ti mostrano come. Sono scritti in Python Sage. Pensalo come Python con molte funzioni pre-fatte per la matematica.

Sage Esempio 1-Trama della fase.
Sage Esempio 2-Trama traiettoria.