void test1(){
  float ScaleFactor = 3.;

  TVirtualFitter::SetDefaultFitter("Minuit2");

  const int N = 5;
  float x[N] = { 40., 50., 55.,  60., 70.};
  float y[N] = { 67., 69., 69.8, 70., 74.};
  
  float *dy;
  float dyU[N] = { 2.0, 1.9, 1.5, 1.7, 2.0};
  float dyL[N] = { 1.5, 1.4, 1.5, 1.2, 1.5};
  
  for(int e=0; e<N; e++){
    dyU[e] *= ScaleFactor;
    dyL[e] *= ScaleFactor;
  }
  
  TH2F *h1 = new TH2F("h1","",1000,30.,80.,1000,40.,90.);
  TGraphAsymmErrors *g1 = new TGraphAsymmErrors(N,x,y,0,0,dyL,dyU);
  TF1 *f1 = new TF1("f1","[0]+[1]*TMath::Log10(x)",35.,75.);
  
  TCanvas *c1 = new TCanvas("c1","CANVAS 1",500,400);
  c1->cd();
  h1->Draw();
  g1->SetMarkerStyle(20);
  g1->Draw("Psame");
  g1->Fit("f1");
  gPad->Update();

  TGraphErrors *gci = new TGraphErrors(N);
  gci->SetFillColor(5);
  for(int i=0; i<N; i++)
    gci->SetPoint(i,g1->GetX()[i],0);
  (TVirtualFitter::GetFitter())->GetConfidenceIntervals(gci,0.95);
  c1->cd();
  gci->Draw("e3 same");
  g1->Draw("p");

  double xx[N] = { 40., 50., 55., 60., 70.};
  double *ci;
  ci = new double[N];
  (TVirtualFitter::GetFitter())->GetConfidenceIntervals(N,1,&xx[0], &ci[0], 0.95);
  printf("\n");
  for(int e=0; e<N; e++){
    printf("%2d --> %8.5E\n",e+1,ci[e]);
  }
}