# Differential Equations

You can use the

*Mathematica* function

DSolve to find symbolic solutions to ordinary and partial differential equations.

Solving a differential equation consists essentially in finding the form of an unknown function. In

*Mathematica*, unknown functions are represented by expressions like

. The derivatives of such functions are represented by

,

, and so on.

The

*Mathematica* function

DSolve returns as its result a list of rules for functions. There is a question of how these functions are represented. If you ask

DSolve to solve for

, then

DSolve will indeed return a rule for

. In some cases, this rule may be all you need. But this rule, on its own, does not give values for

or even

. In many cases, therefore, it is better to ask

DSolve to solve not for

, but instead for

itself. In this case, what

DSolve will return is a rule which gives

as a pure function, in the sense discussed in

"Pure Functions".

If you ask

DSolve to solve for

, it will give a rule specifically for

.

Out[1]= | |

The rule applies only to

itself, and not, for example, to objects like

or

.

Out[2]= | |

If you ask

DSolve to solve for

, it gives a rule for the object

on its own as a pure function.

Out[3]= | |

Now the rule applies to all occurrences of

.

Out[4]= | |

Substituting the solution into the original equation yields

True.

Out[5]= | |

DSolve[eqn,y[x],x] | solve a differential equation for |

DSolve[eqn,y,x] | solve a differential equation for the function y |

Getting solutions to differential equations in different forms.

In standard mathematical notation, one typically represents solutions to differential equations by explicitly introducing "dummy variables" to represent the arguments of the functions that appear. If all you need is a symbolic form for the solution, then introducing such dummy variables may be convenient. However, if you actually intend to use the solution in a variety of other computations, then you will usually find it better to get the solution in pure-function form, without dummy variables. Notice that this form, while easy to represent in

*Mathematica*, has no direct analog in standard mathematical notation.

DSolve[{eqn_{1},eqn_{2},...},{y_{1},y_{2},...},x] |

| solve a list of differential equations |

Solving simultaneous differential equations.

This solves two simultaneous differential equations.

Out[6]= | |

*Mathematica* returns two distinct solutions for

in this case.

Out[7]= | |

You can add constraints and boundary conditions for differential equations by explicitly giving additional equations such as

.

This asks for a solution which satisfies the condition

.

Out[8]= | |

If you ask

*Mathematica* to solve a set of differential equations and you do not give any constraints or boundary conditions, then

*Mathematica* will try to find a

*general solution* to your equations. This general solution will involve various undetermined constants. One new constant is introduced for each order of derivative in each equation you give.

The default is that these constants are named

C[n], where the index

n starts at

for each invocation of

DSolve. You can override this choice, by explicitly giving a setting for the option

GeneratedParameters. Any function you give is applied to each successive index value

n to get the constants to use for each invocation of

DSolve.

The general solution to this fourth-order equation involves four undetermined constants.

Out[9]= | |

Each independent initial or boundary condition you give reduces the number of undetermined constants by one.

Out[10]= | |

You should realize that finding exact formulas for the solutions to differential equations is a difficult matter. In fact, there are only fairly few kinds of equations for which such formulas can be found, at least in terms of standard mathematical functions.

The most widely investigated differential equations are linear ones, in which the functions you are solving for, as well as their derivatives, appear only linearly.

This is a homogeneous first-order linear differential equation, and its solution is quite simple.

Out[11]= | |

Making the equation inhomogeneous leads to a significantly more complicated solution.

Out[12]= | |

If you have only a single linear differential equation, and it involves only a first derivative of the function you are solving for, then it turns out that the solution can always be found just by doing integrals.

But as soon as you have more than one differential equation, or more than a first-order derivative, this is no longer true. However, some simple second-order linear differential equations can nevertheless be solved using various special functions from

"Special Functions". Indeed, historically many of these special functions were first introduced specifically in order to represent the solutions to such equations.

This is Airy's equation, which is solved in terms of Airy functions.

Out[13]= | |

This equation comes out in terms of Bessel functions.

Out[14]= | |

This requires Mathieu functions.

Out[15]= | |

And this Legendre functions.

Out[16]= | |

Occasionally second-order linear equations can be solved using only elementary functions.

Out[17]= | |

Beyond second order, the kinds of functions needed to solve even fairly simple linear differential equations become extremely complicated. At third order, the generalized Meijer G function

MeijerG can sometimes be used, but at fourth order and beyond absolutely no standard mathematical functions are typically adequate, except in very special cases.

Here is a third-order linear differential equation which can be solved in terms of generalized hypergeometric functions.

Out[18]= | |

This requires more general Meijer G functions.

Out[19]= | |

For nonlinear differential equations, only rather special cases can usually ever be solved in terms of standard mathematical functions. Nevertheless,

DSolve includes fairly general procedures which allow it to handle almost all nonlinear differential equations whose solutions are found in standard reference books.

First-order nonlinear differential equations in which

does not appear on its own are fairly easy to solve.

Out[20]= | |

This Riccati equation already gives a significantly more complicated solution.

Out[21]= | |

This Bernoulli equation, however, has a fairly simple solution.

Out[22]= | |

An

order Bernoulli equation typically has

distinct solutions.

Out[23]= | |

This Abel equation can be solved, but only implicitly.

Out[24]= | |

In practical applications, it is quite often convenient to set up differential equations that involve piecewise functions. You can use

DSolve to find symbolic solutions to such equations.

This equation involves a piecewise forcing function.

Out[25]= | |

Here the solution is explicitly broken into three cases.

Out[26]= | |

Beyond ordinary differential equations, one can consider

*differential-algebraic equations* that involve a mixture of differential and algebraic equations.

This solves a differential-algebraic equation.

Out[27]= | |

DSolve[eqn,y[x_{1},x_{2},...],{x_{1},x_{2},...}] |

| solve a partial differential equation for |

DSolve[eqn,y,{x_{1},x_{2},...}] | solve a partial differential equation for the function y |

Solving partial differential equations.

DSolve is set up to handle not only

*ordinary differential equations* in which just a single independent variable appears, but also

*partial differential equations* in which two or more independent variables appear.

This finds the general solution to a simple partial differential equation with two independent variables.

Out[28]= | |

Here is the result represented as a pure function.

Out[29]= | |

The basic mathematics of partial differential equations is considerably more complicated than that of ordinary differential equations. One feature is that whereas the general solution to an ordinary differential equation involves only arbitrary

*constants*, the general solution to a partial differential equation, if it can be found at all, must involve arbitrary

*functions*. Indeed, with

independent variables, arbitrary functions of

arguments appear.

DSolve by default names these functions

C[n].

Here is a simple PDE involving three independent variables.

Out[30]= | |

The solution involves an arbitrary function of two variables.

Out[31]= | |

Here is the one-dimensional wave equation.

Out[32]= | |

The solution to this second-order equation involves two arbitrary functions.

Out[33]= | |

For an ordinary differential equation, it is guaranteed that a general solution must exist, with the property that adding initial or boundary conditions simply corresponds to forcing specific choices for arbitrary constants in the solution. But for partial differential equations this is no longer true. Indeed, it is only for linear partial differential and a few other special types that such general solutions exist.

Other partial differential equations can be solved only when specific initial or boundary values are given, and in the vast majority of cases no solutions can be found as exact formulas in terms of standard mathematical functions.

Since

and its derivatives appear only linearly here, a general solution exists.

Out[34]= | |

This weakly nonlinear PDE turns out to have a general solution.

Out[35]= | |

Here is a nonlinear PDE which has no general solution.

Out[36]= | |