void init_mass_FDF (double X, void * sim_params,  double * f, double * df) {
	Strategy * p = ( Strategy * ) sim_params;
	*f = p->TotalMass(X) - p->offspring_mass();
	*df = p->dTotalMass_dm(X) * 1.0;
}
double init_mass_dF(double X, void *sim_params) {
	Strategy * p = ( Strategy * ) sim_params;
	return p->dTotalMass_dm(X) * 1.0;
}