University of Florida/Egm6341/s10.team2.niki/HW3

Problem Statement

Pg. 17-2

See also HW1-problem 8[[1]] Using the error estimates of

• Taylor Series
• Composite Trapezoidal rule
• Composite Simpson's rule

estimate $𝒏$ such that ${𝑬}_{𝒏}\mathit{=}𝑰\mathit{-}{𝑰}_{𝒏}\mathit{=}𝑶\mathit{(}10^{\mathit{-}6}\mathit{)}$

Problem Solution

Solution: Taylor series

Error is defined as

$E_n=I-I_n=\underset{a}{\overset{b}{\int }}[f(x)-f_n(x)]dx$

For the Taylor series, from the discussion on p6-4, we know that the error is nothing but the remainder of the Taylor series integrated over the given interval i.e

$\underset{a}{\overset{b}{\int }}{R}_{n+1}(x)=\underset{a}{\overset{b}{\int }}[\frac{(x-x_0{)}^{n+1}}{(n+1)!}{f}^{(n+1)}(\xi )]$ where $\xi \epsilon [a,b]$

(1 p2-3)

For the given function the error is

$E_n=\underset{0}{\overset{1}{\int }}\underset{w(x)}{\underbrace{\frac{x^n}{(n+1)!}}}\underset{g(x)}{\underbrace{e^{\xi (x)}}}dx$

(2)

Using the Integral Mean Value Theorem, we have

$E_n=e^{\xi (x=\alpha )}\underset{0}{\overset{1}{\int }}\frac{x^n}{(n+1)!}dx$

(3)

integrating we get,

$E_n=e^{\xi (x=\alpha )}\frac{1}{(n+1)!(n+1)}$

(4)

This function has a minimum when $\alpha =0$ and maximum when $\alpha =1$ thus we have

$\frac{1}{(n+1)!(n+1)}\le E_n\le \frac{e}{(n+1)!(n+1)}$

(5)

Setting the upper bound of the error to $10^{-6}$ we have

$E_n\le \frac{e}{(n+1)!(n+1)}=10^{-6}$

(6)

Solving this equation for n we get

$(n+1)!(n+1)=e*10^6\Rightarrow n\approx 7.92=7$

(A)

Solution:Composite Trapezoidal rule

From the discussion on page 16-3 we have

$\left|E_n^1\right|\le \frac{(b-a)}{12n^2}{M}_2$

(1)

where $M_2=max|f^{(2)}(\zeta )|$ for $\zeta \epsilon [a,b]$

For the given function $f(x)=\frac{e^x-1}{x}$, we have

$f^{(2)}(x)=\frac{e^x[x^2-2x+2]-2}{x^3}$

(2)

Evaluating over the given interval it is seen that the function $f^{(2)}(x)$ has maximum value at $x=1$

$M_2=f^{(2)}(x=1)=\frac{e[1+2-2]-2}{1}=0.718282$

(3)

Setting the error to the $10^{-6}$ and solving for $n$ we get

$\frac{(1-0{)}^3}{12n^2}(0.718282)=10^{-6}\Rightarrow n\approx 244.657=245$

(B)

Solution :Composite Simpson's Rule

We have from p 17-2, the error estimate of the Composite Simpson's rule as

$\left|E_n^2\right|\le \frac{(b-a{)}^5}{2880n^4}{M}_4$

(1)

where $M_4=max|f^{(4)}(\zeta )|$ for $\zeta \epsilon [a,b]$

For the given function $f(x)=\frac{e^x-1}{x}$, we have

$f^{(4)}(x)=\frac{e^x[x^4-4x^3+12x^2-24x]-24}{x^5}$

(2)

Evaluating over the given interval it is seen that the function $f^{(2)}(x)$ has maximum value at $x=1$

$M_4=f^{(4)}(x=1)=\frac{e[1-4+12-24+24]-24}{1}=9(e)-24=0.464536$

(3)

Setting the error to the $10^{-6}$ and solving for $n$ we get

$\frac{(1-0{)}^5}{2880n^4}(0.464536)=10^{-6}\Rightarrow n\approx 3.35735=3$

(c)

Part 2:Numerical determination of power of h

In order to verify the power of $𝒉$ in the error, data from Problem 8 of HW1 is used.

In the case of a Semilog plot (log(y) vs x), an equation of the form

$𝒍𝒐𝒈\mathit{(}𝒀\mathit{)}\mathit{=}𝒂𝑿\mathit{+}𝒃$

such that

$𝒀\mathit{=}{𝒆}^{\mathit{(}𝒂𝑿\mathit{+}𝒃\mathit{)}}\mathit{=}{𝒆}^{𝒃}\mathit{(}{𝒆}^{𝒂𝑿}\mathit{)}\mathit{=}𝒌{𝒆}^{𝒂𝑿}$

From the above equation it is seen that if the plot on a semilog graph is a straight line then the relationship between the two variables is exponential.

A log-log plot(log(y) vs log(x)) is a similar plot, for which the equation is of the form

$𝒍𝒐𝒈\mathit{(}𝒀\mathit{)}\mathit{=}𝒂\mathit{[}𝒍𝒐𝒈\mathit{(}𝑿\mathit{)}\mathit{]}\mathit{+}𝒃$

such that, $𝒀\mathit{=}𝒌{𝒙}^{𝒂};𝒌\mathit{=}{𝒆}^{𝒃}$

This discussion is used in the interpretation of the graphs given below.

It is readily seen that the slope of the line in the log-log graph is the power of the x variable.

Composite trapezoidal Rule

Given below is the data from the numerical evaluation of the given function using Composite Trapezoidal Rule.

Template:Center topComposite Trapezoidal RuleTemplate:Center bottom
No. of terms $\mathit{(}𝒏\mathit{)}$	Absolute Error $\mathit{(}\left|E_n^1\right|\mathit{)}$	$𝒉$
2	1.3282917278	0.5
4	1.3205046195	0.25
8	1.3185530869	0.125
16	1.3180649052	0.0625
32	1.3179428411	0.03125
64	1.3179123240	0.015625
128	1.3179046946	0.0078125
256	1.3179027872	0.00390625
512	1.3179023104	0.001953125
1024	1.3179021912	0.000976563

First we plot a semilog graph for the data. The graph is shown below:

File:Niki-trap-semilog.png

The graph is not a straight line which implies that the relationship between $𝒉$ and $𝒍𝒐𝒈\mathit{(}𝒆𝒓𝒓𝒐𝒓\mathit{)}$ is not exponential. Hence we plot a log-log graph as below:

File:Niki-trap-log.png

This is seen to be linear. A straight line if fitted to the data the equation of which is given above. From the discussion above, we see that the slope of line is 2.1447 which is very close to the analytical value of 2.

Composite Simpson's Rule

Given below is the data obtained from HW1 for the COmposite Simpsons Rule.

Template:Center topComposite Simpson's RuleTemplate:Center bottom
No. of terms $\mathit{(}𝒏\mathit{)}$	Absolute Error $\mathit{(}\left|E_n^2\right|\mathit{)}$	$𝒉$
2	1.318008666	0.5
4	1.317908917	0.25
8	1.317902576	0.125
16	1.317902178	0.0625

Plotting a semilog graph of $𝒍𝒐𝒈\mathit{(}𝒆𝒓𝒓𝒐𝒓\mathit{)}$ against $𝒉$ we see that it is non-linear as in the case of the Composite Trapezoidal Rule.

File:Niki-CSR-semilog.png

Thus, plotting the log-log graph as below,

File:Niki-CSR-log.png

we see that the slope of the line is 4.0881 which is very close to the analytically determined value of 4.

Taylor Series

Given below is the data for the Taylor series method. Eventhough $𝒉$ is not used in the Taylor series method,defining $𝒉$ based on the number of terms of the series, it is seen that the error and h are not related exponentially or by a power relation i.e. both semilog and log log plots are not linear.

Template:Center topTaylor SeriesTemplate:Center bottom
No. of terms $\mathit{(}𝒏\mathit{)}$	Absolute Error $\mathit{(}\left|E_n^2\right|\mathit{)}$	$𝒉$
2	1.2500000	0.5
4	1.3159722222	0.25
8	1.3179018152	0.125
16	1.3179021515	0.0625
32	1.3179021515	0.03125

--Egm6341.s10.team2.niki 08:10, 17 February 2010 (UTC)