int application::main(kit* tk,int argc,char* argv[])
  {
    simulation* sim = new simulation();
    histogram* h1 = new histogram(".h1","-title curve");
    histogram* h2 = new histogram(".h2","-title line");
    analysis* a1 = new analysis(h1);
    analysis* a2 = new analysis(h2);
    for (int i=0;i<10;i++)
    {
      double val1 = 1.0*(i+0.5)*(i+0.5) - 2.5*(i+0.5) + 19.1;
      double val2 = 8.5*(i+0.5) + 9.1;
      h1 -> sample((i+0.5),val1);    // sample first curve
      h2 -> sample((i+0.5),val2);    // sample second curve
    }
    a1 -> quadratic();            // fit curves
    a1 -> linear();
    a2 -> quadratic();
    a2 -> linear();
    cout << (*h1);
    cout << (*h2);
    delete a1;
    delete a2;
    delete h1;
    delete h2;
    delete sim;
    return 0;
  }