/* Header file HW03_03.c */

#define G_univ 6.67428e-20  /* km**2/(kg sec**2) */
#define max_body 1000       /* Maximum number of bodies allowed 
                               (saves pre-reading file */

int num_body ; // Number of bodies

double dur ; // Duration of integration days
double tol ; // Tolerance on integaton fractional

double mass[max_body] ; // Masses (kg)
double posi[2][max_body] ; // Initial XY coordinates (km)
double veli[2][max_body] ; // Initial XY velocities (km/s)
double closest = 1.e20 ;  //  Closest distance any body reached (km)
int close_body[2] ;  // Two closest bodies

char infile[128] ; // Name of input file with mass pos vel of each body