Double_t tf(Double_t *x, Double_t *par){
  return x[0]*par[0]+par[1];
}