Terminology and Properties of Linear Differential Equations
The highest derivative in a linear differential equation determines its order. The term b(x) is known as the nonhomogeneous term and can be a function of x. When b(x) equals zero, the equation is homogeneous. The corresponding homogeneous equation is derived by setting b(x) to zero. Solutions to a linear differential equation are functions that satisfy the equation, and the solution set of a homogeneous equation forms a vector space. This vector space has a finite dimension equal to the order of the equation. The complete solution to a linear differential equation is constructed by adding a particular solution to the general solution of the associated homogeneous equation.Linear Differential Operators
Linear differential operators are instrumental in analyzing linear differential equations. These operators are linear combinations of differentiation operators, acting on differentiable functions to produce their derivatives. For a single variable, a linear differential operator is represented as L = a_0(x) + a_1(x)d/dx + ... + a_n(x)d^n/dx^n, with a_0(x), ..., a_n(x) being functions of x. Applying operator L to a function f is denoted by Lf. These operators are linear because they preserve addition and scalar multiplication. They constitute a vector space over the field of real or complex numbers and a free module over the ring of differentiable functions. The kernel of a linear differential operator, which is the set of solutions to the homogeneous equation Ly = 0, forms a vector space.Characteristic Polynomial and Homogeneous Equations
The characteristic polynomial is a key concept in solving homogeneous linear differential equations with constant coefficients. This polynomial is obtained by substituting a trial solution of the form e^(αx) into the differential equation and solving the resulting algebraic equation for α. The roots of the characteristic polynomial, known as eigenvalues, correspond to the fundamental solutions of the differential equation. When the characteristic polynomial has repeated roots, additional solutions are constructed by multiplying the exponential solution by powers of x. For equations with real coefficients, complex conjugate root pairs can be combined using Euler's formula to form real-valued solutions.Non-Homogeneous Equations and the Method of Variation of Constants
Non-homogeneous linear differential equations with constant coefficients take the form y^(n)(x) + a_1y^(n-1)(x) + ... + a_ny(x) = f(x), where f(x) is a given function and y is the unknown function. Several methods exist for solving these equations, such as the method of undetermined coefficients and the annihilator method. The most versatile technique is the method of variation of constants, which seeks a particular solution by allowing the constants in the homogeneous solution to vary as functions. This approach leads to a system of equations that can be solved to determine these functions, resulting in a general solution that combines the particular solution with the homogeneous solution.First-Order Equations and Systems of Linear Differential Equations
First-order linear differential equations are typically solved by integrating both sides of the equation. For non-homogeneous equations, an integrating factor is employed to convert the equation into an exact differential, which is then integrated to find the general solution. Systems of linear differential equations, which involve multiple unknown functions, can be expressed in matrix form. Solving these systems parallels the approach for single first-order equations but includes the complexity of matrix operations. The general solution involves a matrix exponential, and specific initial conditions can be applied to determine a particular solution.Higher-Order Equations and Differential Galois Theory
Higher-order linear differential equations with variable coefficients are generally not solvable by quadrature, as indicated by Picard–Vessiot theory and differential Galois theory. These theories provide criteria to ascertain which equations are solvable in closed form and offer methods for solving them when possible, though the procedures can be intricate. An exception is the Cauchy–Euler equation, which, despite variable coefficients, is explicitly solvable. These equations have the form x^n y^(n)(x) + a_(n-1) x^(n-1) y^(n-1)(x) + ... + a_0 y(x) = 0, with constant coefficients a_0, ..., a_(n-1).Holonomic Functions and Their Computational Applications
Holonomic functions, also known as D-finite functions, are solutions to homogeneous linear differential equations with polynomial coefficients. They cover a broad spectrum of functions, including polynomials, exponentials, logarithms, and trigonometric functions. Holonomic functions are closed under operations such as addition, multiplication, differentiation, and definite integration, making these operations computationally efficient. Representing holonomic functions through their defining differential equations and initial conditions enables automated calculus operations, such as differentiation, integration, and limit evaluation. This computational framework is particularly useful for confirming identities and studying the behavior of functions near singularities and at infinity.