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

[Analabha Roy]

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,
> daemon                            mailto:address@hidden
>
>
>
> _______________________________________________
> Help-gsl mailing list
> address@hidden
> http://lists.gnu.org/mailman/listinfo/help-gsl
>