grid on function dydt=BonusPointQ2Part2(K,theta) dydt=zeros(4,1); w1=0.26 w2=0.24 w3=0.27 w4=0.25 dydt(1)=w1 + (K./4).* (sin(theta(1)-theta(1)) + sin(theta(2)-theta(1)) + sin(theta(3)-theta(1)) + sin(theta(4)-theta(1))); dydt(...