% Volterra vector field (System 1.9, Kuznetsov, 2nd ed.)

  function v = voltfun(x,gamma) ;
  xdot = x(1) - x(1)*x(2) ;
  ydot = -gamma*x(2) + x(1)*x(2) ;
  v = [xdot ; ydot] ; 
