clear all
s=20;
a=13;
b=1;
m=16;
for I=1:m
       n=a*s+b;
u=mod(n,m);
v=u/m;
Theta=0;
Muestra=(2*v-1)*2;
phi= 1/(pi*(1+Muestra^2));
Theta=Theta+phi/m;
s=u;
end
Theta
Var=0.114035/m
