phase reduction was wrong
reducing angle to -pi, pi and uneven number of bins
RungeKutta, multiplication by timestep was missing