#include <TSystem.h>
#include <TRandom3.h>
#include <TCanvas.h>
#include <TF1.h>
#include <math.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() {
	gSystem->Load("libMathCore");
	
	TF1 * NBD = new TF1("NBD", "Binomial(2,2.53,x)", 0,5);
	TCanvas * c1 = new TCanvas("c1", "Match Fit", 800, 800);
	NBD->Draw();
	c1->SaveAs("/home/jmerges/public_html/tests/EnergyProductionFit.png");
return 0;
}