#include "Rtypes.h"
#include "TH1.h"
#include "TMath.h"
#include "TRandom.h"
#include <iostream>
using namespace std ;

double func(double x) 
{
  return TMath::Gaus(x,0,3,kTRUE) ;
}


void ex2()
{
  TH1* h = new TH1D("h","h",1000,-10,10) ;

  Double_t mean = 0 ;
  double fmax = func(mean) ;
}