在做关于外界刺激与神经细胞电位发放的关系的论文时不知道如何编程得到Hodgkin-Huxley方程的图像。由于这里不能进行公式编辑，方程具体内容参见附件。请各位大虾帮忙，急，谢谢！

= — gNam3h(V—VNa) — gkn4(V—Vk) — gL(V—VL) + Idc + Ip + D (t)

= — (am + bm)m + am

= — (ah + bh)h + ah

= — (an + bn)n + an

其中am = ah = 0.07e—0.05(V+65)

an = bm = 4e—(V+65) / 18

bh = bm = 0.125e—(V+65) / 80

V(t)表示膜电位，对应神经细胞的电位发放。m(t),h(t),n(t)分别描述细胞膜内外的钠离子，钾离子和泄露电流。gNa= 120ms/cm2, gk = 36ms/cm2, gL = 0.3ms/cm2分别对应钠离子，钾离子和泄露电流关于细胞膜的电导系数的最大值，参数电位
VNa = 50mV, Vk = —77mV, VL = —54.5mV分别对应模型(1)中的反向电压。

[此贴子已经被作者于2007-3-20 0:31:49编辑过]

看不懂啊
附件在那里

重新编辑了一遍，可以看到附件了，呵呵

参考一下ode系列算法，很简单的
Example
[t,y]=ode45(@vdp1,[0 20],[2 0]);
plot(t,y(:,1));
solves the system y' = vdp1(t,y)
可以令y=[Dv Dm Dh Dn]'

谢谢hitzhang的提示加上搜索本论坛的相关信息终于画出了这个图，呵呵

请问你是怎么解的啊，能不能告诉我啊！