Jim Belk University of Glasgow

Homework 10

Due Date: Saturday, December 3

Instructions: Feel free to work together with other students in the class, though you must turn in your own copy of the solutions, and you must acknowledge anyone that you worked with. You can turn in your homework assignment by e-mailing me your solutions.

  1. In the Lotka-Volterra predator-prey model, the population of prey tends to decrease when there are more predators, while the population of predators tends to increase when there are more prey. However, another common situation in ecology is that of competition between two species, where a greater population of either species has an adverse effect on the population of the other.
    1. Write a simple system of differential equations for the populations \(x\) and \(y\) of two species in competition. You should assume that either population would experience exponential growth in the absence of the other species, but that interactions between the two species are harmful to both. Your model should involve four constants \(A\), \(B\), \(C\), and \(D\).
    2. Find the nonzero equilibrium point \((x^*,y^*)\) for your model in terms of \(A\), \(B\), \(C\), and \(D\).
    3. Compute the eigenvalues of the Jacobian for the equilibrium point you found in part (b) in terms of \(A\), \(B\), \(C\), and \(D\). What type of equilibrium is this?
    4. Use VectorPlot to plot the vector field for your system in part (a) in the case where \(A=B=C=D=1\). Your plot should include the range \(0\leq x \leq 2\) and \(0 \leq y \leq 2\), and should use the option VectorScale -> {Automatic, Automatic, None} to make the arrows visible.
    5. Based on your plot from part (d), what tends to happen to the population of each species over time in the case where \(A=B=C=D=1\)? Does it depend on the initial populations of the species? Explain.
  2. In a diatomic molecule, the two atoms are bound together by a long-range attractive force offset by a short-range repulsive force. One simple model for such a bond (known as the Morse potential model) is given by \[ \frac{d^2 r}{dt^2} \,=\, ae^{-2kr} - be^{-kr}, \] where \(r\) is the distance between the two nuclei, \(k\) is a constant with units of inverse length, and \(a\) and \(b\) are constants satisfying \(a>b\).
    1. Which term in the above equation represents the repulsive force? Which term represents the attractive force? In what sense does the attractive force have a longer range than the repulsive force?
    2. Find the equilibrium value \(r^*\) for the distance between the nuclei.
    3. Let \(u=r -r^*\). Show that \(u\) satisfies the differential equation \[ \frac{d^2 u}{dt^2} \,=\, \frac{b^2}{a}\bigl(e^{-2ku}-e^{-ku}\bigr). \]
    4. Find the linearization of the equation in part (c) in the case where \(u\approx 0\), and use your answer to determine the period of small oscillations of \(r\) around the equilibrium.
    5. In the case where \(a=2\) and \(b=k=1\), use NDSolve and Plot to plot the solution with initial conditions \(r(0) = 0.7\) and \(r'(0)=0\) for \(0\leq t \leq 30\).
    6. Does the period of the solution in part (e) agree with your answer to part (d)? What is the amplitude of the oscillations? Do the oscillations appear to be sinusoidal?
    7. Repeat parts (e) and (f) for the solution with initial conditions \(r(0) = 2\) and \(r'(0)=0\).