function dth = Kuramoto_ode(t,th,om,K)


dth = om + K/length(om) * sum(sin(  th'  - th  ),2);