#include <CLHEP/Matrix/Matrix.h>
#include <CLHEP/Matrix/Vector.h>

using namespace std;

#include "chisqMin.h"
#include "parabolaResidFunc.h"

#include "globalGslFuncs.h" // include C-style global functions

int main() 
{
  unsigned int n = 7;
  parabolaResidFunc resid(n);
  residFuncPtr = &resid; // set the global function pointer
  chisqMin fitter(&resid);
  HepVector start(3);
  start(1) = 1.0; start(2) = 0.0; start(3) = 0.0;
  fitter.setStartParams(start);
  fitter.fit();
}