integrand=function(xeval,x){ #fuction to be evaluated for integration.. # For example the LogNormal distribution.. xdens = dlnorm(xeval, meanlog=mean(x), sdlog=sd(x)) xdens }