Setting Up the ProblemThe first argument given to DSolve is the differential equation, the second argument is the unknown function, and the last argument identifies the independent variable. Here is the input for solving for a first-order linear ODE using DSolve. The variable sol identifies the solution for use in further work. In[423]:=  |
Out[423]=
|
The output of DSolve is a list of solutions for the differential equation. Here, since the equation is of order 1 and is linear, there is only one solution:  . The solution has an undetermined constant C[1] because no initial condition was specified. The solution can be extracted from the list of solutions using a part specification. This extracts the solution. In[424]:=  |
Out[424]=
|
This form of the solution is useful for finding  itself, but not for finding derivatives of  or the value of  at a point. This shows the value of  given by the solution. In[425]:=  |
Out[425]=
|
The solution does not apply to  or  because the solution is a rule for  only. In[426]:=  |
Out[426]=
|
In[427]:=  |
Out[427]=
|
If the solution will be used in further work, it is best to specify the unknown function using y rather than  . This gives the solution using pure functions of the type Function[x, expr]. Here, the unknown function is specified as  . The solution is a pure function. In[428]:=  |
Out[428]=
|
When the solution is in the form of pure functions, expressions can be found for derivatives of  and for the values of  at specific points. This gives expressions for  ,  , and  . In[429]:=  |
Out[429]=
|
In[430]:=  |
Out[430]=
|
In[431]:=  |
Out[431]=
|
When a problem has multiple solutions, you can pick out individual solutions from the solution list or you can work directly with the list. This solves a nonlinear first-order equation. There are two solutions. In[432]:=  |
Out[432]=
|
The solutions can be extracted using part specifications. In[433]:=  |
Out[433]=
|
In[434]:=  |
Out[434]=
|
This returns a list of both expressions. In[435]:=  |
Out[435]=
|
To solve a system of equations, the first argument to DSolve must be a list of the equations and the second argument must be a list of the unknown functions. Here is an example of a system of first-order linear equations with three unknowns. Because this system is linear, there is only one solution. In[436]:=  |
Out[436]=
|
Each solution to the system is a list of replacement rules for the unknown functions. The expressions for the unknown functions can be extracted as in previous examples. This gives a list of the expressions for the unknown functions. Simplify is used to return the expressions in a compact form. In[437]:=  |
Out[437]=
|
If initial conditions are prescribed for the problem, some or all of the undetermined constants can be eliminated. Here the value of the unknown function and its derivative are both prescribed at the initial point. In[438]:=  |
Out[438]=
|
If only the initial value is specified, then the solution still contains an arbitrary constant. In[439]:=  |
Out[439]=
|
For a partial differential equation, the third argument to DSolve is a list of the independent variables for the equation. This solves a PDE with independent variables x and y. C[1] represents an arbitrary function of  . In[440]:=  |
Out[440]=
|
A differential-algebraic equation is specified in the same way as a system of ordinary differential equations. Here is an example of a DAE with an initial condition. In[441]:=  |
Out[441]=
|
Note that it is not always possible to give the solutions for a problem in explicit form. In this case, the solution is given using an unevaluated Solve object or using InverseFunction. The solution to this equation is not available explicitly. The output represents an implicit solution. In[442]:=  |
Out[442]=
|
The solution can be extracted as usual with a part specification. In[443]:=  |
Out[443]=
|
The solutions to this equation are given as InverseFunction objects. In[444]:=  |
Out[444]=
|
Each solution can be rewritten as an implicit equation by eliminating the InverseFunction object as follows. In[445]:=  |
Out[445]=
|
In[446]:=  |
Out[446]=
|
| 
