%% Monte Carlo - tikkataulu

%% ympyrä ja neliö tyylikkäästi
theta = linspace(0,2*pi,256);
phi = linspace(pi/4,2*pi+pi/4,5);
% 'tikkojen' lkm
N = 600000;
% piirtoa
plot(cos(theta),sin(theta),sqrt(2)*cos(phi),sqrt(2)*sin(phi))
grid
axis([-1.5,1.5,-1.5,1.5])

%% arvotaan satunnaispisteet
X = 2*rand(N,1)-1;
Y = 2*rand(N,1)-1;
hold on

% osuneet
in = X.^2+Y.^2<=1;
plot(X(in),Y(in),'bo',X(~in),Y(~in),'ro')

% estimaatti piille
disp(4*(sum(in)/N))
disp(pi)

hold off