We use high-order approximation schemes for the space derivatives in the nonlinear cable equation and investigate the behavior of numerical solution errors by using exact solutions, where available, and grid convergence. The space derivatives are numerically approximated by means of differentiation matrices. A flexible form for the injected current is used that can be adjusted smoothly from a very broad to a narrow peak, which leads, for the passive cable, to a simple, exact solution. We provide comparisons with exact solutions in an unbranched passive cable, the convergence of solutions with progressive refinement of the grid in an active cable, and the simulation of spike initiation in a biophysically realistic single-neuron model.
Model Type: Neuron or other electrically excitable cell
Simulation Environment: MATLAB
Implementer(s): Omurtag, Ahmet [aomurtag at gmail.com]