I have a notebook to solve a differential equation using gsl's rkf45 algorithm. Now I realize that the equation is too stiff and rkf45 is may not be the right algorithm to use. The function I am trying to solve is written in C syntax. I want to try using octave's lsode to solve the same equation since it has a option of solving "stiff" equation. I don't know the jacobian of the equation so I can't use the stiff options in gsl. The problem is that my function is very complicated it took me a long time to type it out. So my question is - is it possible to solve a differential equation in octave where the rhs is written in C syntax.
asked Aug 18 '11Shashank
1570 ● 5 ● 22 ● 55
Before going overboard, how about you try the other algorithms that GSL supports? See http://www.sagemath.org/doc/reference/sage/gsl/ode.html
posted Aug 18 '11Volker Braun
2238 ● 5 ● 22 ● 52
Asked: Aug 18 '11
Seen: 108 times
Last updated: Aug 18 '11
powered by ASKBOT version 0.7.22
Copyright Sage, 2010. Some rights reserved under creative commons license.