s2rout5.m

Industrial Mathematics

% Routine 2.5
N = 5000;                % number of passerbys
p = 100/N;               % prob that is a potential customer
d = 365;                 % experiment for a year
papers = 90:105;         % experiment: buy 90 to 105 papers
prof = 90:105;           % size the profit vector
prof = 0*prof;           % initialize to 0
sold = prof;             % size # sold vector
                         %
for i=1:d                % start the year
 sold = 0*sold;          % none sold at start of day
  for k=1:N              % people begin walking by
   if rand <= p          % if he/she is a customer,    
    sold = sold+(sold<=papers); % sell a paper if any are left
    end;                % endif   
   end;                 % day ends
 prof=prof+.75*sold-.5*papers;  % add in day's profit
end;                    % end year
prof = prof/d;          % average daily profit 
plot(papers,prof)       % plot average daily profit