v0         = -> Ca, 1
v1         = -> Ca, \
                 Bstep*k1
    
Bstep      = 0.4
#Bstep     = forcing(step)
k1         = 7.3
B          = 0.4
t_stimulus = 1.0
    
export     = Ca ->  , 10 ..
leak       = CaComp -> Ca, 1 ..
    
v2         = Ca -> CaComp, \
                  65 * Ca**2 / (1+Ca**2)
    
v3         = CaComp -> Ca, \
                  500*CaComp**2/(CaComp**2+4) * Ca**4/(Ca**4 + 0.6561)
    
init       = state(Ca = 0.1, CaComp = 0.63655)
