function [] = MC_pi(n)
  x=4*(1-rand(1,n).^2).^(1/2);
  X=cumsum(x)./(1:n);
  var_x=cumsum(x.^2)./(1:n) - X.^2;
  
  y=(4*(rand(1,n).^2+rand(1,n).^2<=1));
  Y=cumsum(y)./(1:n);
  var_y=cumsum(y.^2)./(1:n) - Y.^2;
 	
  z=(3*2*(rand(1,n).^2+rand(1,n).^2+rand(1,n).^2<=1));
  Z=cumsum(z)./(1:n);
  var_z=cumsum(z.^2)./(1:n) - Z.^2;
 
  PI=%pi*ones(1,n);
  
  xset('window',0)
  xbasc(0);
  plot2d([1:n],PI,1)

  plot2d([1:n],X,3)
  plot2d([1:n],X+1.96*sqrt(var_x)./sqrt([1:n]),2) //Intervalle de confiance
  plot2d([1:n],X-1.96*sqrt(var_x)./sqrt([1:n]),2) //à 95%

  xset('window',1)
  xbasc(1);
  plot2d([1:n],PI,1)

  plot2d([1:n],Y,4)
  plot2d([1:n],Y+1.96*sqrt(var_y)./sqrt([1:n]),2) //Intervalle de confiance
  plot2d([1:n],Y-1.96*sqrt(var_y)./sqrt([1:n]),2) //à 95%
  
  xset('window',2)
  xbasc(2);
  plot2d([1:n],PI,1)

  plot2d([1:n],Z,5)
  plot2d([1:n],Z+1.96*sqrt(var_z)./sqrt([1:n]),2) //Intervalle de confiance
  plot2d([1:n],Z-1.96*sqrt(var_z)./sqrt([1:n]),2) //à 95%
  
endfunction