#include <iostream>
#include <string>
#include <vector>
#include <TH1D.h>
#include <TRandom.h>
#include <TCanvas.h>
#include <TMath.h>
#include "PLUtil.h"

using namespace std;

double PLUtil::BernsteinPoly(double x, int k, int n){
  double kd = static_cast<double>(k);
  double nminuskd = static_cast<double>(n-k);
  double f = TMath::Binomial(n,k)*pow(x,k)*pow(1.-x,nminuskd);
  return f;
}