clear all
m=16;
for I=1:m
v=rand(1);
Theta=0;
Muestra=(2*v-1)*2;
phi1= 1/(pi*(1+Muestra^2));
Muestras=(2*(1-v)-1)*2;
phi2= 1/(pi*(1+Muestras^2));
Theta=(Theta+phi1+phi2)/2*m;
end
Theta
