function s = sdof_stiff2_free_function(t,w,flag,b)
%
s = zeros(2,1);
s(1) = w(2);
s(2) = -b*(w(1)^2)*sign(w(1));    