This is documentation for Mathematica 6, which was
based on an earlier version of the Wolfram Language.
View current documentation (Version 11.2)

Second-Order PDEs

The general form of a linear second-order PDE is
Here u=u (x, y), and a, b, c, d, e, f, and g are functions of x and y only—they do not depend on u. If g=0, the equation is said to be homogeneous.
The first three terms containing the second derivatives are called the principal part of the PDE. They determine the nature of the general solution to the equation. In fact, the coefficients of the principal part can be used to classify the PDE as follows.
The PDE is said to be elliptic if b2-4ac<0. The Laplace equation has a=1, b=0, and c=1 and is therefore an elliptic PDE.
The PDE is said to be hyperbolic if b2-4ac>0. The wave equation has a=1, b=0, and c=-1 and is therefore a hyperbolic PDE.
The PDE is said to be parabolic if b2-4ac=0. The heat equation has a=1, b=0, and c=0 and is therefore a parabolic PDE.
DSolve can find the general solution for a restricted type of homogeneous linear second-order PDEs; namely, equations of the following form
Here a, b, and c are constants. Thus, DSolve assumes that the equation has constant coefficients and a vanishing non-principal part.
Following are some examples of the three basic types (elliptic, hyperbolic, and parabolic) and an explanation of their significance.
Here is the general solution for Laplace's equation, an elliptic PDE.
In[1]:=
Click for copyable input
In[2]:=
Click for copyable input
Out[2]=
This general solution contains two arbitrary functions, C[1] and C[2]. The arguments of these functions, y+x and y-x, indicate that the solution is constant along the imaginary straight line y=-x+ when C[2]=0 and along y=x+ when C[1]=0 . These straight lines are called characteristic curves of the PDE. In general, elliptic PDEs have imaginary characteristic curves.
Here is another elliptic PDE.
In[3]:=
Click for copyable input
Out[3]=
In[4]:=
Click for copyable input
Note the imaginary characteristic curves for the equation.
In[5]:=
Click for copyable input
Out[5]=
The solution is verified as follows.
In[6]:=
Click for copyable input
Out[6]=
This finds the general solution of the wave equation, a hyperbolic PDE. The constant c in the wave equation represents the speed of light and is set to 1 here for convenience.
In[7]:=
Click for copyable input
In[8]:=
Click for copyable input
Out[8]=
The characteristic lines for the wave equation are x=k+t and x=k-t where k is an arbitrary constant. Hence the wave equation (or any hyperbolic PDE) has two families of real characteristic curves. If initial conditions are specified for the wave equation, the solution propagates along the characteristic lines. Also, any fixed pair of characteristic lines determine the null cone of an observer sitting at their intersection.
Here is another example of a hyperbolic PDE.
In[9]:=
Click for copyable input
Out[9]=
In[10]:=
Click for copyable input
Notice that the equation has two families of real characteristics.
In[11]:=
Click for copyable input
Out[11]=
The solution can be verified as follows.
In[12]:=
Click for copyable input
Out[12]=
Finally, here is an example of a parabolic PDE.
In[13]:=
Click for copyable input
Out[13]=
In[14]:=
Click for copyable input
In[15]:=
Click for copyable input
Out[15]=
The equation has only one family of real characteristics, the lines y=5x+. In fact, any parabolic PDE has only a single family of real characteristics.
The solution can be verified as follows.
In[16]:=
Click for copyable input
Out[16]=
The heat equation is parabolic, but it is not considered here because it has a nonvanishing non-principal part, and the algorithm used by DSolve is not applicable in this case.
This concludes the discussion of partial differential equations.