randcorr.sample.sink function

Sample from the (unnormalized) distribution sin(x)^k, 0 < x < pi, k >= 1