\[\dot\theta_n = \omega_n +\frac{K}{N}\sum_m \sin(\theta_m-\theta_n). \]

based on Complexity Explorable:
" Ride my Kuramotocycle! -
Synchronization of Phase-Coupled Oscillators
"