// FILE TwoBody.h

#include "TwoBody.h"


void TwoBody::LoadProblemParam( int InitCond )
{
  SetT0( 0 );
  SetTend( 20 );
  
  INTERVAL_VECTOR Y(4);
  REAL eps = 0.9;
  
  Y(1) = 1.0 - eps;
  Y(2) = 0.0;
  Y(3) = 0.0;
  Y(4) = sqrt( (1.0+eps)/(1-eps) );
  
  SetInitCond(Y);
}