MAR function

C++ function to estimate Pitman-Yor univariate mixtures via marginal sampler - LOCATION SCALE