COMMENT //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% // // NOTICE OF COPYRIGHT AND OWNERSHIP OF SOFTWARE // // Copyright 2007, The University Of Pennsylvania // School of Engineering & Applied Science. // All rights reserved. // For research use only; commercial use prohibited. // Distribution without permission of Maciej T. Lazarewicz not permitted. // mlazarew@seas.upenn.edu // //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ENDCOMMENT NEURON { SUFFIX icapr USEION ca WRITE ica RANGE gca, ica } UNITS { (mA) = (milliamp) (mV) = (millivolt) (mS) = (millisiemens) } PARAMETER { gca = 10 (mS/cm2) eca = 80 (mV) } ASSIGNED { ica (mA/cm2) v (mV) sinf (1) taus (ms) } STATE { s } INITIAL { rates(v) s = sinf } BREAKPOINT { SOLVE states METHOD cnexp ica = (1e-3) * gca * s^2 * (v-eca) } DERIVATIVE states { rates(v) s' = (sinf-s)/taus } PROCEDURE rates(v(mV)) { LOCAL a,b a = fun2(v, 5, 1.6, -1/0.072) b = fun3(v, -8.9, 0.02, 5) sinf = a/(a+b) taus = 1.0/(a+b) } INCLUDE "aux_fun.inc"