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 .
In[1]:=
Click for copyable input
Out[1]=
The rule applies only to itself, and not, for example, to objects like or .
In[2]:=
Click for copyable input
Out[2]=
If you ask DSolve to solve for , it gives a rule for the object on its own as a pure function.
In[3]:=
Click for copyable input
Out[3]=
Now the rule applies to all occurrences of .
In[4]:=
Click for copyable input
Out[4]=
Substituting the solution into the original equation yields True.
In[5]:=
Click for copyable input
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[{eqn1,eqn2,...},{y1,y2,...},x]
solve a list of differential equations

Solving simultaneous differential equations.

This solves two simultaneous differential equations.
In[6]:=
Click for copyable input
Out[6]=
Mathematica returns two distinct solutions for in this case.
In[7]:=
Click for copyable input
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 .
In[8]:=
Click for copyable input
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.
In[9]:=
Click for copyable input
Out[9]=
Each independent initial or boundary condition you give reduces the number of undetermined constants by one.
In[10]:=
Click for copyable input
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.
In[11]:=
Click for copyable input
Out[11]=
Making the equation inhomogeneous leads to a significantly more complicated solution.
In[12]:=
Click for copyable input
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.
In[13]:=
Click for copyable input
Out[13]=
This equation comes out in terms of Bessel functions.
In[14]:=
Click for copyable input
Out[14]=
This requires Mathieu functions.
In[15]:=
Click for copyable input
Out[15]=
And this Legendre functions.
In[16]:=
Click for copyable input
Out[16]=
Occasionally second-order linear equations can be solved using only elementary functions.
In[17]:=
Click for copyable input
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.
In[18]:=
Click for copyable input
Out[18]=
This requires more general Meijer G functions.
In[19]:=
Click for copyable input
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.
In[20]:=
Click for copyable input
Out[20]=
This Riccati equation already gives a significantly more complicated solution.
In[21]:=
Click for copyable input
Out[21]=
This Bernoulli equation, however, has a fairly simple solution.
In[22]:=
Click for copyable input
Out[22]=
An ^(th) order Bernoulli equation typically has distinct solutions.
In[23]:=
Click for copyable input
Out[23]=
This Abel equation can be solved, but only implicitly.
In[24]:=
Click for copyable input
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.
In[25]:=
Click for copyable input
Out[25]=
Here the solution is explicitly broken into three cases.
In[26]:=
Click for copyable input
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.
In[27]:=
Click for copyable input
Out[27]=
DSolve[eqn,y[x1,x2,...],{x1,x2,...}]
solve a partial differential equation for
DSolve[eqn,y,{x1,x2,...}]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.
In[28]:=
Click for copyable input
Out[28]=
Here is the result represented as a pure function.
In[29]:=
Click for copyable input
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.
In[30]:=
Click for copyable input
Out[30]=
The solution involves an arbitrary function of two variables.
In[31]:=
Click for copyable input
Out[31]=
Here is the one-dimensional wave equation.
In[32]:=
Click for copyable input
Out[32]=
The solution to this second-order equation involves two arbitrary functions.
In[33]:=
Click for copyable input
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.
In[34]:=
Click for copyable input
Out[34]=
This weakly nonlinear PDE turns out to have a general solution.
In[35]:=
Click for copyable input
Out[35]=
Here is a nonlinear PDE which has no general solution.
In[36]:=
Click for copyable input
Out[36]=
New to Mathematica? Find your learning path »
Have a question? Ask support »