function dx = my_system( t, x )

  global a b f m
  dx = zeros(2,1);
  dx(1) = x(1)*(1-x(1)-a*x(2));
  dx(2) = x(2)*(1-x(2)- (b+m*sin(f*t)) *x(1));