Fix to set critical points in integration: start and end of current pulse

Author: pgleeson

Tutorial/Source/HodgkinHuxley.py

@@ -246,7 +246,7 @@ def simulate(self, init_values=[-64.99584, 0.05296, 0.59590, 0.31773]):
         """
 
         # init_values are the steady state values for v,m,h,n at zero current injection
-        X = odeint(self.dALLdt, init_values, self.t, args=(self,))
+        X = odeint(self.dALLdt, init_values, self.t, args=(self,), tcrit=[self.I_inj_delay, self.I_inj_duration + self.I_inj_delay])
         V = X[:,0]
         m = X[:,1]
         h = X[:,2]