       Help me

• To: mathgroup at yoda.physics.unc.edu
• Subject: Help me
• From: neurox.cu!bosch at ceniai.cu (Jorge Bosch)
• Date: Thu, 11 Feb 93 11:30:55 -0500

```Hello, Mathgroup:

I am new at Mathematica.

I'm introducing this framework in my scientific group but I'm having
trouble with it.

I have installed 2.1 version of Mathematica under Windows, on an Acer 1100sx
computer (16 MHz, 4 MBytes Extended Memory, 80387 Numerical Coprocessor),
MSDOS Operating System.

I am trying to solve a differential equation system (the Hoodgking and Huxley
model for a neuron) but different solutions appear at different times and
sometimes the computer crashes. Additionally, I can't reproduce the solution
that I got with the same system programmed in Turbo Pascal.

I enclose a copy of the non autonomous and non linear differential
equation system (four equations) as I am using. I will really appreciate if
someone could help me to find out what is wrong in my system.

Jorge Bosch

a1[v_]:=0.1 (25-v)/(Exp[(25-v)/10]-1)
b1[v_]:=4 Exp[-v/18]
a2[v_]:=0.07 Exp[-v/20]
b2[v_]:=1/(Exp[(30-v)/10]+1)
a3[v_]:=0.01 (10-v)/(Exp[(10-v)/10]-1)
b3[v_]:=0.125 Exp[-v/80]

c1oo[v_]:=a1[v]/(a1[v]+b1[v])
c2oo[v_]:=a2[v]/(a2[v]+b2[v])
c3oo[v_]:=a3[v]/(a3[v]+b3[v])

i[t_]:=If[t>=0.5 && t<=7.55,130,0]

rm=35.4;
cm=1;
gnax=120;
gkx=36;
glx=0.3;
vna=115;
vk=-12;
vl=11;

Plot[i[t],{t,0,10},PlotPoints->100]

NDSolve[{v'[t]==-(gnax m[t]^3 h[t] (v[t]-vna) +
gkx n[t]^4 (v[t]-vk) + glx (v[t]-vl) + v[t]/rm - i[t]),
m'[t]==a1[v[t]](1-m[t])-b1[v[t]] m[t],
n'[t]==a2[v[t]](1-n[t])-b2[v[t]] n[t],
h'[t]==a3[v[t]](1-h[t])-b3[v[t]] h[t],
v==0,
m==c1oo,
n==c2oo,
h==c3oo},
{v[t],n[t],m[t],h[t]},{t,0,10}]

Plot[Evaluate[v[t]/.%],{t,0,10},PlotPoints->100];

```

• Prev by Date: Re: NeXT and Cannon
• Next by Date: QR for upper Hessenberg matrices
• Previous by thread: NeXT and the NeXT of us!
• Next by thread: RE: Help me