void fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
{
  int check;
  cout << "in fcn" << endl;

  f=func(par);
  cout << "after func" << endl;                  
  if(call%1000==0) cout<<"chisq in fcn: "<<f;
 
  
  if(call==1) cout<<"\n\n Initial chisq: "<<chisq<<"\n\n";
  outdata.open ("chiS.dat",ios::app);
  assert (!outdata.fail());
  if(call%1000==0) outdata <<f<<endl;
  outdata.close();
  cout << "end of fcn" << endl;
}