Rcpp::NumericVector _a(a), _b(b), _N(N), _x(x), _y=rcauchy(1, _x[0], pow(_a[0], -0.5)), _q=runif(1, 0., 1.); RNGScope scope; // use the same RNG state that R uses double x2=pow(_x[0], 2.), y2=pow(_y[0], 2.), p=pow((1.+y2)/(1.+x2), 0.5*_N[0]-1.5) * exp(_b[0]*(_y[0]*sqrt(1.+y2)-_x[0]*sqrt(1.+x2))-0.5*_a[0]*(y2-x2)); if (_q[0]