function xdot = f(t,x)
  xdot = [x(2);
          (1/5)*(sin(t)-4*x(1)-7*x(2))];
endfunction