Hi,
Thank you for the interesting question. It depends a bit on what your intended outcome looks like for including two separate decay components. A potential issue with the approach you suggest is that the fast component will come to dominate the slow, and you will not see effects of the latter. Instead, it might be necessary to include a third state variable, so that fast and slow components can be tracked separately. A simple test script that I coaxed out of an LLM provided the following result. Would this work for you? The script is attached.
With kind regards,Charl
> Dear experts,
>
> I am writing to request your feedback about the modification of the
> standard NMDA conductance β‑function to include two separate decay
> components.
> Here our implementation of the related kernel in nestml:
>
> kernel g_NMDA' = g_NMDA$ - NMDA_A1_decay1 * (g_NMDA / NMDA_Tau_decay1)
> - NMDA_A2_decay2 * (g_NMDA /NMDA_Tau_decay2),
> g_NMDA$' = NMDA_A_rise * (-g_NMDA$ / NMDA_Tau_rise)
>
> Do you think that this writing might be the right way to implement it?
>
> Thank you very much for your insight.
>
> Kind regards,
> MG
> _______________________________________________