# filename: coupled_oscillators.ode
# Kuramoto model
#

p K=0

table w % 10 0  9  1+((t-49.5)/49.5)^3
table kk % 11 0 10 1
gam(u,v)=sin(u-v)
special r=fconv(periodic,10,5,kk,u0,u0,gam)


u[0..9]'=w([j]) + K/10*r([j])

aux rr[0..9]=r([j])


@ autoeval=0, meth=cvode,nplot=5