function dvdt = flinsys(t,v,p)
   x = v(1);
   y = v(2);

   omega = p(1);
   a     = p(2);
   e     = p(3);

   dvdt = [ y; -e*y - x - e*x^3 + e*a*cos(omega*t)];