function [spi]=simulace_pi(n)
% funkce aproximuje cislo pi 
% pomoci ctyrnasobku relativni cetnosti uspechu
% uspech nastane, kdyz soucet kvadratu 
% dvou kladnych cisel mensich nez 1 je mensi nez 1
% vstupni parametr: n - pocet pokusu
% vystupni parametr: spi - aproximace cisla pi
u = unifrnd(0,1,n,2);
s = u.^2;
spi = 4 * sum(s(:,1) + s(:,2) < 1) / n;

