#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <fstream>
#include <string>
#include <vector>
#include <sstream>
#include <cmath>
//#include "/home/irfulx157/mnt/fbelloni/workspace/library/lib.cpp"
#include "TColor.h"
#include "TStyle.h"
#include <TFile.h>
#include <TCanvas.h>
#include <TH2F.h>
#include <TSystem.h>
#include "TPaletteAxis.h"
#include "TPaveStats.h"
#include "TMatrixD.h"
#include "TMatrixDSym.h"
#include "TMinuit.h"
#include "TFitter.h"
#include "TVirtualFitter.h"
#include "TROOT.h"
#include "TF2.h"
using namespace std;

void run()
{
     gStyle->SetOptDate(0);
     gROOT->SetStyle("Plain");
     gStyle->SetPalette(1,0);
     TF2*tr = new TF2("Prova","0.5*[0]*( TMath::Erf((x-[1]+[2])/(sqrt(2)*[3]))- TMath::Erf((x-[1]-[2])/(sqrt(2)*[3]))) *( TMath::Erf((y-[4]+[5])/(sqrt(2)*[6])) - TMath::Erf((y-[4]-[5])/(sqrt(2)*[6]))) ", -4,4,-4,4);
     tr->SetParameters(0,2.00489e+04);
     tr->SetParameters(1,-0.0866694);//2*-4.33347e-02
     tr->SetParameters(2,2.43426e-01);
     tr->SetParameters(3,6.53876e-01);
     tr->SetParameters(4,-0.202782);//2*-1.01391e-01
     tr->SetParameters(5,3.13602e-01);
     tr->SetParameters(6,7.04909e-01);
     tr->Draw("surf");
}