Function for reverse computing mean squared diffusion
function for reproducing spiking behavior of cortical neurons from count data of synaptic surface receptors

The postspike() function takes count data of synaptic surface receptors
over time. Using count data from these discrete timepoints, the
function will inject a standard upstream input and compute the
postmembrane potential. This function is adapted from a model by E.
Izhikevich that reproduces spiking and bursting behavior of known types
of cortical neurons. The model combines the biologically plausibility
of Hodgkin-Huxley-type dynamics and the computational efficiency of
integrate-and-fire neurons.