{
 "id": "epmkjw1c38",
 "slug": "variation-of-parameters",
 "title": "Variation of parameters",
 "updated": "2026-10-15",
 "topic_path": [
  {
   "id": "physical",
   "label": "Physical world and mathematics",
   "api_url": "https://www.edgepedia.org/api/v1/topics/physical"
  },
  {
   "id": "physical.math",
   "label": "Mathematics and statistics",
   "api_url": "https://www.edgepedia.org/api/v1/topics/physical.math"
  },
  {
   "id": "physical.math.analysis",
   "label": "Analysis and mathematical models",
   "api_url": "https://www.edgepedia.org/api/v1/topics/physical.math.analysis"
  },
  {
   "id": "physical.math.analysis.ordinary_differential_equations",
   "label": "Ordinary differential equations",
   "api_url": "https://www.edgepedia.org/api/v1/topics/physical.math.analysis.ordinary_differential_equations"
  }
 ],
 "geo": [
  {
   "id": "geo.nongeo.t1000.physical",
   "label": "Non-geographic · 1000 to 1799: Physical world and mathematics",
   "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1000.physical",
   "path": [
    {
     "id": "geo.nongeo",
     "label": "Non-geographic",
     "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo"
    },
    {
     "id": "geo.nongeo.t1000",
     "label": "Non-geographic · 1000 to 1799",
     "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1000"
    },
    {
     "id": "geo.nongeo.t1000.physical",
     "label": "Physical world and mathematics",
     "api_url": "https://www.edgepedia.org/api/v1/geo/geo.nongeo.t1000.physical"
    }
   ]
  }
 ],
 "excerpt": "Variation of parameters is a method for solving nonhomogeneous linear differential equations, replacing the constants in the homogeneous solution with functions chosen to satisfy the full equation.",
 "snippet": "Variation of parameters is a method for solving nonhomogeneous linear differential equations, replacing the constants in the homogeneous solution with functions chosen to satisfy the full equation.",
 "node": "physical.math.analysis.ordinary_differential_equations",
 "markdown": "# Variation of parameters\n\nVariation of parameters is a method for solving a nonhomogeneous linear differential equation by replacing the arbitrary constants in the general solution of the associated homogeneous equation with unknown functions, chosen so that the result also satisfies the nonhomogeneous equation.\n\n| Key fact | Detail |\n|---|---|\n| What it produces | A particular solution of a nonhomogeneous linear ODE or system, formed from a fundamental set of homogeneous solutions.[3](https://terpconnect.umd.edu/~lvrmr/2008-2009-S/Classes/MATH246/NOTES/HigherLin4.pdf) |\n| Second-order formula | with \\( W \\) the Wronskian of \\( y_{1}, y_{2} \\).[5](https://www.math.utah.edu/~gustafso/2250variation-of-parameters.pdf) |\n| Pinning condition | The constraint fixes the unknown functions; it is the envelope condition \\( \\partial Y/\\partial \\xi = 0 \\) at tangency of a family of trajectories.[4](https://math.berkeley.edu/~arash/54/notes/n4_6.pdf)[6](https://projecteuclid.org/journalArticle/Download?urlId=10.35834%2Fmjms%2F1316092232&isResultClick=False) |\n\n| Scope | Works for any continuous (even piecewise continuous) forcing, including \\( \\lvert x\\rvert \\), \\( e^{x^{2}} \\), where undetermined coefficients fails.[5](https://www.math.utah.edu/~gustafso/2250variation-of-parameters.pdf)[8](https://www.cfm.brown.edu/people/dobrush/am33/Mathematica/ch4/varpar.html) |\n| Cost | Requires the complementary solution and two integrals that may have no elementary closed form.[9](https://tutorial.math.lamar.edu/Classes/DE/VariationofParameters.aspx) |\n\n## How it works\n\nThe method starts from the general solution of the homogeneous equation, where \\( y_{1}, y_{2} \\) form a fundamental set. Substituting this ansatz into the nonhomogeneous equation gives one condition on the two unknown derivatives, so a second condition is free to choose. The standard choice is\n\nwhich makes the second-derivative terms cancel and leaves a solvable system. This seemingly arbitrary hypothesis has a geometric meaning: it is the envelope condition at the point of tangency of a one-parameter family of fictitious trajectories.[6](https://projecteuclid.org/journalArticle/Download?urlId=10.35834%2Fmjms%2F1316092232&isResultClick=False)\n\nWith this constraint, the unknowns \\( u_{1}', u_{2}' \\) satisfy a 2 × 2 linear system whose coefficient determinant is exactly the [Wronskian](https://www.edgepedia.org/wronskian). The Wronskian itself obeys Abel's identity, and two homogeneous solutions are independent if and only if \\( W(x) \\neq 0 \\).[5](https://www.math.utah.edu/~gustafso/2250variation-of-parameters.pdf)\n\n## How it is done\n\nFor the normalized second-order equation with \\( p, q, g \\) continuous, the procedure is:[11](https://courses.math.wichita.edu/math555/ch3/6/VariationOfParameters.html)[12](https://debook.tomcuchta.com/ch2-higherorderodes-8.html)\n\n1. Solve the homogeneous equation.\n2. Put the equation in standard form with leading coefficient 1; this step is necessary because the formula assumes it.\n3. Compute the Wronskian \\( W \\).\n4. Assume \\( u_{1}' \\cdot y_{1} + u_{2}' \\cdot y_{2} = 0 \\).\n5. Substitute and solve by [Cramer's rule](https://www.edgepedia.org/cramers-rule): \\( u_{2}' = \\frac{y_{1} \\cdot g}{W} \\).[11](https://courses.math.wichita.edu/math555/ch3/6/VariationOfParameters.html)\n6. Integrate, absorb any homogeneous terms into \\( y_{h} \\), and apply initial conditions if given.[12](https://debook.tomcuchta.com/ch2-higherorderodes-8.html)\n\nEquivalently, for the normalized second-order equation the particular solution is \\( y_{p}=-y_{1}\\int \\frac{y_{2}g}{W}\\,dx+y_{2}\\int \\frac{y_{1}g}{W}\\,dx \\), with the integrals taken on an interval where the coefficients and the fundamental solutions are defined.\n\nA worked example: \\( y'' + y = \\sec x \\) gives \\( y = c_{1}\\cos x + c_{2}\\sin x + x\\sin x + \\cos(x)\\ln\\lvert\\cos x\\rvert \\) with \\( W = 1 \\).[5](https://www.math.utah.edu/~gustafso/2250variation-of-parameters.pdf)\n\nFor an nth-order equation, the derivatives satisfy \\( u_{i}' = W_{i}/W \\), where \\( W_{i} \\) is the Wronskian with column i replaced by \\( (0, \\ldots, 0, f(x)) \\).[14](https://faculty.fiu.edu/~rukim/MAC2302/ODEu-lect2.pdf) For the first-order system \\( x' = P(t)x + g(t) \\) with fundamental matrix \\( X(t) \\) (\\( \\det X(t) \\neq 0 \\)), a particular solution is \\( x_{p}(t) = X(t)\\int X^{-1}(t)g(t)\\,dt \\), and the general solution is \\( x(t) = X(t)c + x_{p}(t) = \\sum_{i=1}^{n}c_{i}x_{i}(t) + x_{p}(t) \\), where the columns of \\( X \\) are a fundamental set of n solutions.[7](https://www.its.caltech.edu/~zuev/teaching/2012Spring/Math245-Lecture-21.pdf)\n\n## Origin\n\nThe extension of the method to fractional calculus rests on the conformable fractional derivative, introduced by R. Khalil and colleagues in 2014 in the Journal of Computational and Applied Mathematics.[17](https://doi.org/10.1016/j.cam.2014.01.002)\n\n## Variants\n\nFor the inhomogeneous linear system, the Cauchy formula\n\nThe solution of the initial-value problem can likewise be written with a solution operator, a statement known as Duhamel's principle, which reduces to variation of parameters for scalar linear inhomogeneous ODEs.[10](https://ar5iv.labs.arxiv.org/html/1909.04916)\n\nThe particular solution can also be written as a Green-function integral, where \\( G \\) is fixed by initial-value problems for the operator alone; the variation-of-parameters and Green-function methods are essentially equivalent, and \\( G \\) does not depend on which fundamental set was used to build it.[3](https://terpconnect.umd.edu/~lvrmr/2008-2009-S/Classes/MATH246/NOTES/HigherLin4.pdf) For second-order equations the two integrals combine into \\( G(x,\\xi) = \\frac{y_{2}(x)y_{1}(\\xi) - y_{1}(x)y_{2}(\\xi)}{y_{1}(\\xi)y_{2}'(\\xi) - y_{2}(\\xi)y_{1}'(\\xi)} \\), known as a [Green's function](https://www.edgepedia.org/greens-function).[8](https://www.cfm.brown.edu/people/dobrush/am33/Mathematica/ch4/varpar.html) When the operator has constant coefficients, \\( G(t,s) = g(t-s) \\), so the particular solution becomes a convolution; for \\( y'' + y = f(t) \\) the general solution is \\( y(t) = c_{1}\\cos t + c_{2}\\sin t + \\int_{0}^{t}f(s)\\sin(t-s)\\,ds \\).[3](https://terpconnect.umd.edu/~lvrmr/2008-2009-S/Classes/MATH246/NOTES/HigherLin4.pdf)[12](https://debook.tomcuchta.com/ch2-higherorderodes-8.html)\n\nThe idea also extends beyond linear equations: the variation-of-constants formula holds as an integral equation relating the perturbed and unperturbed systems.[1](https://encyclopediaofmath.org/wiki/Variation_of_constants) In fractional calculus, the method has been extended to nonhomogeneous linear fractional differential equations using the conformable fractional derivative.[16](https://isr-publications.com/jmcs/articles-1032-variation-of-parameters-for-local-fractional-nonhomogenous-lineardifferential-equations)[17](https://doi.org/10.1016/j.cam.2014.01.002) That extension defines an \\( \\alpha \\)-Wronskian \\( W_{\\alpha}[y_{1}, y_{2}] \\) built from conformable derivatives of order \\( \\alpha \\), gives \\( c_{1}' = -f(x)y_{2}/W_{\\alpha} \\) and \\( c_{2}' = f(x)y_{1}/W_{\\alpha} \\) for the second-order case, and yields \\( c_{m}^{(\\alpha)}(x) = f(x)W_{\\alpha}^{m}/W_{\\alpha} \\) by Cramer's rule in the nth-order case, with a particular-solution kernel weighted by \\( t^{1-\\alpha}\\,dt \\).[16](https://isr-publications.com/jmcs/articles-1032-variation-of-parameters-for-local-fractional-nonhomogenous-lineardifferential-equations)\n\n## Applications\n\nThe method's main application is solving linear nonhomogeneous ODEs whose forcing terms defeat other techniques. It handles inputs such as \\( \\lvert x\\rvert \\) and \\( e^{x^{2}} \\), where undetermined coefficients fails, assuming continuity of the coefficients and \\( a(x) \\neq 0 \\).[5](https://www.math.utah.edu/~gustafso/2250variation-of-parameters.pdf) Through the Green-function and Duhamel forms it underlies the input-output description of linear systems, where constant-coefficient operators give convolution integrals.[3](https://terpconnect.umd.edu/~lvrmr/2008-2009-S/Classes/MATH246/NOTES/HigherLin4.pdf)\n\n## Limitations and alternatives\n\nThe method has two practical disadvantages: the complementary solution is absolutely required, and the method requires evaluating integrals for which no closed form may exist.[9](https://tutorial.math.lamar.edu/Classes/DE/VariationofParameters.aspx) Finding a fundamental set is itself difficult when \\( P(t) \\) is not constant, and the computations are often tedious.[7](https://www.its.caltech.edu/~zuev/teaching/2012Spring/Math245-Lecture-21.pdf) It needs many integrations and explicit knowledge of a fundamental set of homogeneous solutions.[8](https://www.cfm.brown.edu/people/dobrush/am33/Mathematica/ch4/varpar.html)\n\nThe [method of undetermined coefficients](https://www.edgepedia.org/method-of-undetermined-coefficients) fails when the forcing is outside its trial-function class, for example \\( t^{-1} \\) or \\( \\ln t \\) on suitable intervals, and variation of parameters is then one possible alternative.[13](https://math.libretexts.org/Bookshelves/Analysis/Supplemental_Modules_(Analysis)/Ordinary_Differential_Equations/3%3A_Second_Order_Linear_Differential_Equations/3.5%3A_Variation_of_Parameters)\n\n## References\n\n---\n*Topic: Encyclopedia › Physical world and mathematics › Mathematics and statistics › Analysis and mathematical models › Ordinary differential equations*\n\n*Initially written Sep 29, 2026 · Edited: Oct 15, 2026*\n\n*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*\n\nLicense: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license\n",
 "same_as": [],
 "url": "https://www.edgepedia.org/variation-of-parameters",
 "markdown_url": "https://www.edgepedia.org/variation-of-parameters.md",
 "license": {
  "name": "Edgepedia Community License 1.0",
  "url": "https://www.edgepedia.org/edgepedia/license",
  "summary": "Free with credit, commercial use included. AI training is open to everyone. For other uses, organizations over USD 100M in revenue or 100M monthly users license separately.",
  "spdx": "LicenseRef-Edgepedia-Community-1.0"
 },
 "credit": "\"Variation of parameters\", Edgepedia (EdgeChat), https://www.edgepedia.org/variation-of-parameters. Edgepedia Community License 1.0.",
 "credit_md": "\"[Variation of parameters](https://www.edgepedia.org/variation-of-parameters)\", Edgepedia (EdgeChat), [https://www.edgepedia.org/variation-of-parameters](https://www.edgepedia.org/variation-of-parameters). [Edgepedia Community License 1.0](https://www.edgepedia.org/edgepedia/license).",
 "credit_html": "\"<a href=\"https://www.edgepedia.org/variation-of-parameters\">Variation of parameters</a>\", Edgepedia (EdgeChat), <a href=\"https://www.edgepedia.org/variation-of-parameters\">https://www.edgepedia.org/variation-of-parameters</a>. <a href=\"https://www.edgepedia.org/edgepedia/license\">Edgepedia Community License 1.0</a>.",
 "speakable": "Variation of parameters is a method for solving nonhomogeneous linear differential equations, replacing the constants in the homogeneous solution with functions chosen to satisfy the full equation."
}