
 resolution = 8; // Number of points per decade
    min_rate = 10; // Min number of samples
   
   
   
   int n_step = 10000000; // Number of time steps //used  10000000
    dt = 0.000001; // Time step size used 0.001
int n_part = 100; 

    for f=0
        int n_step = 100000000; // Number of time steps //used  10000000
    dt = 0.00001; // Time step size used 0.001
    int n_part = 100; 