On 5 July 2013 09:44, Berend Hasselman <bhh@xs4all.nl> wrote:

>
> On 05-07-2013, at 09:53, Jannetta Steyn <jannetta@henning.org> wrote:
>
> >
> >
> > >
> > > I don't quite know how to explain the "doesn't work" in more detail
> without
> > > any visual aid.
> >
> > You said that R got into an indefinite loop, whatever that maybe.
> >
> >
> > > When I change the solver to ode45 the script never stops running. I
> have even left it over night at one stage but the next day it was still
> busy. Is there a way to see what it is doing and to determine why it seems
> to be in an
> > > "infinite loop"?
> >
> > The script just ran using ode45!! For the first time ever.
> >
>
> Please keep the conversation on R-help.
> Don't reply to me personally.
>


Apologies for that. It was meant to go to the list.


>
> Which script?
> With the corrections suggested by me and Peter? (gK_axon_AB=52.5-3 ==>
> gK_axon_AB=52.5e-3)
>
>
The R script. With the correction of AB=52.5e-3



> For what it's worth: lsdoda seems quicker.
>

Yes, lsoda is quicker but I want to use ode45 because that is what they
used in the paper I am referencing and I am trying to get the same results
as they did.


> Variable v_acon_AB now converges to -60 (the equilibrium state?)
>

I'm trying to determine now whether this is correct.

Regards
Jannetta


-- 

===================================
Web site: http://www.jannetta.com
Email: jannetta@henning.org
===================================

	[[alternative HTML version deleted]]

