[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## Re: Re[2]: [Help-gsl] solving linear n-order ode

**From**: |
Analabha Roy |

**Subject**: |
Re: Re[2]: [Help-gsl] solving linear n-order ode |

**Date**: |
Sat, 18 Feb 2006 15:47:41 -0600 |

On 2/18/06, daemon <address@hidden> wrote:
>
>* Hello Analabha,*
>
>* Saturday, February 18, 2006, 11:03:47 PM, you wrote:*
>
>* > See chapter on ODEs in gsl manual. Reduce ur ODE to*
>* > y_i`(x)=g_i(x,y_k) for all i*
>
>* If I've understood you right, I need to reduce my ode _MANUALLY_, but*
Well you can just do that once, right? (I mean, write down a general
transformation rule from your form to theat understandable by gsl)
it's not that I really need. I need to make a solver wich will solve
>* any linear ode. And only information I'll need to pass to solver is an*
>* order of equation, array of coefficients and a function g(x). Do you*
>* have any suggestions to solve this problem?*
Yeah, so you create a subroutine for g(x), & a subroutine that integrates
the ODE. This suboutine will accept the order, the array of coefficients and
then applies your transformation rule and make it gsl-understandable func()
jac() funcions (you'll have to write those), then initialize the ODE system
& evolve_apply.
I do something similar when I solve Schroedinger's equation with complex
variables instead of real #s in N-space representation or using FEM .
--
>* Best regards,*
>* daemon mailto:address@hidden*
>
>