Listing 1: Using the toolbox.

CSolverSettings settings;
// ... set member variables of settings as appropriate...
CSolver* pSolver = 
        CSolverFactory::CreateSolver(CSolverFactory::MILSTEIN, settings);
if (NULL != pSolver)
{
    CDynSystem* pSystem = new CGBMDynSystem();
    pSolver->RegisterDynSystem(pSystem);
    if (!pSolver->Run())
    {
        fprintf(stderr, "%s\n", pSolver->GetLastError().c_str());
        exit(1);
  }
}