On 5 July 2013 09:44, Berend Hasselman <b...@xs4all.nl> wrote: > > On 05-07-2013, at 09:53, Jannetta Steyn <janne...@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: janne...@henning.org =================================== [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
