%=====================================================
%
% MOS21 - 2012/2013
% 
% Monte-Carlo Methods
% Calcul approche du nombre pi
%
%=====================================================

close all;
clear all;
clc;

n=10000;

tab_x=zeros(n,1);
tab_y=zeros(n,1);
tab_isinside=zeros(n,1);
compteur=0;

for i=1:n
	tab_x(i,1)=rand();
	tab_y(i,1)=rand();
	xi=tab_x(i,1);
	yi=tab_y(i,1);
	r=sqrt((xi-0.5)^2+(yi-0.5)^2);
	if(r<=0.5)
		compteur=compteur+1;
		tab_isinside(i,1)=1;
	end
end

fprintf('* il y a %d / %d points dedans\n',compteur,n);
pi_MC=4.*compteur/n;
fprintf('* pi voisin de %f (erreur %f)\n',pi_MC,(pi_MC-pi)/pi);

figure(1)
hold on
for i=1:n
	if(tab_isinside(i,1)==1)
		couleur=[0,0,1];%RGB, red green blue
	else
		couleur=[1,0,0];
	end
	plot(tab_x(i,1),tab_y(i,1),'marker','o','color',couleur);
end
axis equal;
hold off