#include <iostream.h>
#include <fstream.h>
#include <math.h>
#include <cfortran/cfortran.h>
#include <cfortran/hbook.h>
void get_data(int*, double*, double*, double*);
const int MAX_BINS = 100;

void main(){

  int numBins;
  double n[MAX_BINS], nu_theory1[MAX_BINS], nu_theory2[MAX_BINS];
  get_data(&numBins, n, nu_theory1, nu_theory2);
  for (int i=0; i<numBins; i++){
    cout << n[i] << "  " << nu_theory1[i] << "  " << nu_theory2[i] << endl;
  }  

}

void get_data (int* numBins, double* n, double* nu_theory1, 
               double* nu_theory2) 
{

  int num;
  double xmin[MAX_BINS], xmax[MAX_BINS];

  ifstream dataIn ("data_1.dat");
  dataIn >> num;
  *numBins = num;
  for(int i=0; i<num; i++)
  {
    dataIn >> xmin[i] >> xmax[i] >> n[i];
  }
  dataIn.close();

  ifstream theory1In ("theory_1.dat");
  theory1In >> num;
  *numBins = num;
  for(int i=0; i<num; i++)
  {
    theory1In >> xmin[i] >> xmax[i] >> nu_theory1[i];
  }
  theory1In.close();

  ifstream theory2In ("theory_2.dat");
  theory2In >> num;
  *numBins = num;
  for(int i=0; i<num; i++)
  {
    theory2In >> xmin[i] >> xmax[i] >> nu_theory2[i];
  }
  theory2In.close();

}