problem has the form
      <pre>
        du/dt =         del u + ( 1 - u ) * u        - v * h(u/alpha)

        dv/dt = delta * del v     - gamma * v + beta * v * h(u/alpha)
      </pre>
      with initial conditions:
      <pre>
        u(x,0) = u0(x)
        v(x,0) = v0(x)
      </pre>
      and boundary conditions at the left and right endpoints [A,B]:
      <pre>
        du/dx = 0
        dv/dx = 0
      </pre>
      The Type II functional response employed here is
      <pre>
        h(eta) = eta / ( 1