next up previous
Next: Notation Up: Numerical O.D.E.'s Previous: Numerical O.D.E.'s

Background

Given an ordinary differential equation of the form (5.1), it is natural (or even essential) that the question of whether such a solution, y(t), even exists. In other words, there are two directions of pursuit given an equation of the form (5.1). There are first the questions of when an analytic solution exists and when the solution is unique. An analytic solution refers to a closed-form formula which can be used for computing the value of the solution function y at any point t. In contrast, there are subsequent issues associated with finding a numerical solution; the generation of a table of approximate values of the solution function at a finite set of points.

    To address the issues of existence and uniqueness, consider the simplest form of the first-order ordinary differential equation, namely

  equation839

with f(t) a given function. A natural and intuitive approach to solving this problem comes from calculus. Referred to as the ``Second Fundamental Theorem of Calculus'' in some calculus texts, equation (5.2) is solved completely from the analytical vantage point through integration:

displaymath3468

Using the initial condition tex2html_wrap_inline3490 , the constant of integration is uniquely determined and we obtain the closed-form formula

  equation844

That is, analytically our quest at finding the solution is complete, with the sole caveat that the function f(t) is integrable in some neighborhood of tex2html_wrap_inline3494 .

example848

    In looking to the solution of the more general format of the first-order ordinary differential equation, (5.1), the approach taken is very similar in spirit and the analytical solution of the differential equation becomes a question of judiciously reformulating the problem into an equivalent integration problem. For example, the technique of integrating factors can be used to solve a slightly more general right-hand side of the differential equation than given in (5.2). Consider the ordinary differential equation

  eqnarray854

Moving terms around, the equivalent problem becomes

  equation857

The technique of integrating factors is motivated by the question of whether the above equation could be multiplied through by some function - appropriately referred to as the integrating factor - which would cast the left-hand side of the above expression into the derivative of a common expression. Most commonly, this expression takes the form of a derivative of a product of two functions. If such an integrating factor can be found, then solving the expression is as simple as integrating both sides.

    In detail, consider the task of finding an integrating factor, v(t), which when multiplies (5.5), transforms the left-hand side of the equation to the derivative of the product v(t)y(t). Thus, such an integrating factor would permit

  eqnarray861

Thus, integration of equation (5.6) yields

eqnarray865

If division by v(t) is possible, i.e. tex2html_wrap_inline3508 , then the analytic derivation of solution y(t) is complete:

  equation870

Looking back at the motivation of what we were looking for in an integrating factor and the implication it had on equation (5.6), it is clear that in order for these results to hold together, such an integrating factor needs to satisfy the property

displaymath3472

This can be seen by observing that equality must hold between the two left-hand sides of equation (5.6) above. This is a first-order ordinary differential equation as well, one that is separable. Again, the condition tex2html_wrap_inline3508 that was necessary in obtaining the formula (5.9) appears here in the separable equation technique for solving v'(t) =-a(t)v(t). Recall the separation of variables technique from calculus for solving such an equation, which involves placing all terms involving v on one side of the equation and all non-v terms on the other. In the problem at hand, we assume that tex2html_wrap_inline3508 so that we can divide through by v(t) and obtain the equation

displaymath3473

Using the chain rule, integration of both sides of this expression with respect to t yields

eqnarray882

Hence

displaymath3474

    Now, it's time to wrap up some loose ends in our above derivation. Our assumption that tex2html_wrap_inline3508 here appears to be well-founded inasmuch as the exponential function on the right-hand side of the above calculation is never zero provided a(t) is finite. Moreover, the choice of sign ( tex2html_wrap_inline3530 ) is ultimately irrelevant due to the nature of the problem, (see equation 5.9). Consequently, the choice of the integrating factor v(t) is complete:

  equation891

example895

    At the risk of giving the incorrect impression that all differential equations are uniquely solvable, the problem of when unique solutions to equation (5.1) exist should be made precise. The result from the area of differential equations which addresses exactly this point is embodied below.

  theorem920

    In Theorem 5.1, `` tex2html_wrap_inline3566 '' refers to the partial derivative of the function f(t,z) with respect to the variable z. If partial derivatives are new to you, a partial derivative of a function with respect to the variable z is loosely defined as the regular derivative of the function with respect to the variable z when you consider all other variables (such as t) to be constants.

    To illustrate the ideas presented in Theorem 5.1, we conclude this section with the following example.

example929


next up previous
Next: Notation Up: Numerical O.D.E.'s Previous: Numerical O.D.E.'s
Paul Gray

Wed Oct 28 11:42:13 EST 1998