function dy = msd(t,y)
global z w

% ddy + 2 * z * w * dy + w * w * y = 0

dy = zeros(2,1);

dy(1) = y(2);
dy(2) = - 2 * z * w * y(2) - w * w * y(1);
