clear all
m=16;
Thever=0.1476;
for I=1:m
v=rand(1);
Theta=0;
Mu=pi*(v-0.5);
Muestra=tan(Mu);
if Muestra >= 2
    phi=1;
else
    phi=0;
end
u=rand(1);
if u <= 1/2
    psi=1;
else 
    psi=0;
end
Theta=Theta+phi/m+0.2956*(psi-1/2)/m;
end
Theta
