def METHOD_NAME(frame_number):
    updated_artists = []
    for (dots, arrow, mean_dot), K, Theta in zip(artists, K_values, theta_values):
        t = frame_delay*frame_number
        ts = timestep(t, defaultclock.dt)
        x, y = to_x_y(Theta.T[ts])
        dots.set_offsets(np.vstack([x, y]).T)
        r, phi = calc_coherence_and_phase(x, y)
        arrow.set_data(dx=r*np.cos(phi), dy=r*np.sin(phi))
        mean_dot.set_data(r*np.cos(phi), r*np.sin(phi))
        updated_artists.extend([dots, arrow, mean_dot])
    return updated_artists