public BetaDistribution(double a, double b)
 {
     alpha  = a;
     beta   = b;
     gamma1 = new GammaDistribution(alpha, 1);
     gamma2 = new GammaDistribution(beta, 1);
 }
 public BetaDistribution()
 {
     alpha  = 1;
     beta   = 1;
     gamma1 = new GammaDistribution(alpha, 1);
     gamma2 = new GammaDistribution(beta, 1);
 }