[R] R and MatLab implementations of the same model differs

Berend Hasselman bhh at xs4all.nl
Fri Jul 5 10:44:28 CEST 2013

On 05-07-2013, at 09:53, Jannetta Steyn <jannetta at 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.

Which script?
With the corrections suggested by me and Peter? (gK_axon_AB=52.5-3 ==> gK_axon_AB=52.5e-3)

For what it's worth: lsdoda seems quicker.
Variable v_acon_AB now converges to -60 (the equilibrium state?)


More information about the R-help mailing list