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"