The DEplot[interactive] System Library
This example page describes the default models (available on startup) for the DEplot[interactive] differential system tool.
>



Damped Harmonic Oscillator


Suppose a mass m , which is attached to a spring which is fixed to a wall, is resting on a table. Let be the point where the mass is in rest position. If the mass is moved so that the spring is stretched or compressed, a force is applied to the mass. Hooke's law states that the force exerted by the spring on the mass is where is the spring constant, that is, the force is proportional to the displacement of the mass from the rest position. Since , that is, force equals mass times acceleration, we have the following second order differential equation which governs the motion of the mass.
When the mass slides over the table there will be a frictional force applied to the mass in the opposite direction of motion. Assuming it is proportional to the velocity of the mass, we obtain the following equation
which governs the motion of the mass. Here , is the constant of spring, is the mass. The term models the friction with the table where is the damping coefficient. One can convert this linear second order differential equation into a system of two first order differential equations by letting , that is, v is the velocity, obtaining:
,
>


>


>




Metastasis of Malignant Tumors


This is a twocompartment model for the metastasis of cancer, that is, the spread of cancer cells through the boundary of an organ tissue. Let be the number of arrested cells (sitting on the boundary of the tissue in compartment A) and be the number of cells that have invaded the target tissue (compartment B). Cells pass from compartment A to B with a rate of . In addition, cells die from compartment A with the rate and die from compartment B with the rate . We obtain the following system of two first order differential equations representing these relations.
, .
>


>


>




Brusselator


The brusselator is a mathematical model for chemical oscillation. The dynamics of this system are given by:
, ,
where and , are concentrations of two chemicals, and and are positive constants.
For certain values of a and b, for example, a = 0.5, b = 2, the brusselator exhibits natural oscillation, that is, a stable limit cycle.
>


>


>




KermackMcKendrick Epidemic Model


KermackMcKendrick model is an SIR (Susceptible, Infected, Recovered) model for the number of people infected with a contagious illness in a closed population over time.
Here we are using a simplified version of the model with no recovery, R, element. So we have
, ,
where 0 < S < 1 is the proportion of the population which is susceptible to the illness, 0 < I < 1 is the proportion of infectious individuals in the population, is the mean transmissions rate, and is the mean rate at which infected individuals either recover or die from the illness.
>


>


The most interesting observation to make about this model is that the S does not drop to zero, that is, not everyone in the population gets the illness, even if is increased (individuals are more infectious) and is decreased (individuals stay infectious for longer). In the plot above we see that when I=0, that is, when all infectious individuals have either died or recovered, the value of S=0.2, that is, about 20% of the population is susceptible, i.e., has not caught the illness. In the plot below it is about S=0.02 or 2%. This correctly predicts what has been observed in real life.
>


>




Van der Pol


The van der Pol equation describes the selfsustaining oscillations in which the energy is fed into small oscillations and removed from large oscillations. This equation arises in the study of circuits containing vacuum tubes. The Van der Pol equation is given by:
,
where , a positive constant, describes the current in a triode oscillator. When the equation reduces to the simple harmonic oscillator. One can write this second order differential equation as a system of differential equations:
, .
>


>


>




House warming model


The following system of differential equations describes the relationship between temperatures of rooms in a house with a furnace.
A(t) is the temperature for room A, B(t) is the temperature for room B, F is the heating rate of the furnace located in room B, and Am is the outside temperature.
,
The parameters , and k3 describe the rates of exchange for heat between the rooms A,and B,and outside.
>


>




Richardson's Arms Race Model


Consider two competing nations. Both nations are self defensive and fight back to protect their nation. Both nations maintain army and stock weapons. When one nation expands their army the other nation finds it offensive. Let x(t) and y(t) represent the yearly rate of armament expenditures of the two nations in some standard unit. One can model the Mutual Fear factor for each nation. Assume that extensive armament expenditures create a drag on the nation's economy. And also assume that the each nation's mutual fear's rate is directly proportional to the expenditure of the other nation. In addition, assume that there are constants r and s indicating the grievance of one country toward the other country. Note that negative values for s and r represent feelings of good will.
The Model will be:
where a, b, m, and n are positive constant.
>


>



Return to Index for Example Worksheets
