[Top][All Lists]

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

Re: Re[4]: [Help-gsl] solving linear n-order ode

From: Analabha Roy
Subject: Re: Re[4]: [Help-gsl] solving linear n-order ode
Date: Sat, 18 Feb 2006 18:33:30 -0600

On 2/18/06, daemon <address@hidden> wrote:
> Hello Analabha,
> Saturday, February 18, 2006, 11:47:41 PM, you wrote:
> > Well you can just do that once, right? (I mean, write down a general
> > transformation rule from your form to theat understandable by gsl)
>   Yeah, but I'm not deal with "differential equations theory" so much
>   to do this :((( I don't belive that there is no free library in the
>   world, that has abilities to solve n-order linear ode, given in
>   general form?

 It seems to me that you want to do a logical reduction of your generic ODE
system before solving it, & don't wanna do it by hand. I think what you need
is a symbolic computational platform like mathematica to do that. It's
proprietary, though. There is a FOSS symbolic compu platform called "maxima
or macsyma" (Google for it) but I haven't used it so I can't attest to
whether it does what you want (mathematica does).

 Anywho, as far as gsl is concerned, I've been using it to solve nth order
non-linear ODE's without any problems (with general n which I can input thru
scanf and take array elements accordingly). I didn't do anything terribly
fancy; merely read the relevant chapter in the manual , dumbly copied their
example(s) and tweaked it for my problem. Seems to be workin' so far. If you
need to see my working code (Solving Scroedinger's equation for a particle
in a box with a temporally harmonic drive in 40 dimensions), it's here
(ignore the essl/zgeev() part of the code, that's just to diagonalize the
non-hermitian Floquet matrix and is specific to my problem only, & also the
#pragmas, they are for parallelization):

 The main gsl stuff is in the subroutines integrate() func() & jac().

Hope this helps

> Best regards,
> daemon                            mailto:address@hidden
> _______________________________________________
> Help-gsl mailing list
> address@hidden

reply via email to

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