1) trap the case where the model does not solve for any parameter values in the MC sample;
2) fixed initialization of alpha2_stab;