integrand=function(xeval,x){

#fuction to be evaluated for integration..

# For example the Gaussian distribution..
#parameters of the function
#zsd=sd(x)
#zvar=var(x)
#zmean=mean(x)
#fact=1/(sqrt(2*pi)*zsd)
#xdens=fact*exp(-(xeval-zmean)*(xeval-zmean)/(2*zvar))

#OR

xdens = dnorm(xeval, mean=mean(x), sd=sd(x))

xdens

}