fit_skew_normal function

Fit a skew normal distribution to log-density evaluations