#include <iostream>
#include <stdio.h>
#include "TFile.h"
#include "TRandom.h"
#include "TMath.h"
#include "TF1.h"
#include "TH1D.h"
#include <string>
#include <cmath>
#include "gsl/gsl_sf_ellint.h"
#include "Math/WrappedTF1.h"
#include "Math/GaussIntegrator.h"

void Code()
{
 float V;
 TF1 f("f", "TMath::Exp(-1/x)*(1/x^2)*(TMath::comp_ellint_1(x)-TMath::comp_ellint_2(x)/(1-x^2))", 0,0.99);
 ROOT::Math::WrappedTF1 wf(f);
 ROOT::Math::GaussIntegrator ig;
 ig.SetFunction(wf);
 ig.SetRelTolerance(0.001);
 V=ig.Integral(0,0.99);
}