mybuffon.m

algotithm to calculate pi

function y=mybuffon(n)
count=0;
for k=1:n
    a=2*rand(1);
    b=2*rand(1);
    if a>=1&&b<=1 count=count+1;
    elseif a>=1&&b>1;
    elseif a<1&&b>=1 count=count+1;
    else;
    end
end
vpa(2*n/count)