help-gsl
[Top][All Lists]

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

 From: Analabha Roy Subject: Re: [Help-gsl] solving linear n-order ode Date: Sat, 18 Feb 2006 15:41:32 -0600

---------- Forwarded message ----------
Date: Feb 18, 2006 3:03 PM
Subject: Re: [Help-gsl] solving linear n-order ode

See chapter on ODEs in gsl manual. Reduce ur ODE to

y_i`(x)=g_i(x,y_k) for all i, then write functions for evaluating g_i &
delG_i/dely_k (jacobian) (you can put your a_n 's in an array & call them by
reference using pointers in your subroutines), then inititlize the gsl_odeiv
datatypes & use the gsl_odeiv_evolve_apply function to integrate your ODE>

>
>
> Hi. Could anyone help me to find out how to solve n-order linear ode with
> gsl. Is it possible with gsl? In general case I need to solve sush
> equation:
> any(n)(x)+an-1y(n-1)(x)+...+a2y''(x)+a1y'(x)+a0y(x)=bf(x)
>
> --
> View this message in context:
> http://www.nabble.com/solving-linear-n-order-ode-t1145994.html#a3004600
> Sent from the Gnu - Gsl - Help forum at Nabble.com.
> _______________________________________________
> Help-gsl mailing list