#
#	Neuron type:				CA1 Oriens/Alveus (i)2233
#	Model ID:					4-055-4
#	Downloaded from:			http://hippocampome.org/
#	This script runs on XPP:	http://www.math.pitt.edu/~bard/xpp/xpp.html
#

# model parameters
par k=2.4478812416060833
par a=0.08988364891296882
par b=-29.880388216490868
par d=271.0
par Cm=98.0
par vr=-51.804767403318856
par vt=-32.28831673460901
par vmin=-36.46625123387067
par vpeak=4.125090119854121

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

# plotting variables
@ xp=t, yp=v, xlo=0, xhi=100.354, ylo=-61.804767403318856, yhi=14.125090119854121

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

done