dθᵢ/dt = ωᵢ + (K/N)Σⱼsin(θⱼ−θᵢ) — order parameter r = |Σe^(iθⱼ)|/N → 1 at synchronization (K > K_c = 2/πg(0))