nrun  =  
      5;       
distance  =  
      Array[d, nrun];
llimit   =  
      1/(nrun + 
            1);   (*the interval endpts*) 
ulimit  =  1  -  (1/(nrun + 1));
inc  =  
      1/(100);    
good  =  
      1;

speed  =  Array[v, nrun];
v[1]  =  1;
v[2]  =  2;
v[3]  =  3;
v[4]  =  4;
v[5]  =  
        5;

 (*as the only time these speeds work is 1/5,  1/6,  
      which do not have finite decimal expansion,  
      the code will produce no data for this set of speeds*) 


For[time =  0,  time  < 1,  time  +=   .001 ,  
      For[j = 1,  j  ≤  nrun,  j++,  
        d[j]  =  Mod[v[j]*time, 1]];  
      For[k = 1;  good  =  1,  k  ≤  nrun,  k++, 
        good  =  
          good  &&  (d[k]  ≥  
                llimit)  &&  (d[k]  ≤  ulimit)];  
      If[good == 1,  Print[time,  Array[d, nrun]],  0]]