help-gsl
[Top][All Lists]

## 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):

http://www.ph.utexas.edu/~daneel/floquet_particle_box_diag_eigenmatrix.c

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

Hope this helps

--
> Best regards,
>
>
>
> _______________________________________________
> Help-gsl mailing list