function [absoluteArrivalTimes]=OrderArrivalTimes(T) muOrderTime=1.2; interArrivalTimes(1)=exprnd(muOrderTime); absoluteArrivalTimes(1)=interArrivalTimes(1); n=1; while absoluteArrivalTimes(n)<(T-1) interArrivalTimes(n+1)=exprnd(muOrderTime); absoluteArrivalTimes(n+1)=absoluteArrivalTimes(n)+interArrivalTimes(n+1); n=n+1; end fprintf('There were %0.0f orders placed.\n',size(interArrivalTimes,2)) interArrivalTimes=round(100.*interArrivalTimes)/100; hist(interArrivalTimes,20) xlabel('Time Between Order Arrivals') ylabel('Frequency of Interarrival Time') grid on figure plot(0:(size(interArrivalTimes,2)-1),absoluteArrivalTimes) xlabel('Order Number') ylabel('Time of Order Arrival') grid on