### density plots for meltpt of beeswax 

bees <- matrix(scan("H:\\hw3\\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")