Numerical stability of firing rate derivatives

The implementation of the error function integration in the firing rates used to be numerically unstable. This was improved by replacing `_Phi` by `_Phi_neg` and `_Phi_pos`. But the function `_derivative_of_firing_rates_wrt_mean_input` is still using the old version. This probably could be solved in the same way as for the firing rates.