On 2/18/06, daemon <[EMAIL PROTECTED]> 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:[EMAIL PROTECTED] > > > > _______________________________________________ > Help-gsl mailing list > [email protected] > http://lists.gnu.org/mailman/listinfo/help-gsl > _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
