s2rout1.m

Industrial Mathematics

% Routine 2.1
A = 1;           % area of enclosing rectangle 
N = 10000;       % set the number of trials
s = 0;           % initialize the success counter
for  i = 1:N     % begin trials
 x = rand;       % choose a random  x-coordinate
 y = rand;       %  choose a random y-coordinate
 if y  <=  exp(-x^3)  % if below the curve, then
    s = s + 1;   % increment successes
    end;         % end if 
end;             % end of each trial
I = A*s/N        % integral= area*success/trials