Why Differential Equations

INTRODUCTION

MODELING

SOLUTIONS

MODELING EXAMPLE

TORRICELLI EXPERIMENT

SEPARATION OF VARIABLES

Given that [Graphics:1txgr114.gif] and that [Graphics:1txgr115.gif] can be factored into a product of two functions such as [Graphics:1txgr116.gif], we can always separate the variables and attempt the integration. Integrate both sides with respect to [Graphics:1txgr117.gif] gives: [Graphics:1txgr118.gif]. On the left side of this integral equation, we use the chain rule for a single variable along with the change of variable as follows: We let [Graphics:1txgr119.gif]. Then recall that [Graphics:1txgr120.gif], which is the derivative of the inner function,[Graphics:1txgr121.gif] times the derivative of what is inside, [Graphics:1txgr122.gif]. Making these substitutions in the above and dividing both sides by [Graphics:1txgr123.gif] gives us the integral equation written below. This equation is ready for integration. Of course now we must consider if the integration can be done analytically? Often it can't, and we must resort to an numerical integration.

[Graphics:1txgr124.gif]

If the original differential equation is autonomous, that is [Graphics:1txgr125.gif] or if [Graphics:1txgr126.gif] then it is always possible to separate the variables and attempt the integration. We will consider some examples of these types.

Example 1: [Graphics:1txgr127.gif].

This is obviously a product of the dependent and independent variables. So we can separate the variable by inspection and write

[Graphics:1txgr128.gif].
Both sides of this equation can now be integrated. And we have for a solution

[Graphics:1txgr129.gif] or

[Graphics:1txgr130.gif]

We see from these results that there are an infinite number of solutions for this equation depending on the value of "C." If [Graphics:1txgr131.gif], we find that [Graphics:1txgr132.gif]. Then our solution for this specific initial value problem is

[Graphics:1txgr133.gif]

In this particular example, we must keep in mind that initial conditions such as [Graphics:1txgr134.gif] will not work because with the inverse the sin function only takes on values in the interval [-1,1]. This is also apparent from the [Graphics:1txgr135.gif]. It is possible to solve this problem using the Mathematica built in differential equation solver "DSolve." Note carefully the syntax. We have a list of two items enclosed in {}, the differential equation and the initial condition, both with a double equal. Next comes the dependent variable and finally, the independent variable. We see that we got the same results as found above. The solution is the same as found above. Also, note that the plot oscillates between plus and minus one.

	[Graphics:1txgr136.gif]
	[Graphics:1txgr137.gif]
	[Graphics:1txgr138.gif]
	[Graphics:1txgr139.gif]
	[Graphics:1txgr140.gif]

               [Graphics:1txgr141.gif]

Example 2: [Graphics:1txgr142.gif]

This differential equation is autonomous and therefore separable. The left side may be integrated by using the method of partial fractions. We can either do this by hand or we may use a Mathematica function.

[Graphics:1txgr143.gif]

	[Graphics:1txgr144.gif]
	[Graphics:1txgr145.gif]
	[Graphics:1txgr146.gif]
	[Graphics:1txgr147.gif]

So we see with a little help from Mathematica that the solutions are:

     [Graphics:1txgr148.gif].

Now if we convert this logrithmetic function to the exponential form, we arrive at a solution for y.

     [Graphics:1txgr149.gif]

Now "[Graphics:1txgr150.gif]" is the effective integration constant. If we rearrange and solve for y, we have

     [Graphics:1txgr151.gif].

When [Graphics:1txgr152.gif], we have

     [Graphics:1txgr153.gif].

Solving gives [Graphics:1txgr154.gif]. The particular solution for this specific initial value problem is as follows:

     [Graphics:1txgr155.gif].

Let us again use the Mathematica built in function to solve this initial value problem. In this example we will set up the equation along with the initial value in a separate instruction. We arrive at the same answer.

	[Graphics:1txgr156.gif]
	[Graphics:1txgr157.gif]

Next, let us plot this solution.

	[Graphics:1txgr158.gif]
	[Graphics:1txgr159.gif]

               [Graphics:1txgr160.gif] In the plot we see that there is a verticle asymtote. We can set the denominator equal to zero and solve for x to find the asymtote.

	[Graphics:1txgr161.gif]
	[Graphics:1txgr162.gif]
	[Graphics:1txgr163.gif]

Also, looking at the original problem, and the solution [Graphics:1txgr164.gif] we must be aware that the values of y = 0 or y = -4 are not allowed for the initial conditions, because y = 0 results in an undefined log. With [Graphics:1txgr165.gif], we get zero in the denominator. Differential equations often have these sublities and must be scrutinized before accepting any solution. However, notice that these are both satisfactory solutions to the original differential equation. These are equilibrium solutions. Thus the general solution consists of functions of the form [Graphics:1txgr166.gif] along with the equilibrium solutions [Graphics:1txgr167.gif]. We did not get the two equilibrium solutions from the results obtained just from the separation of variables and integration.

EXERSISES

SOLUTIONS TO THE EXERSISES

MAIN INDEX


PEVIOUS -- INDEX -- NEXT