I am trying to simulate the effect of bipolar electric stimulation at a single neuron. I use (1) the 1/distance-rule in order to calculate the extracellular voltage induced by each pole, (2) the built-in function
Code: Select all
The above simulation seems to work fine when using a constant time step (Dt) method with a relatively small Dt=0.05 ms. (Note that the extracellular stimulation is always < 200 Hz.) When I try to use CVODE, Vm doesn't seem to be affected by the extracellular stimuli at all. This, even when I set atol < 1e-7 for all STATE variables in the kinetic mechanisms I use (for the intracellular Ca-mechanism I even use 1e-10).
Based on the above I am wondering whether I am missing something... For any help suggestions I would be most thankful!