Numerical solution of a nonlinear advance-delay-differential equation from nerve conduction theory

A functional differential equation which is nonlinear and involves forward and backward deviating arguments is solved numerically. The equation models conduction in a myelinated nerve axon in which the myelin completely insulates the membrane, so that the potential change jumps from node to node. The equation is of first order with boundary values given at t = +/- infinity. The problem is approximated via a difference scheme which solves the problem on a finite interval by utilizing an asymptotic representation at the endpoints, cubic interpolation and iterative techniques to approximate the delays, and a continuation method to start the procedure. The procedure is tested on a class of problems which are solvable analytically to access the scheme's accuracy and stability, then applied to the problem that models propagation in a myelinated axon. The solution's dependence on various model parameters of physical interest is studied. This is the first numerical study of myelinated nerve conduction in which the advance and delay terms are treated explicitly.