% Program name: random_1.m
% Program to find area under curve by Monte Carlo method
% Consider a quarter of a cirular disk in range of
% 0 < x < 1 and 0 < y < 1. We are interested
% to find its area using Monte Carlod method.
% This generates a movie of 20 frames with
% increasing number of dots thrown at the region.
N = 5;
hit = 0;
mis = 0;
xs = 1;
ys = 0;
rand('state', 0) % this sets the generator to an initial state
nframes = 200; % number of frames for movie
delta = 1.0/N
for i = 1: nframes
i
for j = 1: N
xs = xs - delta;
ys = abs(sqrt(1- (xs*xs)));
plot(xs, ys, '--')
axis([0 1 0 1])
hold on
end
for k = 1: N
x = rand;
y = rand;
plot(x, y, 'r:.')
hold on
end
m(k)= getframe;
end
movie(m, 100) % play movie 30 times.
docsity.com