N = 50
K = 1.0
theta0 = zeros(N, float)
for i in range(N):
    theta0[i] = random.uniform(0, 2 * pi)
odecons = array([K])
odeist = theta0
odestsc = array([0.0])
ode_eg = p.dde()
ode_eg.dde(y=odeist, times=arange(0.0, 300.0, 1.0), func=odegrad, parms=
    odecons, tol=5e-06, dt=1.0, hbsize=0, nlag=0, ssc=odestsc)
print(ode_eg.data)
Solve Differential equation using Python PyDDE solver