To calculate the Discs of β, run the first cell and then set α = β. To calculate γ and η, run their respective notebooks and then in this notebook set α = γ or η. They do not need to be calculated in this notebook to have their values stored. They just need to be calculated in that session of Mathematica  and their values will be stored globally.

p = 3;
β  =  (∑_(i = 1)^7) (1/(p^(i!)));

<< Graphics`MultipleListPlot`

K =  50000;
miny  =  (10^15);
Digits  =  60;
α  =  β;
Table1  =  
        Table[FractionalPart[N[(k^2)*α, Digits]], {k, miny, 
            K + miny - 1}];
Table2  =  Sort[Table1];
Table3  =  
        Table[K*(Table2[[p + 1]] - Table2[[p]]), {p, 
            K - 1}];
Table3  =  
        Append[Table3, 
          K*(1 + Table2[[1]] - 
                Table2[[K]])];
Table4  =  Sort[Table3];

Table5  =  
        Table[{Table4[[k]], (k/K)}, {k, K}];
Table6  =  
        Table[{20*(i/K), 1 - Exp[-20*(i/K)]}, {i, 
            (K/2)}];
MultipleListPlot[Table5, Table6, 
      SymbolStyle -> {RGBColor[1, 0, 0], RGBColor[0, 0, 1]}]

Disc  =  0;
For[i = 1, i < K + 1, i++, 
        Disc =  Max[Disc, 
            Abs[(i/K) - (1 - 
                    Exp[-Table4[[i]]])]]];
Disc

Disc2  =  0;
For[i = 1, 
        i < K + 1, i++, 
Disc2  =  
          Disc2  + 
            Abs[(i/K) - (1 - 
                    Exp[-Table4[[i]]])]];
Disc2  =  (Disc2/K)