function dy = ode3(x, y, mu)

dy = mu*(sin(x) - y);

end