S = 1000000;
x = cauchyrnd(0,1,S);
y = x>2;
thetahat = mean(y)
sy = sqrt(var(y))





