University of Florida/Egm4313/s12.team1/R6

Statement
Find the (smallest) period of $$ \cos n \omega x $$ and $$ \sin n \omega x $$. Show that these functions also have period $$ p $$. Show that the constant $$ a_0 $$ is also a periodic function with period $$ p $$.

Solution
First, we identify the property of periodic functions, as shown on p9-2 of the notes, equation (1)


 * {| style="width:100%" border="0"

$$\displaystyle{f(x+np)=f(x)}$$ (1.0)
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }

So for $$ \cos n \omega x $$, the following determines the smallest period.


 * {| style="width:100%" border="0"

$$\displaystyle{\cos (n\omega x})$$ (1.1) We then add $$ 2 \pi $$ to see a final expression that will be equal to $$\displaystyle{f(x+np)=f(x)} $$
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$\displaystyle{\cos (n\omega x + 2 \pi)}$$ (1.2)
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }

We then pull out a term that leaves an x plus a number (our period) inside the trigonometric function, cosine.
 * {| style="width:100%" border="0"

$$\displaystyle{\cos[n\omega ( x + \frac{2 \pi}{n \omega})]}$$ (1.3) This means the smallest period is $$ \frac{2 \pi}{n \omega} $$.
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }

For $$ \sin n \omega x $$, the following determines the smallest period.


 * {| style="width:100%" border="0"

$$\displaystyle{\sin (n\omega x)}$$ (1.4) We then add $$ 2 \pi $$ to see a final expression that will be equal to $$\displaystyle{f(x+np)=f(x)} $$
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$\displaystyle{\sin (n\omega x + 2 \pi)}$$ (1.5)
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }

We then pull out a term that leaves an x plus a number (our period) inside the trigonometric function, cosine.
 * {| style="width:100%" border="0"

$$\displaystyle{\sin[n\omega ( x + \frac{2 \pi}{n \omega})]}$$ (1.6) This means the smallest period is $$ \frac{2 \pi}{n \omega} $$, just as for its cosine equivalent.
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }

From the notes (page 9-5, (1)), it can be stated that $$ \omega = \frac{2 \pi}{p} $$ where $$ p $$ is the period. Both of these functions have $$ \frac{2 \pi}{n \omega} $$ as their period. If we put in 1 for n, and $$ \frac{2 \pi}{p} $$ for $$\omega $$, the period simplifies to simply $$ \displaystyle p$$.

As for the constant $$ a_0 $$, a constant is, almost by definition, a periodic function. The property of periodic functions as referenced in (1.0) supports this.


 * {| style="width:100%" border="0"

$$\displaystyle{f(x+np)=f(x)}$$ (1.0) For a function of f that is a constant, x is not a factor effecting its value. This means, no matter what is added to x, the function will return the same constant number. This added value can just as well be $$ p $$, and thus, $$ a_o $$ is a periodic function with period $$\displaystyle p $$.
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }

Author
Solved and Typed By -Egm4313.s12.team1.silvestri (talk) 02:19, 7 April 2012 (UTC) Reviewed By - --Egm4313.s12.team1.rosenberg (talk) 01:19, 11 April 2012 (UTC)

Statement
Kreyszig pg 491 pbs. 11,12 Is the given function odd or even or neither even nor odd? Show details of your work. -11)$$ f(x) = x^2 (-1 < x <1), p=2 $$ -12)$$ f(x) = 1- \frac{x^2}{4} (-2 < x <2), p=4 $$ -

Continued)Find the Fourier series expansion for $$ f(x) $$ on p.9-8 as follows: 1. Develop the Fourier series expansion of $$f(\bar{x})$$ Plot $$f(\bar{x})$$ and the truncated Fourier series $$f_n(\bar{x})$$

dsfawf

for n=0,1. Observe the values of $$f(\bar{x})$$ at the points of discontinuities, and the Gibbs phenomenon. Transform the variable so to obtain the Fourier series expansion of $$ f(x) $$. Hint: $$f(\bar{x})$$ is an even function; what would become of the coefficients $$\bar{b}_k$$? 2. Do the same as above, but using $$f(\tilde{x})$$ to obtain the Fourier series expansion of $$ f(x) $$; compare to the result obtained above. Hint: Transform $$f(\tilde{x})$$ into an odd function, and expand this odd function into Fourier series; observe the coefficients of this series.

Solution
'''Part 1: K 2011 p.491 pb. 11'''
 * {| style="width:100%" border="0"

$$  \displaystyle f(x)=x^2, (-1<x<1), p=2 $$     (2.0) The preceding function is an even function. Even functions are functions where $$f(-x)=f(x)\!$$. Odd functions on the other hand are functions where $$f(-x)=-f(x)\!$$. The Fourier series of even and odd functions can be found more easily than the Fourier series of a function that is neither even nor odd. Normal Fourier series contain both sine and cosine terms. Because of the nature of even functions, the coefficient of all sine terms in their series is zero. Likewise, all cosine terms in odd functions have a coefficient of zero. Therefore its Fourier series expansion will be of the form:
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle f(x)=a_0+\sum_{n=1}^{\infty}a_n \cos \frac{n\pi}{L}x $$     (2.1) Where:
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_0=\frac{1}{L} \int_0^Lf(x)dx $$     (2.2) And
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_n=\frac{2}{L} \int_0^Lf(x) \cos \frac{n\pi x}{L}dx $$     (2.3) To get the constant term of the expansion, equation (2.0) will be substituted into equation (2.2) with L=1:
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_0=\frac{1}{1} \int_0^1x^2dx=\left [ \frac{x^3}{3} \right ]^1_0=\frac{1}{3} $$     (2.4) To get the coefficients of the rest of the terms in the expansion, equation (2.0) will be substituted into equation (2.3) with L=1:
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_n=\frac{2}{1} \int_0^1x^2 \cos \frac{n\pi x}{1}dx=2\left [ \frac{(\pi^2n^2x^2-2)\sin (\pi n x)+2\pi n x \cos(\pi nx)}{\pi^3n^3} \right ]^1_0 $$     (2.5) Simplifying, this yields:
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_n=\frac{4}{n^2 \pi^2} \cos n\pi $$     (2.6) Therefore, the Fourier expansion of f(x) is:
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle x^2=\frac{1}{3}+\sum_{n=1}^{\infty}\left (\frac{4}{n^2 \pi^2} \cos n\pi \right ) \cos n\pi x $$ (2.7) '''Part 2: K 2011 p.491 pb. 12'''
 * style="width:95%" |
 * style="width:95%" |
 * 
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle f(x)=1-\frac{x^2}{4}, (-2<x<2), p=4 $$     (2.8) The preceding function is an even function. Therefore, equations (2.1)-(2.3) are applicable. To get the constant term of the expansion, equation (2.8) will be substituted into equation (2.2) with L=1:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_0=\frac{1}{2} \int_0^21-\frac{x^2}{4}dx= \frac{1}{2}\left[ x-\frac{x^3}{12} \right ]^2_0=\frac{4}{3} $$     (2.9) To get the coefficients of the rest of the terms in the expansion, equation (2.0) will be substituted into equation (2.3) with L=1:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_n=\frac{2}{2} \int_0^2\left (1-\frac{x^2}{4} \right ) \cos \frac{n\pi x}{2}dx=\left [ \frac{(8-\pi^2n^2(x^2-4))\sin\frac{n\pi x}{2}-4n\pi x\cos\frac{n\pi x}{2}}{2\pi^3n^3} \right ]^2_0 $$     (2.10) Simplifying, this yields:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle a_n=\frac{-4\cos(n\pi)}{n^2 \pi^2} $$     (2.11) Therefore, the Fourier expansion of f(x) is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle 1-\frac{x^2}{4}=\frac{4}{3}+\sum_{n=1}^{\infty}\left (\frac{-4\cos(n\pi)}{n^2 \pi^2} \right ) \cos \frac{n\pi x}{2} $$     (2.12) Part 3 The Fourier series expansion of the function f(x) in Figure 1 can be found in one of two ways. The first way involves transforming x in the following way:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \bar x=x-\frac54 $$     (2.13) The resulting function is now an even function that alternates between the constant value A and zero and has period 4. Using the same methods used above, the Fourier series of the function can be found. The constant coefficient of the series is found by:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \bar a_0=\frac12 \int_0^2f(x)dx=\frac12 \int_0^1Adx+\frac12 \int_1^20dx=1/2\left [ Ax \right ]_0^1=\frac A2 $$ (2.14) Note that because the function is discontinuous at the midpoint, the integral must be split into two different integrals. The coefficients of the cosine terms of the series can be found by:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \bar a_n=\frac22 \int_0^2f(x)\cos \frac{n\pi x}{2}dx= \int_0^1A\cos \frac{n\pi x}{2}dxdx+\frac12 \int_1^20dx=\left [ \frac{2A}{n\pi}\sin\frac{n\pi x}{2} \right ]_0^1=\frac{2A}{n\pi}\sin \frac{n\pi}{2} $$     (2.15) Therefore, the Fourier series for the transformed function is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle f(\bar x)=\frac A2 + \sum_{n=0}^{\infty} \left ( \frac{2A}{n\pi}\sin\frac{n\pi}{2} \right )\left ( \cos \frac{n\pi\bar x}{2} \right ) $$     (2.16) Through a change in variables, the Fourier series of the original function, f(x), can be found:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle f(x)=\frac A2 + \sum_{n=0}^{\infty} \left ( \frac{2A}{n\pi}\sin\frac{n\pi}{2} \right )\left ( \cos \frac{n\pi(x-1.25)}{2} \right ) $$     (2.17) The following figure shows a plot of $$f(\bar x)\!$$ with its truncated Fourier series for n=0,1.
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Figure 2 Figure 3 shows the same plot with more values of n:

Note that all plots assume that A=1. It should be noted that the accuracy of the expansion increases with an increase in the number of terms in the expansion. Also, there is strong fluctuation at points of discontinuity. These fluctuations are apparent even at high values of n, but they decrease in size and move closer to the points of discontinuity. This is a property of Fourier series known as the Gibbs phenomenon. The other method of developing the Fourier series is to transform the function into an odd function. This can be done by transforming both the independent variable and overall function as follows:
 * {| style="width:100%" border="0"

$$  \displaystyle \tilde x=x-\frac14 $$     (2.18)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle f(\hat x)=f(\tilde x)-\frac A2 $$ (2.19) An odd function has a Fourier series expansion of the form:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle f(x)=\sum_{n=0}^{\infty}b_n\sin \frac{n\pi x}{L} $$     (2.20) Where:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle b_n=\frac2L\int_{0}^{L}f(x)\sin \frac{n\pi x}{L}dx $$     (2.21) Applying this to the transformed function, with L=2:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \hat b_n=\frac22\int_{0}^{2}f(x)\sin \frac{n\pi \hat x}{2}d\hat x=\int_{0}^{2}\frac A2\sin \frac{n\pi \hat x}{2}d\hat x=\left [ -\frac{A}{n\pi}\cos\frac{n\pi\hat x}{2} \right ]^2_0=\frac{A}{n\pi}(1-\cos n\pi) $$     (2.22) The Fourier series is therefore:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle f(\hat x)=\sum_{n=0}^{\infty}\left (\frac{A}{n\pi}(1-\cos n\pi) \right )\left ( \sin\frac{n\pi\hat x}{2} \right ) $$     (2.23)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle f(\tilde x)=\frac A2+\sum_{n=0}^{\infty}\left (\frac{A}{n\pi}(1-\cos n\pi) \right )\left ( \sin\frac{n\pi\tilde x}{2} \right ) $$     (2.24)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle f(x)=\frac A2+\sum_{n=0}^{\infty}\left (\frac{A}{n\pi}(1-\cos n\pi) \right )\left ( \sin\frac{n\pi(x-0.25)}{2} \right ) $$     (2.25) The following figure shows a plot of $$f(\tilde x)\!$$ with its truncated Fourier series for n=0,1.
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Figure 4 Note that all plots assume A=1. The accuracy of the expansion as shown on the plot is expected given the low number of terms used in the expansion. As Figure 3 shows, if the number of terms used to form the plot in Figure 4 is increased, the accuracy of the expansion will greatly increase (but more slowly at points of discontinuity because of the Gibbs phenomenon). This level of accuracy is shown in Figure 5, which contains the expansion to 100 and 200 terms.

Author
Solved and Typed By - Egm4313.s12.team1.armanious (talk) 07:10, 11 April 2012 (UTC) Reviewed By - Egm4313.s12.team1.silvestri (talk) 18:30, 11 April 2012 (UTC)

Statement
Kreyszig pg 491 pbs 15,17 Is the given function even or odd or neither even not odd? Find its Fourier series. Show details of your work.

- Plot the truncated Fourier series for n=2,4,8

Solution
By definition, a function, $$ f(x)\!$$ is even if $$ f(-x) = f(x)\!$$, and its Fourier series reduces to a Fourier cosine series:
 * {| style="width:100%" border="0"

$$\displaystyle{f(x) = a_{0} + \sum_{n=1}^{\infty}a_{n}cos{\frac{n\pi}{L}}x}$$ (3.0)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

A function, $$ f(x)\!$$ is odd if $$ f(-x) = -f(x)\!$$, and its Fourier series reduces to a Fourier sine series:
 * {| style="width:100%" border="0"

$$\displaystyle{f(x) = \sum_{n=1}^{\infty}b_{n}sin{\frac{n\pi}{L}}x}$$ (3.1) -15) By inspection, it is evident that the given function in problem 15 is odd because for every $$ x\!$$ value that is negative, $$ f(x)\!$$ is negative. Therefore, $$ f(-x) = -f(x)\!$$.
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

The following are the defining functions for the given graph:
 * {| style="width:100%" border="0"

$$\displaystyle{f(x) = \begin{cases} -\pi - x & \text{ if } x= -\pi < x < -\frac{\pi}{2}\\ x & \text{ if } x= -\frac{\pi}{2} < x < \frac{\pi}{2}\\ \pi - x & \text{ if } x= \frac{\pi}{2} < x < \pi \end{cases}}$$ (3.2)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Because the function is odd, we use equation (3.1) to develop the Fourier series, where the period is $$2\pi\!$$ and $$ L\!$$ equals $$\pi\!$$:
 * {| style="width:100%" border="0"

$$\displaystyle{f(x) = \sum_{n=1}^{\infty}b_{n}sin{\frac{n\pi}{L}}x = \sum_{n=1}^{\infty}b_{n}sin(nx)}$$
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">

(3.3)
 * }

Now we need to find $$ b_{n}\!$$, which is of the following form for odd functions:
 * {| style="width:100%" border="0"

$$\displaystyle{b_{n} = \frac{2}{L}\int_{0}^{L}f(x)\sin{\frac{n{\pi}x}{L}}dx}$$ (3.4)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle{b_{n} = \frac{2}{\pi}\int_{0}^{\pi}sin(nx)dx = \frac{2}{\pi}[\int_{0}^{\pi/2}xsin(nx)dx + \int_{\pi/2}^{\pi}(\pi-x)sin(nx)dx]}$$ (3.5)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle{b_{n} = \frac{2}{\pi}[[\frac{sin(nx)-nxcos(nx)}{n^2}]_{0}^{\pi/2}+[\frac{-sin(nx)+n(\pi-x)cos(nx)}{n^2}]_{\pi/2}^{\pi}]}$$     (3.6)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Evaluating and simplifying the definite integral gives the following:


 * {| style="width:100%" border="0"

$$\displaystyle{b_{n} = \frac{2}{n^2\pi}[2sin(\frac{n\pi}{2})-sin(n\pi)]}$$ (3.7)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Evaluating equation (3.7) at $$ n = 1, 2, 3, 4, 5, 6, 7, 8\!$$


 * {| style="width:100%" border="0"

$$\displaystyle{b_{1} = \frac{4}{\pi}, b_{2} = 0, b_{3} = -\frac{4}{9\pi}, b_{4} = 0, b_{5} = \frac{4}{25\pi}, b_{6} = 0, b_{7} = -\frac{4}{49\pi}, b_{8} = 0}$$ (3.8)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Thus:


 * {| style="width:100%" border="1"

$$ f(x) = \frac{4}{\pi}(sinx - \frac{1}{9}sin3x + \frac{1}{25}sin5x - \frac{1}{49}sin7x +...) $$ (3.9)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Plots of truncated Fourier series: n = 2: n = 4: n = 8:

-17) By inspection, it is evident that the given function in problem 17 is even because for every $$ x\!$$ value that is negative, $$ f(x)\!$$ is positive. Therefore, $$ f(-x) = f(x)\!$$.

The following are the defining functions for the given graph:
 * {| style="width:100%" border="0"

$$\displaystyle{f(x)=\begin{cases}x+1 & \text{ if } -1<x<0\\ 1-x & \text{ if } 0<x<1 \end{cases}}$$ (3.10)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Because the function is even, we use equation (3.0) to develop the Fourier series, where the period is $$2\!$$ and $$ L = 1\!$$:


 * {| style="width:100%" border="0"

$$\displaystyle{f(x) = a_{0} + \sum_{n=1}^{\infty}a_{n}cos{\frac{n\pi}{L}x} = a_{0} + \sum_{n=1}^{\infty}a_{n}cos(n\pi x)}$$ (3.11)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Now we need to find $$a_{0}\!$$ using Euler's Formula:
 * {| style="width:100%" border="0"

$$\displaystyle{a_{0} = \frac{1}{2L}\int_{-L}^{L}f(x)dx}$$ (3.12)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle{a_{0} = \frac{1}{2}\int_{-1}^{1}f(x)dx = \frac{1}{2}[\int_{-1}^{0}(x+1)dx + \int_{0}^{1}(1-x)dx]}$$ (3.13)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle{a_{0} = \frac{1}{2}[[\frac{x^2}{2}+x]_{-1}^{0}+[x-\frac{x^2}{2}]_{0}^{1}]}$$     (3.14)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Evaluating the definite integral in equation (3.14) yields $$ a_{0}\!$$:


 * {| style="width:100%" border="0"

$$\displaystyle{a_{0} = \frac{1}{2}}$$ (3.15)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Now we find $$ a_{n}\!$$ using Euler's Formula:


 * {| style="width:100%" border="0"

$$\displaystyle{a_{n} = \frac{1}{L}\int_{-L}^{L}f(x)cos{\frac{n{\pi}x}{L}}dx = \int_{-1}^{1}f(x)cos(n\pi x)dx}$$ (3.16)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle{a_{n} = \int_{-1}^{0}(x+1)cos(n\pi x)dx + \int_{0}^{1}(1-x)cos(n\pi x)dx}$$ (3.17)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle{a_{n} = [\frac{\pi n(x+1)sin(\pi nx) + cos(\pi nx)}{n^2\pi^2}]_{-1}^{0} + [-\frac{\pi n(x-1)sin(\pi nx)+cos(\pi nx)}{n^2\pi^2}]_{0}^{1}}$$ (3.18)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Evaluating the definite integral in equation (3.18) yields the following equation for $$a_n\!$$:


 * {| style="width:100%" border="0"

$$\displaystyle{a_{n} = \frac{2(1-cos(\pi))}{n^2\pi^2}}$$ (3.19)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Evaluating equation (3.19) at $$ n = 1, 2, 3, 4, 5, 6, 7, 8\!$$ yields:


 * {| style="width:100%" border="0"

$$\displaystyle{a_{1} = \frac{4}{\pi^2}, a_{2} = 0, a_{3} = \frac{4}{9\pi^2}, a_{4} = 0, a_{5} = \frac{4}{25\pi^2}, a_{6} = 0, a_{7} = \frac{4}{49\pi^2}, a_{8} = 0}$$ (3.20)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Thus, the Fourier series is:


 * {| style="width:100%" border="1"

$$ f(x) = \frac{1}{2} + \frac{4}{\pi^2}(cos(\pi x) + \frac{1}{9}cos(3\pi x) + \frac{1}{25}cos(5\pi x) + \frac{1}{49}cos(7\pi x) + ...) $$ (3.21)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Plots of truncated Fourier series: n = 2: n = 4: n = 8:

Author
Solved and Typed By ---Egm4313.s12.team1.wyattling (talk) 06:19, 11 April 2012 (UTC) Reviewed By -Egm4313.s12.team1.armanious (talk) 15:18, 11 April 2012 (UTC)

Statement
Consider the L2-ODE-CC (5) p7b-7 with the window function f(x) p.9-8 as excitation:
 * {| style="width:100%" border="0"

$$\displaystyle{y''-3y'+2y=r(x)}$$ $$r(x)=f(x)$$ (4.0) and the initial conditions
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle{y(0)=1, y'(0)=0}$$ (4.1) 1. Find $$ y_n(x) $$ such that:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle{y''_n +ay'_n + by_n=r_n(x)}$$ (4.2) with the same initial conditions (4.1). Plot $$ y_n(x)$$ for n=2,4,8, for x in [0,10]. 2.Use the matlab command ode45 to integrate the L2-ODE-CC, and plot the numerical solution to compare with the analytical solution. Level 1:n=0,1
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Part 1 Solution
First, we shift the excitation f(x) to the left by introducing a new independent variable, t. This allows the period to start at zero.
 * {| style="width:100%" border="0"

$$\displaystyle t=x-\frac{1}{4}$$ (4.3) The piecewise representation of the window function is now (in terms of t) as follows:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle f(t)=A$$ for $$\displaystyle t=\left[ 0,2 \right]$$ $$\displaystyle f(t)=0$$ for $$\displaystyle t=\left[ 2,4 \right]$$ (4.4) To find the Fourier transform, the period of oscillation is determined:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle p=2L=4$$ (4.5) And the frequency of oscillation:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle \omega =\frac{2\pi }{p}=\frac{\pi }{2}$$ (4.6) The general form of a Fourier transform is
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle f(t)={{a}_{0}}+\sum\limits_{1}^{n}{\left[ {{a}_{n}}\cos (n\omega t)+{{b}_{n}}\sin (n\omega t) \right]}$$ (4.7) The following equations are given values of the constants in (4.7), evaluated for the function in (4.4):
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle {{a}_{0}}=\frac{1}{2L}\int\limits_{0}^{2L}{f(t)dt}=\frac{1}{4}\int\limits_{0}^{2}{Adt}=\frac{A}{4}(2-0)=\frac{A}{2}$$ (4.8)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{a}_{n}}=\frac{1}{L}\int\limits_{0}^{2L}{f(t)\cos (n\omega t)dt} =\frac{1}{2}\int\limits_{0}^{2}{A\cos (\frac{n\pi }{2}t)dt} =\frac{A}{n\pi }(\sin n\pi -\sin 0)$$ (4.9)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{b}_{n}}=\frac{1}{L}\int\limits_{0}^{2L}{f(t)\sin (n\omega t)dt}=\frac{1}{2}\int\limits_{0}^{2}{A\sin (\frac{n\pi }{2}t)dt}=-\frac{A}{n\pi }(\cos n\pi -\cos 0)=\frac{A}{n\pi }(1-\cos n\pi )$$ (4.10) In their simplest forms:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle {{a}_{n}}=0$$ (4.11)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{b}_{n}}=\left\{ 0,\frac{A}{\pi },0,\frac{A}{2\pi },0,\frac{A}{3\pi },0,\frac{A}{4\pi },...,0,\frac{2A}{n\pi } \right\}$$ (4.12) Plugging in r(x)=f(t) in (4.0),
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle y''-3y'+2y=\frac{1}{2}+\frac{2}{n\pi }\sin (\frac{n\pi }{2}t)$$ (4.13) From the general form, a particular solution to (4.13) and its derivatives are as follows:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{n}}={{A}_{n}}\cos (\frac{n\pi }{2}t)+{{B}_{n}}\sin (\frac{n\pi }{2}t)$$ (4.14)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{n}}'=-{{A}_{n}}\frac{n\pi }{2}\sin (\frac{n\pi }{2}t)+{{B}_{n}}\frac{n\pi }{2}\cos (\frac{n\pi }{2}t)$$ (4.15)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{n}}''=-{{A}_{n}}\frac{4}\cos (\frac{n\pi }{2}t)-{{B}_{n}}\frac{4}\sin (\frac{n\pi }{2}t)$$ (4.16)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

By plugging in (4.14), (4.15), and (4.16) into (4.13), we find An and Bn in terms of another constant, Cn:


 * {| style="width:100%" border="0"

$$\displaystyle {{A}_{n}}=6n\pi {{C}_{n}}$$ (4.17)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{B}_{n}}=(8-{{n}^{2}}{{\pi }^{2}}){{C}_{n}}$$ (4.18)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{C}_{n}}={{\left[ \tfrac{1}{8}n\pi ({{(8-{{n}^{2}}{{\pi }^{2}})}^{2}}+{{(6n\pi )}^{2}}) \right]}^{-1}}$$ (4.19) After substituting t with (4.3), the solutions are shown for n=2,4,8.
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{2}}={{y}_{1}}+{{y}_{2}}$$ (4.20)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{4}}={{y}_{1}}+{{y}_{2}}+{{y}_{3}}+{{y}_{4}}$$ (4.21)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{8}}={{y}_{1}}+{{y}_{2}}+{{y}_{3}}+{{y}_{4}}+{{y}_{5}}+{{y}_{6}}+{{y}_{7}}+{{y}_{8}}$$ (4.22)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Part 2 Solution
The excitation for n=0 and n=1 are the same, because $$a_1=b_1=0$$.

The homogeneous solution to (4.0) is
 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{h}}={{C}_{1}}{{\operatorname{e}}^{2t}}+{{C}_{2}}{{\operatorname{e}}^{t}}$$ (4.23) For $$r(x)=a_0$$, the particular solution is
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$\displaystyle {{y}_{p}}={{C}_{3}}$$ (4.24) Evaluating, we find that
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$${{C}_{1}}=-\tfrac{3}{4}$$ $${{C}_{2}}=\tfrac{3}{2}$$ $${{C}_{3}}=\tfrac{1}{4}$$ (4.25) Thus, the complete solution $$y_0$$ is
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$ {{y}_{0}}=\tfrac{1}{4}-\tfrac{3}{4}{{\operatorname{e}}^{2(x-\tfrac{1}{4})}}+\tfrac{3}{2}{{\operatorname{e}}^{x-\tfrac{1}{4}}}$$ (4.26) Integrating $$y_0$$ using MATLAB's ode45 command, the following plot is obtained for y:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Author
Solved and Typed By - Egm4313.s12.team1.essenwein (talk) 00:00, 11 April 2012 (UTC) Reviewed By - --Egm4313.s12.team1.wyattling (talk) 14:02, 11 April 2012 (UTC)

Statement
Continuation of R4.2, R4.3, R4.4. R4.2,p.7c-26:$$y_n(x)=y_{h,n}(x) + y_{p,n}(x)$$ For each value of n=3,5,9, redisplay the expressions for the 3 functions $$ y_{p,n}(x),y_{h,n}(x),y_n(x) $$, and plot these 3 functions separately over the interval $$ [0,20 \pi] $$. Exact solution: $$ y(x) = y_h(x) + y_p(x)$$ Redisplay the expressions for $$y_p(x),y_h(x), y(x)$$ Superimpose each of the above plots with that of the exact solution. R4.3,p.7c-28: Level 1: Understand and run the TA's code to produce a similar plot, but over a larger interval [0,10]. Do zoom-in plots about points $$ x=-0.5,0,+0.5 $$ and comment on the accuracy of different approximations. R4.4,p.7c-29: Level 1: Understand and run the TA's code to produce a similar plot, but over a larger interval [0.9,10], and for n=4,7. Do zoom-in plots about x=1,1.5,2,2.5 and comments on the accuracy of the approximations.

R4.2,p.7c-26
This first part of the problem is to redo 4.2:

In order to find the overall solution for the L2-ODE-CC corresponding to the Taylor series expansion of the sine function, both the homogenous and particular solutions must be found. The homogenous equation can be found through this method:
 * {| style="width:100%" border="0"

$$  \displaystyle y''-3y'+2y=0 $$     (5.0)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \lambda ^2-3\lambda+2=0 \rightarrow \lambda = 1,2 $$     (5.1)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_{h,n}(x)=C_1e^x+C_2e^{2x} $$     (5.2) Next, the excitation must be expanded to the desired n.  The following shows the excitation expanded to n=3, 5, and 9:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle r_3(x)=x-\frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{7!} $$     (5.3)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle r_5(x)=x-\frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{7!}+\frac{x^9}{9!}-\frac{x^{11}}{11!} $$     (5.4)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle r_9(x)=x-\frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{7!}+\frac{x^9}{9!}-\frac{x^{11}}{11!}+\frac{x^{13}}{13!}-\frac{x^{15}}{15!}+\frac{x^{17}}{17!}-\frac{x^{19}}{19!} $$     (5.5)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Next, the particular solution must be found. The particular solution will be of the form:
 * {| style="width:100%" border="0"

$$  \displaystyle y_{p,n}=\sum_{i=0}^{n}k_ix^i $$     (5.6) Using the derivation discussed in R4.1, a matrix equation in the form Ak = d can be found, where k is the matrix containing the coefficients of the particular solution and d is the matrix containing the coefficients of each power of x in the expansion of the excitation. The general formula for A is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \begin{bmatrix} b& a & 2 & &  & &\\ & b & 2a & 6 & & &\\ & & b & 3a & 12 & &\\ & &  & ... & ... & ... &\\ &  &  &  & b & a(n-1) & n(n-1)\\ & &  &  &  &b &an\\ & &  &  &  & &b \end{bmatrix} $$     (5.7) n=3 The matrix A is found as:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle A = \begin{bmatrix} 2 & -3 & 2 & &  &  &  & \\ & 2 & -6 & 6 &  &  &  & \\ &  & 2 & -9 & 12 &  & &  \\ &  &  & 2 &  -12& 20 &  & \\ &  &  &  & 2 & -15 & 30 & \\ &  &  &  &  & 2 & -18 & 42\\ &  &  &  &  &  & 2 &-21 \\ &  &  &  &  &  &  & 2 \end{bmatrix} $$     (5.8) Therefore, the matrix equation is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \begin{bmatrix} 2 & -3 & 2 & &  &  &  & \\ & 2 & -6 & 6 &  &  &  & \\ &  & 2 & -9 & 12 &  & &  \\ &  &  & 2 &  -12& 20 &  & \\ &  &  &  & 2 & -15 & 30 & \\ &  &  &  &  & 2 & -18 & 42\\ &  &  &  &  &  & 2 &-21 \\ &  &  &  &  &  &  & 2 \end{bmatrix} \begin{bmatrix} k_0\\ k_1\\ k_2\\ k_3\\ k_4\\ k_5\\ k_6\\ k_7 \end{bmatrix}= \begin{bmatrix} 0\\ 1\\ 0\\ -\frac{1}{6}\\ 0\\ \frac{1}{120}\\ 0\\ -\frac{1}{5040} \end{bmatrix} $$     (5.9)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle \begin{bmatrix} k_0\\ k_1\\ k_2\\ k_3\\ k_4\\ k_5\\ k_6\\ k_7 \end{bmatrix}= \begin{bmatrix} -0.1992\\ -0.3984\\ -0.3984\\ -0.0990\\ -0.0078\\ -0.0031\\ -0.0010\\ -9.9206\times 10^{-5} \end{bmatrix} $$     (5.10) Therefore the particular solution is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_{p,3}(x)=-9.9206\times 10^{-5}x^7-0.0010x^6-0.0031x^5-0.0078x^4-0.0990x^3-0.3984x^2-0.3984x-0.1992 $$     (5.11) And the overall solution for n=3 is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_{3}(x)=C_1e^x+C_2e^{2x}-9.9206\times 10^{-5}x^7-0.0010x^6-0.0031x^5-0.0078x^4-0.0990x^3-0.3984x^2-0.3984x-0.1992 $$     (5.12) Using the initial conditions:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle 1=C_1+C_2-0.1992$$ &  $$ 0=C_1+2C_2-0.3984 $$     (5.13) Solving yields $$C_1=2 $$  &  $$ C_2=-0.8008 \!$$. Therefore the overall solution is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle y_{3}(x)=2e^x-0.8008e^{2x}-9.9206\times 10^{-5}x^7-0.0010x^6-0.0031x^5-0.0078x^4-0.0990x^3-0.3984x^2-0.3984x-0.1992 $$     (5.14) n=5 The matrix equation is: Therefore the particular solution is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_{p,5}(x)=-1.2526\times 10^{-8}x^{11}-2.0668\times 10^{-7}x^{10}-1.0334\times 10^{-6}x^9-4.6503\times 10^{-6}x^8-1.1781\times 10^{-4}x^7-$$ $$0.0011x^6-0.0033x^5-0.0083x^4-0.0999x^3-0.3999x^2-0.3999x-0.2000 $$     (5.15) And the overall solution for n=5 is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_{5}(x)=C_1e^x+C_2e^{2x}-1.2526\times 10^{-8}x^{11}-2.0668\times 10^{-7}x^{10}-1.0334\times 10^{-6}x^9-4.6503\times 10^{-6}x^8$$ $$-1.1781\times 10^{-4}x^7-0.0011x^6-0.0033x^5-0.0083x^4-0.0999x^3-0.3999x^2-0.3999x-0.2000 $$     (5.16) Using the initial conditions:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Solving yields $$C_1=2.0001 \& C_2=-0.8001 \!$$. Therefore the overall solution is:
 * {| style="width:100%" border="1"

$$  \displaystyle y_{5}(x)=2.0001e^x-0.8001e^{2x}-1.2526\times 10^{-8}x^{11}-2.0668\times 10^{-7}x^{10}-1.0334\times 10^{-6}x^9-4.6503\times 10^{-6}x^8$$ $$-1.1781\times 10^{-4}x^7-0.0011x^6-0.0033x^5-0.0083x^4-0.0999x^3-0.3999x^2-0.3999x-0.2 $$     (5.18) n=9 The matrix equation is: Therefore the particular solution is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_{p,9}(x)=-4.1103\times 10^{-18}x^{19}-1.1714\times 10^{-16}x^{18}-1.0543\times 10^{-15}x^{17}-8.9615\times 10^{-15}x^{16}-4.5405\times 10^{-13}x^{15}-9.1408\times 10^{-12}x^{14}$$ $$-6.3985\times 10^{-11}x^{13}-4.1590\times 10^{-10}x^{12}-1.5021\times 10^{-8}x^{11}-2.2040\times 10^{-7}x^{10}-1.1020\times 10^{-6}x^9-4.9591\times 10^{-6}x^8$$ $$-1.1904\times 10^{-4}x^7-0.0011x^6-0.0033x^5-0.0083x^4-0.1000x^3-0.4000x^2-0.4000x-0.2000 $$     (5.19) And the overall solution for n=9 is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_{9}(x)=C_1e^x+C_2e^{2x}-4.1103\times 10^{-18}x^{19}-1.1714\times 10^{-16}x^{18}-1.0543\times 10^{-15}x^{17}-8.9615\times 10^{-15}x^{16}-4.5405\times 10^{-13}x^{15}-9.1408\times 10^{-12}x^{14}$$ $$-6.3985\times 10^{-11}x^{13}-4.1590\times 10^{-10}x^{12}-1.5021\times 10^{-8}x^{11}-2.2040\times 10^{-7}x^{10}-1.1020\times 10^{-6}x^9$$ $$-4.9591\times 10^{-6}x^8-1.1904\times 10^{-4}x^7-0.0011x^6-0.0033x^5-0.0083x^4-0.1000x^3-0.4000x^2-0.4000x-0.2000 $$     (5.20) Using the initial conditions:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle 1=C_1+C_2-0.2 \& 0=C_1+2C_2-0.4 $$     (5.21) Solving yields $$C_1=2 \& C_2=-0.8 \!$$. Therefore the overall solution is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="1"

$$  \displaystyle y_{9}(x)=2e^x-0.8e^{2x}-4.1103\times 10^{-18}x^{19}-1.1714\times 10^{-16}x^{18}-1.0543\times 10^{-15}x^{17}-$$ $$8.9615\times 10^{-15}x^{16}-4.5405\times 10^{-13}x^{15}-9.1408\times 10^{-12}x^{14}-6.3985\times 10^{-11}x^{13}$$ $$-4.1590\times 10^{-10}x^{12}-1.5021\times 10^{-8}x^{11}-2.2040\times 10^{-7}x^{10}-1.1020\times 10^{-6}x^9$$ $$-4.9591\times 10^{-6}x^8-1.1904\times 10^{-4}x^7-0.0011x^6-0.0033x^5-0.0083x^4-0.1000x^3-0.4000x^2-0.4000x-0.2000 $$     (5.22) The following plot shows the overall solution of the ODE at n=3,5,9 over the domain [0,4π]: Figure 4.2-2
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

The near perfect overlap between all three graphs shows that the approximations converge quickly for very low values of n.

Next, the exact ODE will be calculated to find the accuracy of the above approximations. The homogeneous solution is the same as the homogeneous solution above, only with different coefficients. The particular solution will be in the form:
 * {| style="width:100%" border="0"

$$  \displaystyle y_p(x)=K\cos x + M\sin x $$ (5.23)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_p'(x)=-K\sin x + M\cos x $$ (5.24)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_p''(x)=-K\cos x - M\sin x $$ (5.25) Plugging these values into the ODE:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle -K\cos x - M\sin x + 3K\sin x - 3M\cos x + 2K\cos x + 2M\sin x =\sin x $$ (5.26) Separating sine and cosine terms yield two linear equations which can be used to solve for the unknown coefficients:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle K-3M=0 $$ &  $$  3K+M=1 $$     (5.27) Solving these equations yield K = 0.3 and M = 0.1.
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y_p(x)=0.3\cos x + 0.1\sin x $$ (5.28) Therefore the exact overall solution is:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y(x)=C_1e^x + C_2e^{2x}+0.3\cos x + 0.1\sin x $$ (5.29)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle y'(x)=C_1e^x + 2C_2e^{2x}-0.3\sin x + 0.1\cos x $$ (5.30) Using the initial conditions yield the equations:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$  \displaystyle 1=C_1 + C_2 + 0.3 $$  &  $$  0=C_1 + 2C_2 + 0.1 $$     (5.31)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Solving these two equations yield $$C_1=1.5 \& C_2=-0.8 \!$$. Therefore the exact overall solution is:
 * {| style="width:100%" border="1"

$$  \displaystyle y(x)=1.5e^x - 0.8e^{2x}+0.3\cos x + 0.1\sin x $$ (5.32)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

The following figure shows the plot of the exact solution over the plot in Fig. 4.2-2: Figure 4.2-3

R4.3,p.7c-28
In this problem we are given the TA's Matlab code but must change the interval on which we are plotting. This give us the following code:

x=-1:0.01:1.5; y=log(1+x); x1=-1.5:0.01:1.5; y1=zeros(1,301); for i=1:4 for j=1:301 y1(j)=y1(j)-((-x1(j))^i)/i; end end y2=zeros(1,301); for i=1:7 for j=1:301 y2(j)=y2(j)-((-x1(j))^i)/i; end end y3=zeros(1,301); for i=1:11 for j=1:301 y3(j)=y3(j)-((-x1(j))^i)/i; end end y4=zeros(1,301); for i=1:16 for j=1:301 y4(j)=y4(j)-((-x1(j))^i)/i; end end h=plot(x,y); orange = [1 0.5 0.2]; set(h,'Color',orange); hold on; plot(x1,y1,'r'); plot(x1,y2,'g'); plot(x1,y3,'b'); plot(x1,y4,'c'); axis([0 10 -10 10]); legend('log(1+x)','T_4','T_7','T_1_1','T_1_6'); grid on; Zooming at -0.5 This graph shows fairly accurate results with slight deviation, with largest deviation occurring at the T4 term. Zooming about point 0 In this graph, only T16 shows, thus meaning that it is the best approximation. Zooming at 0.5 All lines appear to be fairly good approximations with the exception of T4 which deviates from the rest.

R4.4,p.7c-29
Again we must adjust the TA's Matlab code:

syms x f=log(x+1); fT1=taylor(f,5,1); fT2=taylor(f,8,1); fT3=taylor(f,12,1); X=.9:0.1:10; Y(:,1)=subs(fT1,'x',X); Y(:,2)=subs(fT2,'x',X); Y(:,3)=subs(fT3,'x',X); Y(:,4)=log(1+X); figure plot(X,Y); axis([.9,10,-10,10]);

Zooming at 1 Since the three lines are all stacked, both n=4 and n=7 are good approximations at x=1. Zooming at 1.5 Again, since the three lines are all stacked, both n=4 and n=7 are good approximations at x=1.5. Zooming at 2 Here, both can again be said to be good approximations, however, you can see some very slight deviation in n=4. Zooming at 2.5 In this plot, n=7 is still a pretty strong approximation, while n=4 is beginning to deviate at a fairly rapid rate and thus will soon be a poor approximation.

Author
Solved and Typed By - --Egm4313.s12.team1.rosenberg (talk) 01:19, 11 April 2012 (UTC) Reviewed By ---Egm4313.s12.team1.durrance (talk) 19:12, 11 April 2012 (UTC)

Statement
Verify (2) p.10-7 as follows: Given the ODE:
 * {| style="width:100%" border="0"

$$  \displaystyle y''_{p}+4y'_{p}+13y_{p}=2e^{-2x}cos(3x) $$     (6.1) 1)Simplify the 1st term $$y''_p$$ on the lhs of (2) p.10-4. 2)Simplify the 2nd term $$4y'_p $$, and combine with the simplified 1st term. 3)Finally, add the 3rd term $$13y_p$$ 4)Find the final expression for $$y_p(x)$$
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Solution
From [https://elearning2.courses.ufl.edu/access/content/group/UFL-EGM4313-5641-12012/Lecture%20Notes/iea.s12.sec10-1.djvu Lecture 10 Pg. 9]

To find the particular solution of an ODE we can use the method of undetermined coefficients. In report 3 we used the table 2.1 below to find the standard form of the particular solution given the form of the excitation $$r(x)\!$$. Table 2.1 So the general form of the particular solution is:
 * {| style="width:100%" border="0"

$$  \displaystyle y_{p}(x)=xe^{-2x}[Mcos3x+Nsin3x] $$     (6.2)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Task 1 asks to simplify they first term of the given ODE which is $$y_{p}''(x)\!$$. This is done by taking the second derivative of $$y_{p}(x)\!$$ and plugging it into equation 6.1
 * {| style="width:100%" border="0"

$$  \displaystyle y''_{p}(x)=\frac{\partial }{\partial x}\frac{\partial (xe^{-2x}[Mcos3x+Nsin3x])}{\partial x} $$ (6.3)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

$$y_{p}''(x)\!$$ simplifies to:
 * {| style="width:100%" border="0"

$$  \displaystyle y''_{p}(x)=e^{-2x}(sin(3x)[6M(2x-1)-N(5x+4)]-cos(3x)[M(5x+4)+6N(2x-1)]) $$     (6.4)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Now we must do the same for the second term $$4y_{p}'(x)\!$$.
 * {| style="width:100%" border="0"

$$  \displaystyle 4y'_{p}(x)=\frac{\partial (4xe^{-2x}[Mcos3x+Nsin3x])}{\partial x} $$ (6.5)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$  \displaystyle 4y'_{p}(x)=4e^{-2x}[sin(3x)(-3Mx-2Nx+N)+cos(3x)(-2Mx+M+3Nx)] $$     (6.6)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Add the second term to the first to continue building the original ODE (equation 6.1).
 * {| style="width:100%" border="0"

$$  \displaystyle y''_{p}(x)+4y'_{p}(x)=\frac{\partial }{\partial x}\frac{\partial (xe^{-2x}[Mcos3x+Nsin3x])}{\partial x}+\frac{\partial (4xe^{-2x}[Mcos3x+Nsin3x])}{\partial x} $$ (6.7)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$  \displaystyle y''_{p}(x)+4y'_{p}(x)=4e^{-2x}[sin(3x)(-3Mx-2Nx+N)+cos(3x)(-2Mx+M+3Nx)]$$ $$+e^{-2x}(sin(3x)[6M(2x-1)-N(5x+4)]-cos(3x)[M(5x+4)+6N(2x-1)]) $$     (6.8)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Finally add the third term $$13y_{p}(x)\!$$ to the other two to have the lhs of the ODE 6.1
 * {| style="width:100%" border="0"

$$  \displaystyle 13y_{p}(x)=13xe^{-2x}[Mcos3x+Nsin3x] $$     (6.9)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$  \displaystyle y''_{p}(x)+4y'_{p}(x)+13y_{p}(x)=e^{-2x}(sin(3x)[6M(2x-1)-N(5x+4)]-cos(3x)[M(5x+4)+6N(2x-1)])+4e^{-2x}[sin(3x)(-3Mx-2Nx+N)$$ $$+cos(3x)(-2Mx+M+3Nx)]+13xe^{-2x}[Mcos3x+Nsin3x] $$     (6.10)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$  \displaystyle y_{p}(x)=xe^{-2x}[Mcos3x+Nsin3x]=6e^{-2x}[Ncos(3x)-Msin(3x)] $$     (6.11)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

This simplifies to the particular solution below:
 * {| style="width:100%" border="1"

$$  \displaystyle y_{p}(x)=6e^{-2x}[Ncos(3x)-Msin(3x)] $$     (6.12)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Author
Solved and Typed By ---Egm4313.s12.team1.stewart (talk) 23:50, 10 April 2012 (UTC) Reviewed By - Egm4313.s12.team1.essenwein (talk) 18:12, 11 April 2012 (UTC)

Statement
Find the separated ODEs for the heat equation:
 * {| style="width:100%" border="0"

$$ \frac{\partial u}{\partial t}=\kappa \frac{\partial^2 u}{\partial x^2} $$     (7.0) Where $$\kappa$$ is the coefficient of heat conductivity.
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Solution
Assume:
 * {| style="width:100%" border="0"

$$ u(x,t)=F(x) \cdot G(t) $$     (7.1) So, by take partial derivatives of Eq. (7.1) with respect to the individual variables:
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }
 * {| style="width:100%" border="0"

$$ \frac{\partial u(x,t)}{\partial t}=F(x) \cdot \dot G(t) $$     (7.2)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }


 * {| style="width:100%" border="0"

$$ \frac{\partial^2 u(x,t)}{\partial x^2}=F''(x) \cdot G(t) $$     (7.3)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

By substituting the partial derivatives back into Eq. (7.0):
 * {| style="width:100%" border="0"

$$ F(x) \cdot \dot G(t)=\kappa F''(x) \cdot G(t) $$     (7.4)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Separating the equations to similar variables yields:
 * {| style="width:100%" border="0"

$$ \frac{\dot G(t)}{G(t)}=\kappa \frac{F''(x)}{F(x)}=C \text{ (constant)} $$     (7.4a)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Therefore the separated ODEs are:
 * {| style="width:100%" border="1"

$$ \dot G(t)-CG(t)=0$$ $$ F''(x)-\frac{C}{\kappa}F(x)=0 $$     (7.5)
 * style="width:95%" |
 * style="width:95%" |
 * <p style="text-align:right">
 * }

Author
Solved and Typed By - --Egm4313.s12.team1.durrance (talk) 02:13, 9 April 2012 (UTC) Reviewed By - --Egm4313.s12.team1.stewart (talk) 23:48, 10 April 2012 (UTC)