An investigation into the influence of boundary condition specification in finite difference methods on the behaviour...

Neuronal models provide a major aid to understanding the behaviour of individual neurons and networks of neurons. The solution of the model equations by finite difference methods is widespread because of the inherent simplicity of the technique. Error in the finite difference approach due to spatial and temporal discretisation is shown to be equivalent to a mis-specification of membrane current density. The effect of this mis-specification on the accuracy of the solution to the model equations is shown to depend on the structure of the model and its input, as well as the size of the discretisation intervals themselves. Through a theoretical analysis, illustrated by a number of examples on passive and active dendrites, this article demonstrates that the accuracy with which core current is implemented numerically at segment end-points in elementary models influences the behaviour of the numerical solution of these models, and consequently any physiological conclusions drawn from them.