integrate the solutions of a differential operator or equation
list containing two names
homogeneous linear differential equation
The input L is a differential operator. This procedure computes an operator M of minimal order such that any solution of L has an antiderivative which is a solution of M.
If the order of L equals the order of M then the output is a list [M, r] such that r(f) is an antiderivative of f and also a solution of M for every solution f of L. If the order of L is not equal to M then only M is given in the output. In this case M equals L⁢Dt where Dt is the derivation.
The argument domain describes the differential algebra. If this argument is the list Dt,t, then the differential operators are notated with the symbols Dt and t. They are viewed as elements of the differential algebra C⁡t Dt where C is the field of constants.
If the argument domain is omitted then the differential specified by the environment variable _Envdiffopdomain is used. If this environment variable is not set then the argument domain may not be omitted.
Instead of a differential operator, the input can also be a linear homogeneous differential equation having rational function coefficients. In this case the second argument must be the dependent variable.
This function is part of the DEtools package, and so it can be used in the form integrate_sols(..) only after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[integrate_sols](..).
L ≔ Dz+1z
The result was only the operator M. For an example where both M and r are obtained, consider:
L ≔ Dz2−3⁢Dz2⁢z⁢z−1−z−2⁢z2−z+1z⁢z−12
Regarding the meaning of the second element in the output of integrate_sols, consider the following second order ODE and its solution.
ode ≔ 27⁢x2+4⁢ⅆⅆx⁢ⅆⅆx⁢y⁡x+27⁢x⁢ⅆⅆx⁢y⁡x−3⁢y⁡x
sol ≔ dsolve⁡ode,y⁡x
The ODE satisfied by the integral of sol is given by the first element of the output of integrate_sols while the formula for computing its solution (as a function of sol) is given by the second element.
int_s ≔ integrate_sols⁡ode,y⁡x:
ode2 ≔ int_s1
integral_of_sol ≔ int_s2
So the solution to ode2 is obtained by substituting sol into integral_of_sol.
sol2 ≔ y⁡x=simplify⁡eval⁡integral_of_sol,sol
That sol2 is the integral of sol can also be verified as follows: Differentiate sol2 and verify that it is equal to sol.
Abramov, S.A., and van Hoeij, M. "A method for the Integration of Solutions of Ore Equations." ISSAC '97 Proceedings, pp. 172-175. 1997.
van der Put, M., and Singer, M. F. Galois Theory of Linear Differential Equations, Vol. 328. Springer, 2003. An electronic version of this book is available at http://www4.ncsu.edu/~singer/ms_papers.html.
Download Help Document