(Neural Computation. 2003;15:2129-2146.)
© 2003 The MIT Press
A Simple and Stable Numerical Solution for the Population Density Equation
M. de Kamps
kamps{at}in.tum.de, Section Cognitive Psychology, Faculty of Social Sciences, Leiden University, 2333 AK Leiden, The Netherlands
A population density description of large populations of neurons has generated considerable interest recently. The evolution in time of the population density is determined by a partial differential equation (PDE). Most of the algorithms proposed to solve this PDE have used finite difference schemes. Here, I use the method of characteristics to reduce the PDE to a set of ordinary differential equations, which are easy to solve. The method is applied to leaky-integrate-and-fire neurons and produces an algorithm that is efficient and yields a stable and manifestly nonnegative density. Contrary to algorithms based directly on finite difference schemes, this algorithm is insensitive to large density gradients, which may occur during evolution of the density.
Copyright © 2003 by The MIT Press.