#
#	Neuron type:				CA1 BC CCK+ (i)2232
#	Model ID:					4-079-1
#	Downloaded from:			http://hippocampome.org/
#	This script runs on XPP:	http://www.math.pitt.edu/~bard/xpp/xpp.html
#

# model parameters
par k=1.062999344963705
par a=0.023084662044171957
par b=-26.86555078015733
par d=294.0
par Cm=59.0
par vr=-63.183009204776944
par vt=-38.18266358350243
par vmin=-53.318562041665636
par vpeak=27.443037301495828

# input current
par I=333.0

# model definition
global 1 {v-vpeak}{v=vmin;u=u+d}
v'=(k*(v-vr)*(v-vt)-u+I)/Cm
u'=a*(b*(v-vr)-u)

# initial condition
v(0)=-63.183009204776944
u(0)=0

# plotting variables
@ xp=t, yp=v, xlo=0, xhi=933.72, ylo=-73.18300920477694, yhi=37.44303730149583

# integration parameters
@ meth=rk,dt=0.1, total=933.72, Bounds=1000

done