function s = pendulum_free_ODE(t,w,flag,b)
%
s = zeros(2,1);
s(1) = w(2);
s(2) = -b*sin(w(1));    