x1 = [] #list of 1st iteration coordinates for i in range (num_discs): x_coords1 = xcoords[i] + (velocities[i]*delta_t) + (0.5*((list_of_f[i])/m)*(delta_t**2))