help-octave
[Top][All Lists]

## Re: Octave coding

 From: Shai Ayal Subject: Re: Octave coding Date: Fri, 16 Jun 2006 09:22:19 +0300 User-agent: Thunderbird 1.5.0.4 (Windows/20060516)

```Michael,

```
It seems like you are asking for someone to write code for you. This is beyond the scope of this list, in which *volunteers* offer help to people which have trouble coding. Perhaps if you were willing to pay, someone would be willing to invest the time and write the code for you.
```
Shai

MICHAEL FORREST wrote:
```
```Dear all,

```
I would like Octave code to implement the following equation system and to plot (v) vs. time. However, I am at a complete loss how to do this. Could anyone help? I would be so, so grateful. I am under quite a bit of time pressure and haven't the time to get my Octave programming up to the heights required to do this task in time. Please don't hesitate to write back if you require any more info. I am using Octave as a module of cygwin on a Windows XP machine. I have the GNUPlot module installed as well as Cygwin.
```Once again. Thank you so, so much. All the best,
Michael

NB.
v' denotes dv/dt
b' denotes db/dt
h' denotes dh/dt

-------------

// Master equation

v' = -(ina+ik+ih+il)

// ina

ina  = 0.000060*minf*(v-55)
minf = (1+exp(-(v--53.8)/3))^-1

// ik

ik  = 0.000105*b*(v--85)
b' = ((binf - b)/btau)
binf = (1+exp(-(v--54)/5))^-1
btau = (3000 * (1/(cosh((v--54)/(4*5)))))

// ih

ih  = 0.000200*h*(v--30)
h' = ((hinf - h)/htau)
hinf = (1+exp((v--76.4)/20))^-1
```
htau = (1/((((-2.89 * v) + -445)/(1-exp(((v+-445)/-2.89)/24.02))) + (((27.1 * v) + -1024)/(1-exp(((v+-1024)/27.1)/-17.4)))))
```
// Ileak

il = 0.000100*(v--70)

--------------

_______________________________________________
Help-octave mailing list
```