botev function

Compute the optimal kernel bandwidth