c3 = c2
c2 = c1
c1 = o
o = i
c1 = i
aa = 0.0655*exp(0.05547153*(-36.0 + v))
ik = Gk*(v - vk)*o
a = 0.0555*exp(0.05547153*(-12.0 + v))
b = 0.002357*exp(-0.036588*v)
Gk = 0.013500000000000002*ko^0.59
bi = 1.030072027074169*exp(0.000942*v)*(ko^(-1.0))^0.3
u = (ai*bb)/bi
bb = 0.0029357000000000003*exp(-0.02158*v)
vk = (R*Temp*log(ko/ki))/F
ai = 1.9755*exp(-0.02352*(25.0 + v))/ko
v = piecewise(vtest,50.0<time<=300.0,vhold)
Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.