### density plots for meltpt of beeswax 

bees = matrix(scan("C:\\Users\\Eric A. Suess\\Desktop\\beeswax.dat"),ncol=2,byrow=T)
bees

meltpt = bees[,1]
phydcarb = bees[,2]

par(mfrow=c(2,1))

hist(meltpt)

# density plot for meltpt
# n is the number of equally spaced points at which to estimate the density
# window is the number of equally spaced points at which to estimate the density

plot(density(meltpt, n=50, window="g"),type="b",xlab="Melting Point C",ylab="density")