(* Needs["Histograms`"] *)

(* Testa o Croce *)

Quit[];

normal[μ_,σ_,x_]=PDF[NormalDistribution[μ,σ],x];
Fi[t_]:=Integrate[normal[0,1,x],{x,-Infinity,t}];

n=10;
β=0.9;
α=1.-β;
γ=1.-0.5*α;
sol1=FindRoot[Fi[t]==γ,{t,0.5}][[1]];
∅=t/.sol1;
fel[x_,y_]:=y^2*(n+∅^2)-y*(2*n*x+∅^2)+n*x^2
Do[ix=RandomInteger[];irand[i]=ix,{i,1,n}]
randtable=Table[irand[i],{i,1,n}];
ListPlot[randtable,PlotRange->{-0.5,1.5},Axes->{False,True}]

m=1.*Mean[randtable]
1.*Sqrt[Variance[randtable]];
σ=1.*Sqrt[p*(1-p)];
p9010=ContourPlot[fel[x,y],{x,-0.2,1.2},{y,-0.2,1.2},Contours->{0.},ContourShading->None,ContourStyle->{Thickness[0.01],RGBColor[0,0,0]},Axes->True]

p10=Show[p9010,p9510]
p20=Show[p9020,p9520]
p100=Show[p90100,p95100]

p90=Show[p9010,p9020,p90100]
p95=Show[p9510,p9520,p95100]