[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: matlab to Octave conversion

From: Kevin H. Hobbs
Subject: Re: matlab to Octave conversion
Date: Wed, 14 Jun 2006 15:42:26 -0400

On Wed, 2006-06-14 at 19:49 +0100, MICHAEL FORREST wrote:
> The model is from this paper. The author sent it to me to play around with. 
> The model is for Figure 8 (a) of the paper. It is all pretty complex and I 
> am looking into it - because I'm very interested in purkinje bistability. Do 
> you think that the model is flawed?

No I do not. What scared me at first, was the concatenation of current
clamp and voltage clamp model outputs into the same array and then
plotting them together on the same graph which is what I think is going
on. One would expect current clamp to measure voltage and voltage clamp
to measure current. Plotting them together wouldn't make much sense
unless you want a nice plot of your command signal.

In cell_dynamics_clamp.m you can see that the voltage derivative is set
to 0 ( voltage clamp ) only for times less than 0.05, which is
presumably enough time to switch limit cycles, but then it looks like
the model returns to current clamp. 

The return to current clamp is the second thing that scared me, because
it might throw you off if you weren't expecting it.

>  Thanks so much for taking the time with 
> me. I am so grateful.

No problem. You should try the changes I suggested. It took only a few
seconds and octave produced plots that are very much like those in
figure 8 of that paper.

Attachment: signature.asc
Description: This is a digitally signed message part

reply via email to

[Prev in Thread] Current Thread [Next in Thread]