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

[Berend Hasselman]

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?)

Berend