No BSD License  

Highlights from
Numerical Analysis and Graphic Visualization with MATLAB

f_def(y, t)
% Used in L10_3
function f = f_def(y, t)
a=5; c=20;
f = [y(2),  (-a*abs(y(2))*y(2) - c*y(1))]';

Contact us at files@mathworks.com