2013-06-04 10 views
14

Ho lavorato con sympy e scipy, ma non riesco a trovare o capire come risolvere un sistema di equazioni differenziali accoppiate (non lineare, primo ordine).Un modo per risolvere un sistema di equazioni differenziali accoppiate in python?

Quindi esiste un modo per risolvere equazioni differenziali accoppiate?

Le equazioni sono della forma:

V11'(s) = -12*v12(s)**2 
v22'(s) = 12*v12(s)**2 
v12'(s) = 6*v11(s)*v12(s) - 6*v12(s)*v22(s) - 36*v12(s) 

con condizioni iniziali per v11 (s), V22 (s), V12 (s).

+0

Dai uno sguardo a [saggio] (http://www.sagemath.org/). Offre funzionalità mathmatica con la sintassi python. Potrebbe essere in grado di risolvere i diff eq. – SethMMorton

+0

Stai cercando una soluzione analitica o una soluzione numerica? (Hai menzionato l'utilizzo di sympy, quindi potresti sperare in una soluzione analitica, se ce n'è una.) –

+0

@WarrenWeckesser Una soluzione numerica, simile a NDsolve per mathematica. – faceforest

risposta

12

Per la soluzione numerica di ODE con scipy, consultare la funzione scipy.integrate.odeint o la classe scipy.integrate.ode.

Alcuni esempi sono riportati nello SciPy Cookbook (scorrere fino alla sezione "Equazioni differenziali ordinarie").

+0

ha funzionato perfettamente ed è stato facile da seguire. Grazie! – faceforest

Problemi correlati