#include <TMath.h> 
#include <Math/Polynomial.h> 
#include <complex> 
#include <vector> 
#include <iostream>

void testQuartic() { 
   ROOT::Math::Polynomial poly(4); 
   //cout<<Gamma[i]<<endl; 
   double coeffs[5] = {1,4,6,4,1}; 
   poly.SetParameters(coeffs); 
   vector<complex<double> > roots; 
   roots = poly.FindNumRoots(); 
   std::cout << "number of roots = " << roots.size() << std::endl;
   for (unsigned int j = 0; j < roots.size(); ++j ) 
      std::cout << roots[j] << std::endl;

}