\(\(\(
\)\(
\)\( (*\ 
      The\ gap\ function, \ Epsilon \((t)\), \ 
      is\ read\ from\ the\ output\ data\ file\ *) \)\(
\)\(
\)\(\(gapfuncfile\  = \ 
        ReadList["data.dat", \ {Number, \ Number}];\)

    \(gapfunc\  = \ 
        Table[\(gapfuncfile[\([k]\)]\)[\([2]\)], {k, \ 
            Length[gapfuncfile]}];\)\)\)\)

ListPlot[gapfunc, PlotJoined -> True,  PlotRange -> All,  
        Axes -> False,  Frame -> True,  
        PlotStyle -> { RGBColor[0, 0.2, 0.6]}];

\(\(\(
\)\(
\)\(
\)\( (*\ 
Now\ read\ the\ sample\ output\ from\ a\ longer\ \((10  ps)\)\ simulation\ 
You\ are\ given\ a\ larger\ data\ file\ to\ read\ for\ \(\(this\)\(.\)\)
*) \)\(
\)\(
\)\(\(EpsilonOfTfile\  = \ 
        ReadList["dataLONG.dat", \ {Number, \ 
            Number}];\)

    \(EpsilonOfT = \ 
        Table[\(EpsilonOfTfile[\([k]\)]\)[\([2]\)], {k, \ 
            Length[EpsilonOfTfile]}];\)

    \(noframes\  = \ Length[EpsilonOfTfile];\)\)\)\)

ListPlot[EpsilonOfT, PlotJoined -> True,  
        PlotRange -> All,  Axes -> False,  Frame -> True,  
        PlotStyle -> { RGBColor[0.2, 0. , 0.5]}];

<< Graphics`Graphics`

<< Statistics`DescriptiveStatistics`

Histogram[EpsilonOfT,  Frame -> True];

meanEps  =  
      Sum[EpsilonOfT[[k]], {k, noframes}]/noframes

varianceEps =  
      Sqrt[ (Sum[ (EpsilonOfT[[k]] - meanEps)^2, {k, 
                noframes}])/noframes]