[NEST Users] How to implement double-exponential functions for synapses in C++

Dear NEST Community, is there a good way to implement a double-exponential function for synapses in C++ within a nest module? [cid:image001.jpg@01D87B3E.2D9BE590] As I need to implement my own neuron model in C++, I adapted the code from extraction the information from other neuron models (for example: iaf _cond_beta.cpp) and used the following lines f[S::DG_EXC] = -y[ S::DG_EXC ] / node.P_.tau_synAMPA_on; f[S::G_EXC ] = y[ S::DG_EXC ] - y[ S::G_EXC ] / node.P_.tau_synAMPA_off At first, it seemed to work pretty good. But when I swap the values for rise und decay time, the evolution of g_ADP did not change at all (for example: rise_time = 1ms, decay_time = 100ms). I would be happy to get any feedback on it. Thanks, Pira