dsolve/inttrans  integral transform methods for solutions of ordinary differential equations

Calling Sequence


dsolve(deqns, vars, method=transform)
dsolve(deqns, vars, method=transform, options)


Parameters


deqns



ordinary differential equation in vars, or set of equations and/or initial conditions

vars



variable or set of variables to be solved for

method=transform



literal equation, one of laplace, fourier, fouriercos, fouriersin

options



optional equations





Description


•

Invoking the dsolve command with the options method=transform causes a solution to be found by way of one of the integral transform methods.

•

When adding the optional equation pointwise=true to dsolve, the desired solution should satisfy the initial conditions in a pointwise manner. This guarantees that the solution evaluates to the specified initial conditions at the initial point. By default this is false, which means that the input conditions are satisfied in a distributional manner (that is, applied as specified by the transform itself). This is most relevant when the input differential equations have distributional functions present (Dirac and Heaviside functions). For an example that illustrates the difference between these approaches see the Examples section below.

•

Caution must be used, as some of the transforms (most notably the fouriertype transforms) are only well defined for solutions with certain asymptotic properties.


For example, when solving a second order constant coefficient ode using the fouriersin transform, the solution y(t) must vanish as t > infinity. This restriction is in fact one of the boundary conditions of the problem, so only one other boundary condition can be applied.

•

Some guidelines for using transforms for the solution of secondorder constant coefficient ODE are as follows.


fourier transform: the solution y(t) of the ODE (defined for all t) must be bounded as t approaches +/ infinity. No other boundary conditions can be applied.


fouriersin transform: the solution y(t) of the ODE (defined for t>=0) must be bounded as t approaches +infinity, and y(0) can be specified.


fouriercos transform: the solution y(t) of the ODE (defined for t>=0) must be bounded as t approaches +infinity, and D(y)(0) can be specified.


laplace transform: the solution y(t) of the ODE (defined for t>=0) must be subexponential, and both y(0) and D(y)(0) can be specified.

•

Additionally, care must be taken when choosing integral transform methods to obtain a solution to an ODE or ODE system. For example, use of fouriertype transforms for problems with polynomial coefficients may require the solution of another differential equation of higher differential order (see example below).



Examples


>


 (1) 
>


 (2) 
>


 (3) 
>


 (4) 
>


 (5) 
>


>


 (6) 
>


 (7) 
>


 (8) 
>


 (9) 
However this gives no solution (overspecified):
>


Now a counterproductive transform:
>


 (10) 
>


 (11) 
You would now have an 8th order ODE to solve.
Attempting to use the fourier transform with initial conditions yields no solution:
>


 (12) 
>


The following ODE has different solutions depending on whether 'pointwise'=true has been specified.
>


 (13) 
>


 (14) 
So the solution as provided does not satisfy the initial conditions in a pointwise fashion, as can be seen from
>


 (15) 
To instead enforce the conditions to be pointwise satisfied, the following can be used:
>


 (16) 
>


 (17) 


See Also


dsolve[ck45], dsolve[classical], dsolve[dverk78], dsolve[gear], dsolve[lsode], dsolve[numeric], dsolve[piecewise], dsolve[rkf45], dsolve[taylorseries], inttrans, plots[odeplot]

