library(trapezoid)
library(readr)

png("Example_Historgrams100.png",height=6,width=9,units='in',res=1200)

par(mfrow=c(3,2))
par(mar=c(4,4,2,1),oma=c(0,0,2,0))

par(mfrow=c(2,3))
hist_breaks<-20
s<-10000

#Triangular Sample
a=2
b=4
c=b
d=5
#c=d
m=2
n=2
alpha=1
Sample<-rtrapezoid(s,a,b,c,d,m,n,alpha)

main_title<-substitute("Histogram: Sample size ="~a,list(a=formatC(s,format="f",big.mark=",",digits=0)))
hist(Sample,breaks=hist_breaks,col="firebrick3",freq=FALSE,main=main_title)
#text(2,0.5,"Triangular?",col="deepskyblue4",cex=1.5,pos=4)

#Exponential Sample
lambda=3
Sample<-rexp(s,rate=lambda)

main_title<-substitute("Histogram: Sample size ="~a,list(a=formatC(s,format="f",big.mark=",",digits=0)))
hist(Sample,breaks=hist_breaks,col="firebrick3",freq=FALSE,main=main_title)
#text(1,2,"Exponential?",col="deepskyblue4",cex=1.5,pos=4)


#TSP Sample
a=2
b=2.5
c=b
d=5
#c=d
m=5
n=5
alpha=1
Sample<-rtrapezoid(s,a,b,c,d,m,n,alpha)

main_title<-substitute("Histogram: Sample size ="~a,list(a=formatC(s,format="f",big.mark=",",digits=0)))
hist(Sample,breaks=hist_breaks,col="firebrick3",freq=FALSE,main=main_title)
#text(2.5,1.25,"Asym. Laplace?",col="deepskyblue4",cex=1.5,pos=4)


#Gamma Sample
alpha=3
beta=4
Sample<-rgamma(s,shape=beta,rate=alpha)

main_title<-substitute("Histogram: Sample size ="~a,list(a=formatC(s,format="f",big.mark=",",digits=0)))
hist(Sample,breaks=hist_breaks,col="firebrick3",freq=FALSE,main=main_title)
#text(2.5,0.5,"Gamma?",col="deepskyblue4",cex=1.5,pos=4)

# Weibull Sample
alpha=2
beta=2
Sample<-rweibull(s,shape=beta,scale=alpha)

main_title<-substitute("Histogram: Sample size ="~a,list(a=formatC(s,format="f",big.mark=",",digits=0)))
hist(Sample,breaks=hist_breaks,col="firebrick3",freq=FALSE,main=main_title)
#text(2.5,0.35,"Weibull?",col="deepskyblue4",cex=1.5,pos=4)

# beta Sample

beta<-1.8
alpha<-0.5
a<-alpha*beta
b<-(1-alpha)*beta
Sample<-rbeta(s,shape1=a,shape2=b)

main_title<-substitute("Histogram: Sample size ="~a,list(a=formatC(s,format="f",big.mark=",",digits=0)))
hist(Sample,breaks=hist_breaks,col="firebrick3",freq=FALSE,main=main_title)
#text(0.5,1.0,"Uniform?",col="deepskyblue4",cex=1.5,pos=3)
dev.off()