Skip to content
Snippets Groups Projects
Commit 84ec4f53 authored by Sam Yates's avatar Sam Yates
Browse files

Use name `i_e` for Stim current density

parent 104e3897
No related branches found
No related tags found
No related merge requests found
...@@ -35,10 +35,10 @@ end ...@@ -35,10 +35,10 @@ end
immutable Stim immutable Stim
t0 # start time of stimulus t0 # start time of stimulus
t1 # stop time of stimulus t1 # stop time of stimulus
j # stimulus current density i_e # stimulus current density
Stim() = new(0s, 0s, 0A/m^2) Stim() = new(0s, 0s, 0A/m^2)
Stim(t0, t1, j) = new(t0, t1, j) Stim(t0, t1, i_e) = new(t0, t1, i_e)
end end
vtrap(x,y) = x/(exp(x/y) - 1.0) vtrap(x,y) = x/(exp(x/y) - 1.0)
...@@ -101,7 +101,7 @@ function f(t, state; p=HHParam(), stim=Stim()) ...@@ -101,7 +101,7 @@ function f(t, state; p=HHParam(), stim=Stim())
# calculate current density due to stimulus # calculate current density due to stimulus
if t>=stim.t0 && t<stim.t1 if t>=stim.t0 && t<stim.t1
itot -= stim.j itot -= stim.i_e
end end
# calculate the voltage dependent rates for the gating variables # calculate the voltage dependent rates for the gating variables
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment