Can anyone see what is wrong with this implementation of the Pinsky-Rinzel model in Brian2?
If anyone ever finds this and has the same problem: The original paper has a flaw in the paper:

alphac=(exp((Vd/mV-10)/11)-exp((Vd/mV-6.5)/27))/(18.975)* (Vd&lt;=50*mV)/ms+2*exp((6.5-Vd/mV)/27)*(Vd&gt;50*mV)/ms : Hz


An errata published later changes this to:

alphac=(exp((Vd/mV-10)/11-(Vd/mV-6.5)/27)/18.975)*(Vd&lt;=50*mV)/ms+2*exp((6.5-Vd/mV)/27)*(Vd&gt;50*mV)/ms : Hz


After making this change it acts exactly how it should.