Menu
Home
People
Places
Arts
History
Plants & Animals
Science
Life & Culture
Technology
Reference.org
Method of undetermined coefficients
open-in-new
<h2 id="description-of-the-method">Description of the method</h2> <p>Consider a linear non-homogeneous ordinary differential equation of the form </p> ∑ i = 0 n c i y ( i ) + y ( n + 1 ) = g ( x ) {\displaystyle \sum _{i=0}^{n}c_{i}y^{(i)}+y^{(n+1)}=g(x)} where y ( i ) {\displaystyle y^{(i)}} denotes the i-th derivative of y {\displaystyle y} , and c i {\displaystyle c_{i}} denotes a function of x {\displaystyle x} . <p>The method of undetermined coefficients provides a straightforward method of obtaining the solution to this ODE when two criteria are met:<a class="footnote-ref" id="fnref:2" href="#fn:2"><sup>2</sup></a> </p> <ol><li> c i {\displaystyle c_{i}} are constants.</li> <li><i>g</i>(<i>x</i>) is a constant, a polynomial function, exponential function e α x {\displaystyle e^{\alpha x}} , sine or cosine functions sin β x {\displaystyle \sin {\beta x}} or cos β x {\displaystyle \cos {\beta x}} , or finite sums and products of these functions ( α {\displaystyle {\alpha }} , β {\displaystyle {\beta }} constants).</li></ol> <p>The method consists of finding the general <a href="/facts/Homogeneous_differential_equation/UmTJCNEh">homogeneous</a> solution y c {\displaystyle y_{c}} for the complementary linear <a href="/facts/Homogeneous_differential_equation/UmTJCNEh">homogeneous differential equation</a> </p> ∑ i = 0 n c i y ( i ) + y ( n + 1 ) = 0 , {\displaystyle \sum _{i=0}^{n}c_{i}y^{(i)}+y^{(n+1)}=0,} <p>and a particular integral y p {\displaystyle y_{p}} of the linear non-homogeneous ordinary differential equation based on g ( x ) {\displaystyle g(x)} . Then the general solution y {\displaystyle y} to the linear non-homogeneous ordinary differential equation would be </p> y = y c + y p . {\displaystyle y=y_{c}+y_{p}.} <a class="footnote-ref" id="fnref:3" href="#fn:3"><sup>3</sup></a> <p>If g ( x ) {\displaystyle g(x)} consists of the sum of two functions h ( x ) + w ( x ) {\displaystyle h(x)+w(x)} and we say that y p 1 {\displaystyle y_{p_{1}}} is the solution based on h ( x ) {\displaystyle h(x)} and y p 2 {\displaystyle y_{p_{2}}} the solution based on w ( x ) {\displaystyle w(x)} . Then, using a <a href="/facts/Superposition_principle/GKGs9C92">superposition principle</a>, we can say that the particular integral y p {\displaystyle y_{p}} is<a class="footnote-ref" id="fnref:4" href="#fn:4"><sup>4</sup></a> </p> y p = y p 1 + y p 2 . {\displaystyle y_{p}=y_{p_{1}}+y_{p_{2}}.} <h2 id="typical-forms-of-the-particular-integral">Typical forms of the particular integral</h2> <p>In order to find the particular integral, we need to 'guess' its form, with some coefficients left as variables to be solved for. This takes the form of the first derivative of the complementary function. Below is a table of some typical functions and the solution to guess for them. </p> <table><tbody><tr><th>Function of <i>x</i></th><th>Form for <i>y</i></th></tr><tr><td> k e a x {\displaystyle ke^{ax}\!} </td><td> C e a x {\displaystyle Ce^{ax}\!} </td></tr><tr><td> k x n , n = 0 , 1 , 2 , … {\displaystyle kx^{n},\;n=0,1,2,\ldots \!} </td><td><p> ∑ i = 0 n K i x i {\displaystyle \sum _{i=0}^{n}K_{i}x^{i}\!} </p></td></tr><tr><td> k cos ( a x ) or k sin ( a x ) {\displaystyle k\cos(ax){\text{ or }}k\sin(ax)\!} </td><td><p> K cos ( a x ) + M sin ( a x ) {\displaystyle K\cos(ax)+M\sin(ax)\!} </p></td></tr><tr><td> k e a x cos ( b x ) or k e a x sin ( b x ) {\displaystyle ke^{ax}\cos(bx){\text{ or }}ke^{ax}\sin(bx)\!} </td><td><p> e a x ( K cos ( b x ) + M sin ( b x ) ) {\displaystyle e^{ax}(K\cos(bx)+M\sin(bx))\!} </p></td></tr><tr><td> ( ∑ i = 0 n k i x i ) cos ( b x ) or ( ∑ i = 0 n k i x i ) sin ( b x ) {\displaystyle \left(\sum _{i=0}^{n}k_{i}x^{i}\right)\cos(bx){\text{ or }}\ \left(\sum _{i=0}^{n}k_{i}x^{i}\right)\sin(bx)\!} </td><td><p> ( ∑ i = 0 n Q i x i ) cos ( b x ) + ( ∑ i = 0 n R i x i ) sin ( b x ) {\displaystyle \left(\sum _{i=0}^{n}Q_{i}x^{i}\right)\cos(bx)+\left(\sum _{i=0}^{n}R_{i}x^{i}\right)\sin(bx)} </p></td></tr><tr><td> ( ∑ i = 0 n k i x i ) e a x cos ( b x ) or ( ∑ i = 0 n k i x i ) e a x sin ( b x ) {\displaystyle \left(\sum _{i=0}^{n}k_{i}x^{i}\right)e^{ax}\cos(bx){\text{ or }}\left(\sum _{i=0}^{n}k_{i}x^{i}\right)e^{ax}\sin(bx)\!} </td><td><p> e a x ( ( ∑ i = 0 n Q i x i ) cos ( b x ) + ( ∑ i = 0 n R i x i ) sin ( b x ) ) {\displaystyle e^{ax}\left(\left(\sum _{i=0}^{n}Q_{i}x^{i}\right)\cos(bx)+\left(\sum _{i=0}^{n}R_{i}x^{i}\right)\sin(bx)\right)} </p></td></tr></tbody></table> <p>If a term in the above particular integral for <i>y</i> appears in the homogeneous solution, it is necessary to multiply by a sufficiently large power of <i>x</i> in order to make the solution independent. If the function of <i>x</i> is a sum of terms in the above table, the particular integral can be guessed using a sum of the corresponding terms for <i>y</i>.<a class="footnote-ref" id="fnref:5" href="#fn:5"><sup>5</sup></a> </p> <h2 id="examples">Examples</h2> <h3>Example 1</h3> <p>Find a particular integral of the equation </p> y ″ + y = t cos t . {\displaystyle y''+y=t\cos t.} <p>The right side <i>t</i> cos <i>t</i> has the form </p> P n e α t cos β t {\displaystyle P_{n}e^{\alpha t}\cos {\beta t}} <p>with <i>n</i> = 2, <i>α</i> = 0, and <i>β</i> = 1. </p><p>Since <i>α</i> + <i>iβ</i> = <i>i</i> is <i>a simple root</i> of the characteristic equation </p> λ 2 + 1 = 0 {\displaystyle \lambda ^{2}+1=0} <p>we should try a particular integral of the form </p> y p = t [ F 1 ( t ) e α t cos β t + G 1 ( t ) e α t sin β t ] = t [ F 1 ( t ) cos t + G 1 ( t ) sin t ] = t [ ( A 0 t + A 1 ) cos t + ( B 0 t + B 1 ) sin t ] = ( A 0 t 2 + A 1 t ) cos t + ( B 0 t 2 + B 1 t ) sin t . {\displaystyle {\begin{aligned}y_{p}&=t\left[F_{1}(t)e^{\alpha t}\cos {\beta t}+G_{1}(t)e^{\alpha t}\sin {\beta t}\right]\\&=t\left[F_{1}(t)\cos t+G_{1}(t)\sin t\right]\\&=t\left[\left(A_{0}t+A_{1}\right)\cos t+\left(B_{0}t+B_{1}\right)\sin t\right]\\&=\left(A_{0}t^{2}+A_{1}t\right)\cos t+\left(B_{0}t^{2}+B_{1}t\right)\sin t.\end{aligned}}} <p>Substituting <i>y</i><i>p</i> into the differential equation, we have the identity </p> t cos t = y p ″ + y p = [ ( A 0 t 2 + A 1 t ) cos t + ( B 0 t 2 + B 1 t ) sin t ] ″ + [ ( A 0 t 2 + A 1 t ) cos t + ( B 0 t 2 + B 1 t ) sin t ] = [ 2 A 0 cos t + 2 ( 2 A 0 t + A 1 ) ( − sin t ) + ( A 0 t 2 + A 1 t ) ( − cos t ) + 2 B 0 sin t + 2 ( 2 B 0 t + B 1 ) cos t + ( B 0 t 2 + B 1 t ) ( − sin t ) ] + [ ( A 0 t 2 + A 1 t ) cos t + ( B 0 t 2 + B 1 t ) sin t ] = [ 4 B 0 t + ( 2 A 0 + 2 B 1 ) ] cos t + [ − 4 A 0 t + ( − 2 A 1 + 2 B 0 ) ] sin t . {\displaystyle {\begin{aligned}t\cos t&=y_{p}''+y_{p}\\&=\left[\left(A_{0}t^{2}+A_{1}t\right)\cos t+\left(B_{0}t^{2}+B_{1}t\right)\sin t\right]''+\left[\left(A_{0}t^{2}+A_{1}t\right)\cos t+\left(B_{0}t^{2}+B_{1}t\right)\sin t\right]\\&=\left[2A_{0}\cos t+2\left(2A_{0}t+A_{1}\right)(-\sin t)+\left(A_{0}t^{2}+A_{1}t\right)(-\cos t)+2B_{0}\sin t+2\left(2B_{0}t+B_{1}\right)\cos t+\left(B_{0}t^{2}+B_{1}t\right)(-\sin t)\right]\\&\qquad +\left[\left(A_{0}t^{2}+A_{1}t\right)\cos t+\left(B_{0}t^{2}+B_{1}t\right)\sin t\right]\\&=[4B_{0}t+(2A_{0}+2B_{1})]\cos t+[-4A_{0}t+(-2A_{1}+2B_{0})]\sin t.\end{aligned}}} <p>Comparing both sides, we have </p> { 1 = 4 B 0 0 = 2 A 0 + 2 B 1 0 = − 4 A 0 0 = − 2 A 1 + 2 B 0 {\displaystyle {\begin{cases}1=4B_{0}\\0=2A_{0}+2B_{1}\\0=-4A_{0}\\0=-2A_{1}+2B_{0}\end{cases}}} <p>which has the solution </p> A 0 = 0 , A 1 = B 0 = 1 4 , B 1 = 0. {\displaystyle A_{0}=0,\quad A_{1}=B_{0}={\frac {1}{4}},\quad B_{1}=0.} <p>We then have a particular integral </p> y p = 1 4 t cos t + 1 4 t 2 sin t . {\displaystyle y_{p}={\frac {1}{4}}t\cos t+{\frac {1}{4}}t^{2}\sin t.} <h3>Example 2</h3> <p>Consider the following linear nonhomogeneous differential equation: </p> d y d x = y + e x . {\displaystyle {\frac {dy}{dx}}=y+e^{x}.} <p>This is like the first example above, except that the nonhomogeneous part ( e x {\displaystyle e^{x}} ) is <i>not</i> linearly independent to the general solution of the homogeneous part ( c 1 e x {\displaystyle c_{1}e^{x}} ); as a result, we have to multiply our guess by a sufficiently large power of <i>x</i> to make it linearly independent. </p><p>Here our guess becomes: </p> y p = A x e x . {\displaystyle y_{p}=Axe^{x}.} <p>By substituting this function and its derivative into the differential equation, one can solve for <i>A</i>: </p> d d x ( A x e x ) = A x e x + e x {\displaystyle {\frac {d}{dx}}\left(Axe^{x}\right)=Axe^{x}+e^{x}} A x e x + A e x = A x e x + e x {\displaystyle Axe^{x}+Ae^{x}=Axe^{x}+e^{x}} A = 1. {\displaystyle A=1.} <p>So, the general solution to this differential equation is: </p> y = c 1 e x + x e x . {\displaystyle y=c_{1}e^{x}+xe^{x}.} <h3>Example 3</h3> <p>Find the general solution of the equation: </p> d y d t = t 2 − y {\displaystyle {\frac {dy}{dt}}=t^{2}-y} <p> t 2 {\displaystyle t^{2}} is a polynomial of degree 2, so we look for a solution using the same form, </p> y p = A t 2 + B t + C , {\displaystyle y_{p}=At^{2}+Bt+C,} <p>Plugging this particular function into the original equation yields, </p> 2 A t + B = t 2 − ( A t 2 + B t + C ) , {\displaystyle 2At+B=t^{2}-(At^{2}+Bt+C),} 2 A t + B = ( 1 − A ) t 2 − B t − C , {\displaystyle 2At+B=(1-A)t^{2}-Bt-C,} ( A − 1 ) t 2 + ( 2 A + B ) t + ( B + C ) = 0. {\displaystyle (A-1)t^{2}+(2A+B)t+(B+C)=0.} <p>which gives: </p> A − 1 = 0 , 2 A + B = 0 , B + C = 0. {\displaystyle A-1=0,\quad 2A+B=0,\quad B+C=0.} <p>Solving for constants we get: </p> y p = t 2 − 2 t + 2 {\displaystyle y_{p}=t^{2}-2t+2} <p>To solve for the general solution, </p> y = y p + y c {\displaystyle y=y_{p}+y_{c}} <p>where y c {\displaystyle y_{c}} is the homogeneous solution y c = c 1 e − t {\displaystyle y_{c}=c_{1}e^{-t}} , therefore, the general solution is: </p> y = t 2 − 2 t + 2 + c 1 e − t {\displaystyle y=t^{2}-2t+2+c_{1}e^{-t}} <ul><li>Boyce, W. E.; DiPrima, R. C. (1986). <a href="https://archive.org/details/elementarydiffer4th00boyc"><i>Elementary Differential Equations and Boundary Value Problems</i></a> (4th ed.). John Wiley & Sons. <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 0-471-83824-1.</li> <li>Riley, K. F.; Bence, S. J. (2010). <a href="https://archive.org/details/mathematicalmeth00rile"><i>Mathematical Methods for Physics and Engineering</i></a>. Cambridge University Press. <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-521-86153-3.</li> <li>Tenenbaum, Morris; Pollard, Harry (1985). <a href="https://archive.org/details/ordinarydifferen00tene_0"><i>Ordinary Differential Equations</i></a>. Dover. <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-486-64940-5.</li> <li>de Oliveira, O. R. B. (2013). "A formula substituting the undetermined coefficients and the annihilator methods". <i>Int. J. Math. Educ. Sci. Technol</i>. 44 (3): 462–468. <a href="/facts/ArXiv_(identifier)/H6EtgnBe">arXiv</a>:<a href="https://arxiv.org/abs/1110.4425">1110.4425</a>. <a href="/facts/Bibcode_(identifier)/9HtdQSGB">Bibcode</a>:<a href="https://ui.adsabs.harvard.edu/abs/2013IJMES..44..462R">2013IJMES..44..462R</a>. <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1080%2F0020739X.2012.714496">10.1080/0020739X.2012.714496</a>. <a href="/facts/S2CID_(identifier)/ldJsHa2Y">S2CID</a> <a href="https://api.semanticscholar.org/CorpusID:55834468">55834468</a>.</li></ul> <h2 id="references">References</h2> <ol> <li id="fn:1"><p>Ralph P. Grimaldi (2000). "Nonhomogeneous Recurrence Relations". Section 3.3.3 of Handbook of Discrete and Combinatorial Mathematics. Kenneth H. Rosen, ed. CRC Press. ISBN 0-8493-0149-1. <a href="/wiki/ISBN_(identifier)" target="_blank">/wiki/ISBN_(identifier)</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li> <li id="fn:2"><p>Zill, Dennis G., Warren S. Wright (2014). Advanced Engineering Mathematics. Jones and Bartlett. p. 125. ISBN 978-1-4496-7977-4.{{cite book}}: CS1 maint: multiple names: authors list (link) <a href="978-1-4496-7977-4" target="_blank">978-1-4496-7977-4</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li> <li id="fn:3"><p>Dennis G. Zill (14 May 2008). A First Course in Differential Equations. Cengage Learning. ISBN 978-0-495-10824-5. <a href="978-0-495-10824-5" target="_blank">978-0-495-10824-5</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li> <li id="fn:4"><p>Dennis G. Zill (14 May 2008). A First Course in Differential Equations. Cengage Learning. ISBN 978-0-495-10824-5. <a href="978-0-495-10824-5" target="_blank">978-0-495-10824-5</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></p></li> <li id="fn:5"><p>Ralph P. Grimaldi (2000). "Nonhomogeneous Recurrence Relations". Section 3.3.3 of Handbook of Discrete and Combinatorial Mathematics. Kenneth H. Rosen, ed. CRC Press. ISBN 0-8493-0149-1. <a href="/wiki/ISBN_(identifier)" target="_blank">/wiki/ISBN_(identifier)</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></p></li> </ol>