Advanced Search

Question: jacobian of system of ODEs

Question:jacobian of system of ODEs

longrob 1409 Maple
ode differential-equations jacobian
January 24 2010
1
Greetings Is there a nice way to obtain the jacobian of a system of ODEs ? For example with restart: ross_x:=diff(x(t),t)=-y(t)-z(t): ross_y:=diff(y(t),t)=x(t)+a*y(t): ross_z:=diff(z(t),t)=b+x(t)*z(t)-c*z(t): rossler_sys:=ross_x,ross_y,ross_z; I can't seem to find an easy way to do it, other than re-write the system altogether s_x:=-y-z: s_y:=x + a*y: s_z:=b+z*(x-c): J:=Matrix([[diff(s_x,x),diff(s_x,y),diff(s_x,z)],[diff(s_y,x),diff(s_y,y),diff(s_y,z)],[diff(s_z,x),diff(s_z,y),diff(s_z,z)]]); .....which seems horrible. Thanks LR
Please Wait...