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

[Ben Bolker]

Berend Hasselman <bhh <at> xs4all.nl> writes:

> 
> 
> On 04-07-2013, at 17:15, Jannetta Steyn <jannetta <at> henning.org> wrote:
> 
> > Hi folks
> > 
> > I have implemented a model of a neuron using Hodgkin Huxley equations in
> > both R and MatLab. My preference is to work with R but R is not giving me
> > the correct results. I also can't use ode45 as it just seems to go into an
> > indefinite loop. However, the MatLab implementation work fine with
> >  the same
> > equations, parameters and initial values and ode45. Below is my R and
> > MatLab implementations.
> > 
> 
> No problem in running your R file. Have plot.
> (Mac mini Core2Duo 2.66Ghz running R 3.0.1 Patched 
> (2013-06-19 r62992) on Mac OS X 10.8.4:  16.5 seconds.)
> 
> Trying to run your Matlab file in Octave.
>  No success yet because of unavailable ode45.

   I'm impressed that you (BH) went to the trouble of checking
on this vague "doesn't work" question.  If you want to go farther
I think you can get ode45 for octave by installing the odepkg
package:

http://octave.sourceforge.net/odepkg/index.html