Δημοσιεύτηκε: 29 Ιαν 2010, 12:45
από Dimitris
Off topic:
Είμαι τελειώς εκτός, το ξέρω, αλλά:

Δεν είναι πολύ δύσκολο να προγραμματίσεις στο octave/scilab/python/maxima ή οτιδήποτε άλλο για να λύσεις τις ode και μετά να κάνεις ένα phase plane plot. (Εγώ θα το κανα σε fortran :D ) Όλα τα παραπάνω έχουν λύτες για ode, και δε χρειάζεσαι γραφικό περιβάλλον για να τις λύσεις.

On the fly για το octave:
function xdot = f(x,t)
xdot=-exp(-t).*x.^2;
endfunction

x=lsode("f",2,(t=linspace(0,5,50)'));
plot(t,x)
dx = f(x,t)
figure
plot(x,dx)

Βάλε κι ένα loop για διαφορες αρχικές συνθηκες και τέλειωσες.