MATLAB: How to do a while loop to estimate pi using monte carlo method

MATLAB Online Serverwhile loopwhileloop

hi guys..sorry very new to matlab. i have been trying to write a while loop code to estimate pi using monte carlo method so it displays number of hits and with a tollerance of about 1e-4. please find attached the code

function [hits]=mypi_tolarance(x,y) 
hits=0;
for n=1:10000
  x=rand;
  y=rand;
  hit=x^2+y^2<=1;
  hits=hits+hit;
  
end
 while 1
      pi_est=hits/10000*4;
     
      if pi_est<3.14||pi_est>3.142
          break
      end
     disp(pi_est) 
     disp(hits)
 end

Best Answer

https://in.mathworks.com/matlabcentral/fileexchange/33543-monte-carlo-method?focused=5201749&tab=function

Related Solutions

MATLAB: Estimating pi using for and while loop for specific tolerance.

function [estpi,n] = mypi(tol)
error=1;
hits=0;
n=0;
while error>tol
   if norm(rand(2,1))<=1
       hits=hits+1;
   end
   n=n+1;
   error=abs(pi-4*hits/n);
end
estpi=4*hits/n;
end

MATLAB: Approximate the value of π using the following algorithm based on geometric probability. This algorithm is an example of the so called monte carlo method

Hi,

The code seems to work fine for computing PI using Hit & Miss Monte Carlo algorithm. In this code the variable 'c' is the number of hits and the variable 's' is the total number of samples i.e. hits + misses.

This code can be optimized further in terms of performance in MATLAB using vectorization. This link mentions how this can be done. However for illustration purpose I have taken the liberty to modify the code as below to show how to compute hits and misses using both vectorization and looping techniques. If this is not what you were looking for then kindly elaborate on the query.

% MC estimate of pi
clear all;
clc;
%%initialize random points
n=10000;
x=rand([1 n]);
y=rand ([1 n]);
%%compute using vectorization
radii = sqrt(x.^2+y.^2);
hits = sum(radii<=1);
misses = n-hits;
pi_mc = 4*(hits/n); 
fprintf('\nUsing Vectorization:: HITS = %d, MISSES = %d, PI = %f\n',hits,misses,pi_mc);
%%compute using a loop
c=0;
s=0;
for i=1:n
    s=s+1;
    if x(i)^2 +y(i)^2 <=1 % inside circle 
        c=c+1;
    end
end
p=c/s;
pi_=p*4;
fprintf('\nUsing Loop:: HITS = %d, MISSES = %d, PI = %f\n',c,s-c,pi_);

Best Answer

Related Solutions

MATLAB: Estimating pi using for and while loop for specific tolerance.

MATLAB: Approximate the value of π using the following algorithm based on geometric probability. This algorithm is an example of the so called monte carlo method

Related Question