simulation of Kuramoto model