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

[Analabha Roy]

---------- Forwarded message ----------
From: Analabha Roy <address@hidden>
Date: Feb 18, 2006 3:03 PM
Subject: Re: [Help-gsl] solving linear n-order ode
To: "st." <address@hidden>

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>

On 2/18/06, st. <address@hidden> wrote:
>
>
> 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)
>
> thanks in advance
> --
> 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
> address@hidden
> http://lists.gnu.org/mailman/listinfo/help-gsl
>