Rgamma (rgammadouble)

Wraps R.rgamma for correct input and output

   double rgammadouble(int a, double b, double c)
   {   Rcpp::NumericVector x = rgamma(a,b,1/c);
       return x(0);
   }