How do I simulate the N Oscillator Kuramoto Model with no loop

How do I simulate this formula with no loop?

This is what I am doing.

N=10;

theta=0:10;

W=0:10;

dydt=omega+(K/N).*(sin(theta(1)-theta(1))+sin(theta(2)-theta(1))+sin(theta(3)-theta(2))+(sin(theta(4)-theta(1))));