## Mathematics & Science
Learning Center |
|||

## Applications of Systems of Differential Equations## Combat Models: Who Will Win The Battle?## Guerilla Combat Model
In this situation we suppose that both of the armies involved are waging guerilla warfare on the other. A possible model for the scenario may be the following system of differential equations: x′(t) = -0.1 x(t) y(t)
y′(t) = -0.2 x(t) y(t)
x(0) = x_{o}, y(0)
=y_{o}
where
In a combat model,
Now for our analysis of the model using Switch to Mathematica by clicking on the button at left. This will
open up a fresh notebook for you. Then return here immediately for further instructions.-
Let's begin by considering a situation in which both forces have an equal number of troops, say *x*(0) = 25,*y*(0) = 25, where both are measured in thousands. Use the NDSolve command, just as you did in the previous two labs, to solve the above system of differential equations, along with the initial conditions just given. Use the interval 0 ≤ t ≤ 30, and read the solution into the variable guerilla1. -
That shouldn't have taken too long. Now make a phase-plane diagram by using ParametricPlot, just like you did when analysing predator-prey problems in the previous lab. Remember, you'll be plotting {x[t],y[t]}/.guerilla1[[1]] this time, and your *t*interval is from 0 to 30. Include a PlotRange->All option in your command, and also an AxesOrigin->{0,0} option to force the axes to cross in the location you are used to, and assign the graph to the variable guerillaplot1. Hmm, just about the most boring graph that we've seen in this course. Anyway, now for the important question: Who won? Remember that your initial condition was (25, 25). This means that you start your solution (with respect to time) in the top right corner of the graph, and as time passes you move down and to the left along the solution line. You can determine who won by observing where this solution "curve" ended up. It looks like it stops at about (12.5, 0) on my plot. That would be 12,500 Xanudians, and zero Yarovians. I'd say that the Xanudians won! This is in agreement with our initial observation that the Xanudians are better fighters. By the way, a couple of remarks: The model ceases to make sense if your solution curve ventures out of the first quadrant. (Why?) And, in general, what does the intersection of the solution curve with an axis (in standard position) mean in the "real world" for this model? Well, maybe we should let the Yarovians have more troops to start with. Redo steps 1 and 2 above, but this time with initial conditions *x*(0) = 25,*y*(0) = 50. Of course, name your results guerilla2, and guerillaplot2 this time.-
Who won the battle now? (You should have found that the two forces completely wiped each other out, kind of along the lines of USA/USSR nuclear arms policy during the cold war.) This result probably didn't surprise you. Let's redo steps 1 and 2 again, this time with initial conditions of *x*(0) = 25,*y*(0) = 75. In other words we're giving the Yarovian force an even larger army than in the last scenario. Rename your variables guerilla3, and guerillaplot3 this time. -
At last, the Yarovians score a victory (about 25,000 of their troops survived)! What seems to be the basic requirement for a Yarovian victory in this model? It looks like they have to have more than twice the number of troops as the Xanudians, probably due to the fact that the Xanudians fight twice as effectively. -
Finally, use the Show command to display all three graphs on the same coordinate system. (Is there any wonder that guerilla vs. guerilla warfare is said to obey the **linear combat law**?)
OK, that takes care of the guerilla vs. guerilla scenario. Let's move on to a slightly more interesting situation... |
|||

ODE Laboratories: A Sabbatical Project by Christopher A. Barker©2017 San Joaquin Delta College, 5151 Pacific Ave., Stockton, CA 95207, USA e-mail: cbarker@deltacollege.edu |