#include <stdio.h>
#include <stdlib.h>
#define NDIM 100

void macroWriteArraysInFile_ofstream()
{
  
  TRandom3 R(32525);

  const char* filename="nuovidati.txt";
  ofstream myOutFile;
  
  myOutFile.open(filename);

  int Nlines=0;
  double a1, a2, a3, a4;

  a1 = 0.;
  double astep = 1.;
  double errY  = 1.;
  double errX = 0.1;
  while (Nlines <10) 
    {
      a1 += astep+errX*(R.Uniform()-0.5);
      a2 = a1 * 2.56 + a1 + 0.45;
      a2 = a2+R.Gaus(0.,errY);
      a3 = errX;
      a4 = errY;
      cout<<"Writing to file  x = "<<a1<<" y = "<<a2<<" ex = "<<a3<<" ey = "<<a4<<endl;
      myOutFile<<a1<<"  "<<a2<<"  "<<a3<<"  "<<a4<<endl;
      ++Nlines;
    }
  if (Nlines>=NDIM) cout<<"Maximum size of data arrays reached; N = "<<Nlines<<" lines of data available"<<endl;
  cout<<"N = "<<Nlines<<" lines of data are written to file <"<< filename<<">"<<endl;
  myOutFile.close();

  /// now you can manipulate the file 
  
}