\(\(f2[energy_, stepsize_, steps_, scaling_] := 
        Module[
{infile, points, p1, p2, \ 
            failures}, 

If[
            scaling > 0, 
failures\  = \ 
              RunThrough[
                StringJoin["$LABBIN/with-scaling ", 
                  ToString[energy], " ", ToString[stepsize], " ", 
                  ToString[
                    steps], " $LABTMP/2-with"], ""]; 
points = 
              ReadList["$LABTMP/2-with", {Number, 
                  Number}];
, 
failures = 
RunThrough[
                StringJoin["$LABBIN/no-scaling ", 
                  ToString[energy], " ", ToString[stepsize], " ", 
                  ToString[
                    steps], \ " $LABTMP/2-no"], ""]; 
points = 
              ReadList["$LABTMP/2-no", {Number, 
                  Number}];
]; 

size = \(Dimensions[
                points]\)[\([1]\)]; 
If[
            failures ≠ 0, 
            Print["Rescaling fails ", 
              failures, " times."]]; 
p1\  = \ 
            points[\([Range[1, Floor[size/2]]]\)]; 
          p2\  = \ points[\([Range[Floor[size/2] + 1, 
                  size]]\)]; 

Show[
            ListPlot[p1, \ PlotStyle -> \ RGBColor[0, 0, 1], 
              DisplayFunction -> Identity], 
            ListPlot[p2, \ PlotStyle -> \ RGBColor[1, 0, 0], 
              DisplayFunction -> Identity], 
            DisplayFunction -> $DisplayFunction];
];\)\)

energy = 1. /18. ;
stepsize = 1. /16. ;
steps = 500000;
scaling = 0;    

f2[energy, 
      stepsize, steps, scaling];

energy  =  1. /18 ;
stepsize = 1. /16. ;
steps = 500000;
scaling = 1;   

f2[energy, stepsize, steps, scaling];

energy  =  1. /8. ;
stepsize = 1. /16. ;
steps = 500000;
scaling = 0;    

f2[energy, stepsize, steps, scaling];

energy  =  1. /8. ;
stepsize = 1. /16. ;
steps = 500000;
scaling = 1;   

f2[energy, stepsize, steps, scaling];