[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] ODE solver: oscillation around expected value

From: Jeroen Meidam
Subject: [Help-gsl] ODE solver: oscillation around expected value
Date: Mon, 01 Aug 2011 13:08:50 +0100
User-agent: Mozilla/5.0 (Windows; U; Windows NT 6.1; en-GB; rv: Gecko/20110616 Lightning/1.0b2 Thunderbird/3.1.11


I'm using the gsl ode solver to solve a system of 4 coupled differential equations. I've done the exact same thing in python, but needed to rewrite it in c. There is one component of the plotted solutions that shows an unexpected oscillatory behaviour around the value that python gives as aresult. I have no idea what could cause this. I have tried all the available stepping functions and none appear to get rid of this oscillation.

Also, I would expect more differences between the different solving methods. If I plot for example the rk2 , rk4, rkf45 and rk8pd solutions in one graph, I can't see the slightest bit of difference. Is that to be expected?

reply via email to

[Prev in Thread] Current Thread [Next in Thread]