#
# A simple integrate-fire-burst model 
#
# Author: 
#
# Gregory D. Smith                    (757) 221-1989 office
# Department of Applied Science       (757) 221-2050 fax
# The College of William and Mary
# McGlothlin-Street Hall, Rm 305      greg@as.wm.edu
# Williamsburg, VA 23187              http://www.as.wm.edu/Faculty/Smith.html
# 
# Reference: 
#
# Smith GD, Cox CL, Sherman SM, and Rinzel J. Fourier analysis of
# sinusoidally-driven thalamocortical relay neurons and a minimal
# integrate-and-fire-or-burst model. J. Neurophys., 83(1):588-610, 2000.   
#

V(0)=-60
h(0)=0

p i0=-0.5,i1=1,freq=1,phi=0.0 
p c=2, gh=0.07,thm=20,thp=100

p gl=0.035,Vl=-65,Vreset=-50,Vspike=-35, Vh=-60, Vt=120

global 1 V-Vspike {V=Vreset}

i=i0+i1*cos(2*pi*t*freq/1000+phi)

dV/dt=(i-gl*(V-Vl)-gh*heav(V-Vh)*h*(V-Vt))/c
dh/dt=-h/thm*heav(V-Vh)-(h-1)/thp*heav(Vh-V)

@ total=3000,trans=0,method=rungekutta,dt=0.01,axes=2,Xlo=0,Xhi=3000,Ylo=-150,Yhi=50,Yp=V,maxstore=10000000,bounds=100000

done