Ch10_3.py

"""
Eqn (10.3) for calcium parabolic bursting
"""

from __future__ import division
from PyDSTool import *
from Ch10_3_model import *

gentype = 'radau' # 'dopri' # 'vode'

par_args = {'gr': 26, 'gc': 3.25, 'gx': 1.93,
            'tau_r': 5.6, 'tau_c': 125, 'tau_x': 300, 'tau_v': 0.97,
                'Ena': 0.48, 'Er': -0.95, 'Ec': -0.95, 'Ex': 1.4,
                'Iapp': 1.5, 'As': 0, 'f': 50, 'An': 0.0}
ic_args = {'v': -0.82, 'r': 0.4, 'x': 0.04, 'c': 0}


HH = makeHHneuron('HH_Ca', par_args, ic_args, const_I=True,
                  gentype=gentype)
HH.set(tdata=[0,1500])

def test_I(I):
    HH.set(pars={'Iapp': I})
    traj = HH.compute('test')
    pts = traj.sample()
    plt.figure(1)
    plt.clf()
    plt.plot(pts['t'], pts['v'], 'k')
    plt.plot(pts['t'], pts['r'], 'r')
    plt.plot(pts['t'], pts['x'], 'g')
    plt.plot(pts['t'], pts['c'], 'y')
    plt.xlim(-5,max(pts['t']))
    return traj, pts

traj, pts = test_I(1.6)

plt.show()