Density function added to preAnalysis.R need to work on it.