#include <stdio.h>
#include <math.h>
#include <TRandom3.h>
#include <time.h>
#include <Math/PdfFuncMathCore.h>
#include <TCanvas.h>
#include <TH1D.h>
#include <TF1.h>

double Binomial(double k, double mu, double n) {
	double a = (tgamma(n+k)/(tgamma(n+1)*tgamma(k)))*(pow(mu/k, n)/pow(mu/(k+1), n+k));
return a;
}

int main() {
	gRandom->SetSeed(0);

	double A = 0.0;
	printf("Gamma function of...? ");
	scanf("%lf", &A);
	A = tgamma(A);	
	printf("\n%.4lf\n", A);

	printf("Binomial PDF evaluated at k = 2, mu = 2.53 =\n");
	A = Binomial(2, 2.53, 3);
	printf("%lf\n", A);
	A = Binomial(2, 2.53, 5);
	printf("%lf\n", A);
	A = Binomial(2, 2.53, 4);

	printf("%lf\n", A);
	//TCanvas *c1 = new TCanvas();
	//TF1 *f = new TF1("f", "ROOT::Math::negative_binomial_pdf(5,0.5,x)", 0, 10);
	//f->Draw();
	//c1->SaveAs("/home/jmerges/public_html/tests/NBDTest.png");

}