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 nafpr USEION na WRITE ina RANGE gna, ina } UNITS { (mA) = (milliamp) (mV) = (millivolt) (mS) = (millisiemens) } PARAMETER { gna = 30 (mS/cm2) ena = 55 (mV) } ASSIGNED { v (mV) ina (mA/cm2) minf (1) hinf (1) tauh (ms) } STATE { h } INITIAL { rates(v) h = hinf } BREAKPOINT { SOLVE states METHOD cnexp ina = (1e-3) * gna * minf^2 * h * (v-ena) } DERIVATIVE states { rates(v) h' = (hinf-h)/tauh } :ina PROCEDURE rates(v(mV)) { LOCAL a, b a = fun3(v, -46.9, -0.32, -4) b = fun3(v, -19.9, 0.28, 5) minf = a/(a+b) a = fun1(v, -43, 0.128, -18) b = fun2(v, -20, 4, -5) hinf = a/(a+b) tauh = 1.0/(a+b) } INCLUDE "aux_fun.inc"