function [losses1 losses2]=mc_study(regulator_parameters,system,n)

losses1=zeros(1,n);
losses2=zeros(1,n);
[H u]=load_H;
for i=1:n
    system.b=10*rand()-5;
    apriori=
    while (system.b==0)
        system.b=10*rand()-5;
    end
    losses1(i)=ce(regulator_parameters,system,apriori);
    losses2(i)=rizeni(H,u,system,apriori);
    
end
end