library(trapezoid)
library(readr)

a=2
b=2.1
c=2.25
d=2.5
m=4
n=4
alpha=1


s=10000
GT_Sample<-rtrapezoid(s,a,b,c,d,m,n,alpha)
x_points<-c(0:1000)/1000
x_points<-(d-a)*x_points+a
pdf_points<-dtrapezoid(x_points,a,b,c,d,m,n,alpha)

x_points
pdf_points

#help(formatC)

main_title<-substitute("Histogram: Sample size ="~a,list(a=formatC(s,format="f",big.mark=",",digits=0)))

hist(GT_Sample,breaks=50,col="firebrick3",freq=FALSE,main=main_title)
lines(x_points,pdf_points,lwd=2,col="blue")

GT_Sample<-sort(GT_Sample)

PTime<-data.frame(GT_Sample)
colnames(PTime)<-c("PTime_R")

write_csv(PTime,"Table1_R.csv")