Numerical simulation of solitons in the nerve axon using finite differences

A High-order accurate finite difference scheme is derived for a non-linear soliton model of nerve signal propagation in axons. Boundary conditions yielding well-posed problems are suggested and included in the scheme using a penalty technique. Stability is shown using the summation-by-parts framework for a frozen parameter version of the non-linear problem