[Top][All Lists]

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

Re: [Help-gsl] Solving a First Order ODE system using GSL

From: Horvath Krisztian
Subject: Re: [Help-gsl] Solving a First Order ODE system using GSL
Date: Tue, 10 Nov 2009 13:15:47 -0500
User-agent: Thunderbird (X11/20090626)

Earl, Joseph írta:

I am trying to solve a first-order ODE of the form

x' = sqrt(1/3x - k + lx^2)

using the GSL C++ library, where k and l are constants for a particular 

I have no problem at all using GSL to solve 2nd order (such as x'' = -x) and 
higher ODEs, something is just confusing about the first order equations and 
how to put them into the required form.
The Jacobian doesn't seem to make much sense with a 1st order ODE.

Also the current documentation is pretty poor - there is only one example of a 
2nd order ODE which is repeated all over the net. So if anyone wants to add a 
1st order example to the help docs I'm sure that would be of some use to others 
(and me!) in the future.

Any help would be much appreciated (if you feel like writing out a full working 
example that would be awesome).

Joseph Earl
Imperial College London

I am also interested in this topic. I would like to use GSL's ODE solver, but I was not able to implement it for first order equations either. So, I also thank any answer for this question.

Best regards,
Krisztian Horvath

reply via email to

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