#
#	Neuron type:				Medial LIII Complex Pyramidal (e)313300
#	Model ID:					6-006-3
#	Downloaded from:			http://hippocampome.org/
#	This script runs on XPP:	http://www.math.pitt.edu/~bard/xpp/xpp.html
#

# model parameters
par k=0.905351704053685
par a=0.0685965384362544
par b=4.075837348538556
par d=45.0
par Cm=115.0
par vr=-58.980533034841315
par vt=-35.20515631656637
par vmin=-38.942803187721466
par vpeak=21.98677890041239

# input current
par I=193.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)=-58.980533034841315
u(0)=0

# plotting variables
@ xp=t, yp=v, xlo=0, xhi=949.12, ylo=-68.98053303484132, yhi=31.98677890041239

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

done