This directory contains source code used in Fig. 1 of the publication 

Michalikova M, Remme MWH, Kempter R (2016):
Spikelets in Pyramidal Neurons: Action Potentials Initiated in the Axon 
Initial Segment that do not Activate the Soma, PLOS Comp Biol

This model has been adapted from Hu et al.(2009): "Distinct Contributions 
of Na(V)1.6 and Na(V)1.2 in Action Potential Initiation and Backpropagation" 
and simulates a layer V pyramidal neuron stimulated with fluctuating synaptic 
input that generates somatic APs and spikelets.

To run this model, either auto-launch from ModelDB or first compile
the .mod files (e.g. with "nrnivmodl"), then run the script
"Spikelets_Main.hoc" in GUI ("nrngui Spikelets_Main.hoc").  Press
"Init & Run" to start the simulation. After about 12 minutes (on a
2012 MacBook pro) you should see a graph like:

screenshot

The graph shows voltage traces at the soma (black) and at the AIS
(red). At the soma, APs and spikelets are generated. Somatic spikelets
correspond to APs at the AIS.

Changelog
---------
2022-12: Fix 9.0.0 Upcoming error: new_seed used as both variable and function in file Gfluct.mod