[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: lsode question
From: |
Tatsuro MATSUOKA |
Subject: |
Re: lsode question |
Date: |
Wed, 26 Sep 2007 05:00:03 +0900 (JST) |
--- Rupert Swarbrick <address@hidden> wrote:
> Hi!
>
> I'm calculating several functions governed by coupled differential
> equations (5 in total), but the first has delta functions all over the
> place and doesn't depend on the others and so is much easier to sort out
> with a custom routine.
>
> Thus I've now got 4 coupled ODEs plus a function that I've calculated
> numerically.
>
> What I'd like to do is use lsode to calculate the other 4 instead of
> the for loop I wrote originally, but my function which gives the
> derivative of the system somehow needs to get hold of the value of the
> first function. Although I can pass its values through as an array,
> this isn't much help, as I'm given "t" not "k": the position in the
> array of t values we are currently at.
>
> Is there a nice solution? At the moment, I'm using the fact that my
> values of t were calculated via linspace and working out k from that,
> but it feels like quite a kludge.
Can you tell us an example of your code or a problem?
At least for me, I cannot fully understand your problem the above description.
Tatsuro
--------------------------------------
Easy + Joy + Powerful = Yahoo! Bookmarks x Toolbar
http://pr.mail.yahoo.co.jp/toolbar/