Problem 3.5.17

Item<-c("Oats","Peas","Beans","Barley")
cost<-c(1.05,3.17,1.99,.95)
sell<-c(1.29,3.76,2.23,1.65)
maxD<-c(10,8,14,11)
data<-data.frame(Item,cost,sell,maxD)

trueD_90<-c()
TotalC<-c()
TotalP<-c()
TotalR<-c()

for(i in 1:10000) {

trueD_90<-c(sum(round(runif(90,0,data$maxD[1]),0)),sum(round(runif(90,0,data$maxD[2]),0)),sum(round(runif(90,0,data$maxD[3]),0)),sum(round(runif(90,0,data$maxD[4]),0)))

TotalC[i]<-sum(data$cost*trueD_90)
TotalR[i]<-sum(data$sell*trueD_90)
TotalP[i]<-sum(TotalR[i]-TotalC[i])

}

##Histogram of Costs
summary(TotalC)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    2933    3262    3335    3336    3409    3791

hist(TotalC)

#Normal

## Revenue
summary(TotalR)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    3655    4065    4153    4154    4243    4710

hist(TotalR)

## Profit

summary(TotalP)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   722.0   799.7   817.9   817.9   836.0   926.3

hist(TotalP)