function vypocetpi

figure;axes;
box on;
rectangle('Position',[0,0,1,1],'Curvature',[1,1]);
hold on;
kapek=0;    % pocet kapek
vkruhu=0;   % pocet kapek v kruhu
npi=0;      % odhad pi
steps=300  % pocet kroku
plotNum=zeros(steps,1);
plotPi=zeros(steps,1);

for i=1:steps
        str=sprintf('Pocet kapek %u Vypoctene PI je %g',kapek,npi);
        h=title(str);
  %      waitforbuttonpress;
        for j=1:50,
                kapek=kapek+1; 
                x=rand(1,2);    % dve nahodna cisla - bod (x1,x2)
                % jsou uvnitr kruhu se stredem (0.5,0.5) a polomerem 0.5?
                if ((x(1)-0.5)^2+(x(2)-0.5)^2<0.25) 
                        vkruhu=vkruhu+1;    % pricteme je
   %                     plot(x(1),x(2),'o','Color',[1 0 0]);
                else
     %                   plot(x(1),x(2),'o','Color',[0 0 1]);
                end;
        end;
        npi=4*vkruhu/kapek; 
%        delete(h);
        plotNum(i)=kapek;
        plotPi(i)=npi;
	i
end;

%str=sprintf('Pocet kapek %u Vypoctene PI je %g',kapek,npi);
%title(str);



figure
plot(plotNum,plotPi, [0 15000], [pi pi]);
axis([0 15000 3 3.3]);

end
