function H0=add_Hc_tilde0(H0,apriori,regulator_parameters)
H0= H0+(apriori.Hmin + (apriori.Hmax-apriori.Hmin)*rand).*(1-regulator_parameters.H_tilde_on);
end