Difference between revisions of "Integrating Factor"
| Line 15: | Line 15: | ||
<table class="wiki-content-table"> | <table class="wiki-content-table"> | ||
<tr> | <tr> | ||
| − | <td><strong>Proposition 1:</strong> If <span class="math-inline"><math>\frac{dy}{dt} + p(t) y = g(t)</math></span> is a differential equation, then an integrating factor <span class="math-inline"><math>\mu (t)</math></span> of this equation is given by the formula <span class="math-inline"><math>\mu (t) = e^{\int p(t) \ | + | <td><strong>Proposition 1:</strong> If <span class="math-inline"><math>\frac{dy}{dt} + p(t) y = g(t)</math></span> is a differential equation, then an integrating factor <span class="math-inline"><math>\mu (t)</math></span> of this equation is given by the formula <span class="math-inline"><math>\mu (t) = e^{\int p(t) \; dt}</math></span>.</td> |
</tr> | </tr> | ||
</table> | </table> | ||
| Line 21: | Line 21: | ||
<li><strong>Proof:</strong> We want to find <span class="math-inline"><math>\mu (t)</math></span> such that <span class="math-inline"><math>\mu ' (t) = \mu (t) p(t)</math></span>. We can rewrite this equation as as <span class="math-inline"><math>\frac{\mu '(t)}{\mu (t)} = p(t)</math></span> and then:</li> | <li><strong>Proof:</strong> We want to find <span class="math-inline"><math>\mu (t)</math></span> such that <span class="math-inline"><math>\mu ' (t) = \mu (t) p(t)</math></span>. We can rewrite this equation as as <span class="math-inline"><math>\frac{\mu '(t)}{\mu (t)} = p(t)</math></span> and then:</li> | ||
</ul> | </ul> | ||
| − | <div style="text-align: center;"><math>\begin{align} \quad \frac{ \mu '(t)}{ \mu (t)} = p(t) \\ \quad \frac{d}{dt} \ln \mid \mu (t) \mid = p(t) \\ \quad \int \frac{d}{dt} \ln \mid \mu (t) \mid \ | + | <div style="text-align: center;"><math>\begin{align} \quad \frac{ \mu '(t)}{ \mu (t)} = p(t) \\ \quad \frac{d}{dt} \ln \mid \mu (t) \mid = p(t) \\ \quad \int \frac{d}{dt} \ln \mid \mu (t) \mid \; dt = \int p(t) \; dt \\ \quad \ln \mid \mu (t) \mid = \int p(t) \; dt \\ \quad \mu (t) = \pm e^{\int p(t) \; dt} \end{align}</math></div> |
<ul> | <ul> | ||
| − | <li>Since we only need one integrating factor to solve differential equations in the form <span class="math-inline"><math>\frac{dy}{dt} + p(t) y = g(t)</math></span>, we can more generally note that <span class="math-inline"><math>\mu (t) = e^{\int p(t) \ | + | <li>Since we only need one integrating factor to solve differential equations in the form <span class="math-inline"><math>\frac{dy}{dt} + p(t) y = g(t)</math></span>, we can more generally note that <span class="math-inline"><math>\mu (t) = e^{\int p(t) \; dt}</math></span> is an integrating factor of this differential equation. <span class="math-inline"><math>\blacksquare</math></span></li> |
</ul> | </ul> | ||
| − | <p><em>Notice that from proposition 1 that integrating factors <span class="math-inline"><math>\mu (t)</math></span> are not unique. In fact, there are infinitely many integrating factors. This can be see when evaluating the indefinite integral in <span class="math-inline"><math>\mu (t) = e^{\int p(t) \ | + | <p><em>Notice that from proposition 1 that integrating factors <span class="math-inline"><math>\mu (t)</math></span> are not unique. In fact, there are infinitely many integrating factors. This can be see when evaluating the indefinite integral in <span class="math-inline"><math>\mu (t) = e^{\int p(t) \; dt}</math></span> which will result in getting <span class="math-inline"><math>\mu (t) = e^{P(t) + C}</math></span> where <span class="math-inline"><math>P</math></span> is any antiderivative if <span class="math-inline"><math>p</math></span> and where <span class="math-inline"><math>C</math></span> is a constant. We will always use the simplest integrating factor in solving differential equations of this type.</em></p> |
<p>Let's now look at some examples of applying the method of integrating factors.</p> | <p>Let's now look at some examples of applying the method of integrating factors.</p> | ||
===Example 1=== | ===Example 1=== | ||
<p><strong>Find all solutions to the differential equation <span class="math-inline"><math>\frac{dy}{dt} + \frac{2y}{t} = \frac{\sin t}{t^2}</math></span>.</strong></p> | <p><strong>Find all solutions to the differential equation <span class="math-inline"><math>\frac{dy}{dt} + \frac{2y}{t} = \frac{\sin t}{t^2}</math></span>.</strong></p> | ||
<p>We first notice that our differential equation is in the appropriate form <span class="math-inline"><math>\frac{dy}{dt} + p(t) y = g(t)</math></span> where <span class="math-inline"><math>p(t) = \frac{2}{t}</math></span> and <span class="math-inline"><math>g(t) = \frac{\sin t}{t^2}</math></span>. We compute our integrating factor as:</p> | <p>We first notice that our differential equation is in the appropriate form <span class="math-inline"><math>\frac{dy}{dt} + p(t) y = g(t)</math></span> where <span class="math-inline"><math>p(t) = \frac{2}{t}</math></span> and <span class="math-inline"><math>g(t) = \frac{\sin t}{t^2}</math></span>. We compute our integrating factor as:</p> | ||
| − | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) = e^{\int p(t) \ | + | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) = e^{\int p(t) \; dt} = e^{\int \frac{2}{t} \; dt} = e^{2 \ln (t)} = e^{\ln (t^2)} = t^2 \end{align}</math></div> |
<p>Thus we have that for <span class="math-inline"><math>C</math></span> as a constant:</p> | <p>Thus we have that for <span class="math-inline"><math>C</math></span> as a constant:</p> | ||
| − | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) \frac{dy}{dt} + \mu (t) p(t) y = \mu (t) g(t) \\ \quad t^2 \frac{dy}{dt} + 2t y = \sin t \\ \quad \frac{d}{dt} \left ( t^2 y \right ) = \sin t \\ \quad \int \frac{d}{dt} \left ( t^2 y \right ) \ | + | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) \frac{dy}{dt} + \mu (t) p(t) y = \mu (t) g(t) \\ \quad t^2 \frac{dy}{dt} + 2t y = \sin t \\ \quad \frac{d}{dt} \left ( t^2 y \right ) = \sin t \\ \quad \int \frac{d}{dt} \left ( t^2 y \right ) \; dt = \int \sin t \; dt \\ \quad t^2 y = -\cos t + C \\ \quad y = \frac{-\cos t}{t^2} + \frac{C}{t^2} \end{align}</math></div> |
===Example 2=== | ===Example 2=== | ||
<p><strong>Find all solutions to the differential equation <span class="math-inline"><math>\frac{dy}{dt} - \frac{y}{t} + te^{-t} = 0</math></span>.</strong></p> | <p><strong>Find all solutions to the differential equation <span class="math-inline"><math>\frac{dy}{dt} - \frac{y}{t} + te^{-t} = 0</math></span>.</strong></p> | ||
<p>We first rewrite our differential equation as <span class="math-inline"><math>\frac{dy}{dt} - \frac{y}{t} = - te^{-t}</math></span>. We note that in this form we have <span class="math-inline"><math>p(t) = - \frac{1}{t}</math></span> and <span class="math-inline"><math>g(t) = -t e^{-t}</math></span>. We now find an integrating factor:</p> | <p>We first rewrite our differential equation as <span class="math-inline"><math>\frac{dy}{dt} - \frac{y}{t} = - te^{-t}</math></span>. We note that in this form we have <span class="math-inline"><math>p(t) = - \frac{1}{t}</math></span> and <span class="math-inline"><math>g(t) = -t e^{-t}</math></span>. We now find an integrating factor:</p> | ||
| − | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) = e^{\int p(t) \ | + | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) = e^{\int p(t) \; dt} = e^{\int -\frac{1}{t} \; dt} = e^{-\ln t} = e^{\ln \left ( \frac{1}{t} \right)} = \frac{1}{t} \end{align}</math></div> |
<p>Thus we have that for <span class="math-inline"><math>C</math></span> as a constant:</p> | <p>Thus we have that for <span class="math-inline"><math>C</math></span> as a constant:</p> | ||
| − | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) \frac{dy}{dt} + \mu (t) p(t) y = \mu (t) g(t) \\ \quad \frac{1}{t} \frac{dy}{dt} - \frac{1}{t^2} y = -e^{-t} \\ \quad \frac{d}{dt} \left ( \frac{y}{t} \right ) = -e^{-t} \\ \quad \int \frac{d}{dt} \left ( \frac{y}{t} \right ) \ | + | <div style="text-align: center;"><math>\begin{align} \quad \mu (t) \frac{dy}{dt} + \mu (t) p(t) y = \mu (t) g(t) \\ \quad \frac{1}{t} \frac{dy}{dt} - \frac{1}{t^2} y = -e^{-t} \\ \quad \frac{d}{dt} \left ( \frac{y}{t} \right ) = -e^{-t} \\ \quad \int \frac{d}{dt} \left ( \frac{y}{t} \right ) \; dt = -\int e^{-t} \; dt \\ \quad \frac{y}{t} = e^{-t} + C \\ \quad y = te^{-t} + tC \end{align}</math></div> |
Revision as of 14:02, 16 November 2021
The Method of Integrating Factors
Let Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p} and be functions of and consider the following first order differential equation:
If we multiply the both sides of the equation above by the function we get that:
If we can guarantee that , then notice that by applying the product rule for differentiation that we get: and substituting and we get that which is exactly the lefthand side of the equation above. Thus we get that:
The above differential equation can be solved by integrating both sides of the equation with respect to and isolating . The question now arises on how we can find such a function .
| Definition: If is a first order differential equation, then is called an Integrating Factor if for we have that . |
The following proposition will give us a formula for obtaining the integrating factor for differential equations in the form .
| Proposition 1: If is a differential equation, then an integrating factor of this equation is given by the formula Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu (t) = e^{\int p(t) \; dt}} . |
- Proof: We want to find Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu (t)} such that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu ' (t) = \mu (t) p(t)} . We can rewrite this equation as as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\mu '(t)}{\mu (t)} = p(t)} and then:
- Since we only need one integrating factor to solve differential equations in the form Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{dy}{dt} + p(t) y = g(t)} , we can more generally note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu (t) = e^{\int p(t) \; dt}} is an integrating factor of this differential equation. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \blacksquare}
Notice that from proposition 1 that integrating factors Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu (t)} are not unique. In fact, there are infinitely many integrating factors. This can be see when evaluating the indefinite integral in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu (t) = e^{\int p(t) \; dt}} which will result in getting Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mu (t) = e^{P(t) + C}} where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle P} is any antiderivative if Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p} and where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle C} is a constant. We will always use the simplest integrating factor in solving differential equations of this type.
Let's now look at some examples of applying the method of integrating factors.
Example 1
Find all solutions to the differential equation Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{dy}{dt} + \frac{2y}{t} = \frac{\sin t}{t^2}} .
We first notice that our differential equation is in the appropriate form Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{dy}{dt} + p(t) y = g(t)} where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p(t) = \frac{2}{t}} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g(t) = \frac{\sin t}{t^2}} . We compute our integrating factor as:
Thus we have that for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle C} as a constant:
Example 2
Find all solutions to the differential equation Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{dy}{dt} - \frac{y}{t} + te^{-t} = 0} .
We first rewrite our differential equation as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{dy}{dt} - \frac{y}{t} = - te^{-t}} . We note that in this form we have Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p(t) = - \frac{1}{t}} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g(t) = -t e^{-t}} . We now find an integrating factor:
Thus we have that for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle C} as a constant: