Matlab implementation to compute the limit cycle trajectory and infinitesimal phase response curve for an oscillating adaptive exponential integrate-and-fire neuron, as well as the interaction function between two identical neurons coupled through a conductance-based synapse model. The phase response curve is calculated using the adjoint method described in: Ladenbauer, Augustin, Shiau, Obermayer, Impact of adaptation currents on synchronization of coupled exponential integrate-and-fire neurons, PLOS Computational Biology 2012
For additional information and helpful comments see also: Ladenbauer, The Collective Dynamics of Adaptive Neurons: Insights from Single Cell and Network Models, PhD thesis (chapter 3.1)