function losses=mc_study(regulator_parameters,system,apriori,n)
losses=zeros(1,n);
for i=1:n
    system.b=10*rand()-5;
    [H,u,losses(i)]=ce(regulator_parameters,system,apriori);
end
end