\(\(f3[timestep_, \ numsteps_, \ cutoff_, \ distance_, \ method_] := 
        Module[
{dataf, 
            points}, 

dataf = "$LABTMP/out.sim"; 
Run["$LABBIN/argon", timestep, 
            numsteps, cutoff, distance, method, 
            dataf]; 
points = 
            ReadList[
              dataf, {Number, 
                Number}]; 

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

timestep  =  10;
numsteps  =  200000;               
cutoff  =  
        6;                               
switchingDistance  =  6;       
switchingMethod =  0;           
 (* 
      0,  1,  2,  
      or 3 : 
 
        0 discontinuous potential, 
  1 continuous potential, 
  2 continously differentiable potential, 
  3 twice continously differentiable potential *) 

f3[timestep, numsteps,  cutoff,  switchingDistance,  
        switchingMethod];

timestep  =  10;
numsteps  =  200000;               
cutoff  =  
        8;                               
switchingDistance  =  6;       
switchingMethod =  1; 
 (* 0,  1,  2,  
      or 3 : 
 
        0 discontinuous potential, 
  1 continuous potential, 
  2 continously differentiable potential, 
  3 twice continously differentiable potential *) 

f3[timestep, numsteps,  cutoff,  switchingDistance,  
        switchingMethod];

timestep  =  10;
numsteps  =  200000;               
cutoff  =  
        8;                               
switchingDistance  =  6;       
switchingMethod =  2;
 (* 0,  1,  2,  
      or 3 : 
 
        0 discontinuous potential, 
  1 continuous potential, 
  2 continously differentiable potential, 
  3 twice continously differentiable potential *) 

f3[timestep, numsteps,  cutoff,  switchingDistance,  
        switchingMethod];

timestep  =  10;
numsteps  =  200000;               
cutoff  =  
        8;                               
switchingDistance  =  6;       
switchingMethod =  3; 
 (* 0,  1,  2,  
      or 3 : 
 
        0 discontinuous potential, 
  1 continuous potential, 
  2 continously differentiable potential, 
  3 twice continously differentiable potential *) 

f3[timestep, numsteps,  cutoff,  switchingDistance,  
        switchingMethod];