// FILE DemoDetestC3.cc

#include "VNODE.h"    
#include "DetestC3.h" 


int main()
{
  
  VnodeInit(); 
  
  PtrODENumeric ODE = new DetestC3; 
  PtrDataRepr DataRepr = new TAYLOR_EXPANSION<DetestC3TaylGenODE,DetestC3TaylGenVAR>;
  
  ODE->LoadProblemParam( 1 );
  PtrODESolver Solver = new SOLVER_1( ODE, DataRepr );


  cout << "*** Integrating " << ODE->GetName() << endl;
  Solver -> Integrate();
  
  cout << endl << "\n *** Enclosure at T = " 
       << ODE->GetTcur() << endl;
  PrintVec( ODE->GetTightEncl() );
  
   return 0;
}