cbetdev.cpp

Go to the documentation of this file.

#include <fvar.hpp>

double beta_deviate(double a,double b,double x,double eps)
{
  double y=cumd_norm(x);
  y=.9999999*y+.00000005;
  double z=inv_cumd_beta_stable(a,b,y,eps);
  return z;
}