    int n_step = 1000000; // Number of time steps //used  10000000
    dt = 0.0001; // Time step size used 0.001

    int n_part = 100; 