bw_cv_likelihood_calc function

Bandwidth selection by likelihood cross validation