Double_t MyFunction(Double_t *x, Double_t *par) {
  Double_t arg = 0;
  if (par[2]) arg = (x[0] - par[1])/par[2];
  return par[0]*TMath::Exp(-0.5*arg*arg);
}