% This ODE is a goofy nonlinear oscillator.
%     d/dt (x,v) = (v, x - x^3 - mu*v).
%
function dxdt = goofy(state,t,mu)

x = state(1);
v = state(2);

dxdt = [ v, x - x^3 - mu*v ];