Log In Start studying!

Select your language

Suggested languages for you:
Deutsch (DE)
Deutsch (UK)
Deutsch (US)

Americas

Europe

Answers without the blur. Sign up and see all textbooks for free! → Illustration

Q18E

Expert-verified
Fundamentals Of Differential Equations And Boundary Value Problems
Found in: Page 131
Fundamentals Of Differential Equations And Boundary Value Problems

Fundamentals Of Differential Equations And Boundary Value Problems

Book edition 9th
Author(s) R. Kent Nagle, Edward B. Saff, Arthur David Snider
Pages 616 pages
ISBN 9780321977069

Answers without the blur.

Just sign up for free and you're in.

Register for Free I’ll do it later
Illustration

Short Answer

Local versus Global Error. In deriving formula (4) for Euler’s method, a rectangle was used to approximate the area under a curve (see Figure 3.14). With

\({\bf{g(t) = f(t,f(t))}}\) , this approximation can be written as \(\int\limits_{{{\bf{x}}_{\bf{n}}}}^{{{\bf{x}}_{{\bf{n + 1}}}}} {{\bf{g(t)dt}} \approx {\bf{hg(}}{{\bf{x}}_{\bf{n}}}{\bf{)}}} \)where \({\bf{h = }}{{\bf{x}}_{{\bf{n + 1}}}}{\bf{ - }}{{\bf{x}}_{\bf{n}}}\) .

  1. Show that if g has a continuous derivative that is bounded in absolute value by B, then the rectangle approximation has error\(\left( {\bf{O}} \right){{\bf{h}}^{\bf{2}}}\); that is, for some constant M, \(\left| {\int\limits_{{{\bf{x}}_{\bf{n}}}}^{{{\bf{x}}_{{\bf{n + 1}}}}} {{\bf{g(t)dt - hg(}}{{\bf{x}}_{\bf{n}}}{\bf{)}}} } \right| \le {\bf{M}}{{\bf{h}}^{\bf{2}}}\).This is called the local truncation error of the scheme. [Hint: Write \(\int\limits_{{{\bf{x}}_{\bf{n}}}}^{{{\bf{x}}_{{\bf{n + 1}}}}} {{\bf{g(t)dt - hg(}}{{\bf{x}}_{\bf{n}}}{\bf{)}}} {\bf{ = }}\int\limits_{{{\bf{x}}_{\bf{n}}}}^{{{\bf{x}}_{{\bf{n + 1}}}}} {\left[ {{\bf{g(t)dt - g(}}{{\bf{x}}_{\bf{n}}}{\bf{)}}} \right]{\bf{dt}}} \). Next, using the mean value theorem, show that\(\left| {{\bf{g(t)dt - g(}}{{\bf{x}}_{\bf{n}}}{\bf{)}}} \right| \le {\bf{B}}\left| {{\bf{t - }}{{\bf{x}}_{\bf{n}}}} \right|\) . Then integrate to obtain the error bound\(\left( {\frac{{\bf{B}}}{{\bf{2}}}} \right){{\bf{h}}^{\bf{2}}}\).]
  2. In applying Euler’s method, local truncation errors occur in each step of the process and are propagated throughout the further computations. Show that the sum of the local truncation errors in part (a) that arise after n steps is (O)h. This is the global error, which is the same as the [ss1] [m2] convergence rate of Euler’s method.

  1. The required result gets by the mean value theorem i.e., \(\left| {\int\limits_{{x_n}}^{{x_{n + 1}}} {g\,\left( t \right)\,dt - hg\,\left( {{x_n}} \right)} } \right| \le M\,{h^2}\).
  2. The sum of the local truncation error arises after n steps i.e., \(O\left( h \right)\).
See the step by step solution

Step by Step Solution

TABLE OF CONTENTS :
TABLE OF CONTENTS

Step 1: Important hint.

For the solution apply mean value theorem.

Step 2: Find the integral value.

Here \(\int\limits_{{x_n}}^{{x_{n + 1}}} {g\left( t \right)dt \approx hg\left( {{x_n}} \right)} \) where \(h = {x_{n + 1}} - {x_n}\)

Assume that g has a continuous derivative that is bounded by B.

\(\left| {g'\left( y \right)} \right| \le B\).

According to the mean value theorem

\(\begin{array}{c}\frac{{g\left( t \right) - g\left( {{x_n}} \right)}}{{t - {x_n}}} = g'\left( y \right)\\g\left( t \right) - g\left( {{x_n}} \right) = g'\left( y \right).\left( {t - {x_n}} \right)\\\left| {g\left( t \right) - g\left( {{x_n}} \right)} \right| \le B\left| {t - {x_n}} \right|\end{array}\)

Step 3: Apply mode on an integral part.

Now,

\(\begin{array}{c}\left| {\int\limits_{{x_n}}^{{x_{n + 1}}} {g\left( t \right)dt - hg\left( {{x_n}} \right)} } \right| = \left| {\int\limits_{{x_n}}^{{x_{n + 1}}} {\left[ {g\left( t \right)dt - g\left( {{x_n}} \right)} \right]dt} } \right|\\ \le \int\limits_{{x_n}}^{{x_{n + 1}}} {\left| {g\left( t \right)dt - g\left( {{x_n}} \right)} \right|dt} \\ \le \int\limits_{{x_n}}^{{x_{n + 1}}} {B\left| {t - {x_n}} \right|dt} \\ = \left[ {\frac{{B\left( {t - {x_n}} \right)\left| {t - {x_n}} \right|}}{2}} \right]_{{x_n}}^{{x_{n + 1}}}\end{array}\)

\(\begin{array}{c} = - \left[ {\frac{{B\left( {{x_n} - {x_{n + 1}}} \right)\left| {{x_n} - {x_{n + 1}}} \right|}}{2}} \right]\\ = \frac{{B{h^2}}}{2}\end{array}\)

Let \(M = \frac{B}{2}\) then

\(\left| {\int\limits_{{x_n}}^{{x_{n + 1}}} {g\left( t \right)dt - hg\left( {{x_n}} \right)} } \right| \le M{h^2}\)

Step 4: use the concept of part (a).

Let \(a = {x_n},b = {x_{n + 1}}\)

From part (a) the local truncation error of the scheme is \(\frac{{\left( B \right){h^2}}}{2}\) and \(h = \frac{{b - a}}{n}\).

Step 5: Finding the result apply the result on n steps.

After n steps the result is

\(\begin{array}{c} = h\frac{{B\left( {b - a} \right)}}{2}\\ = O\left( h \right)\end{array}\)

Hence, the sum of the local truncation error arises after n steps i.e., \(O\left( h \right)\).

Most popular questions for Math Textbooks

An RL circuit with a 5 resistor and a 0.05-H inductor carries a current of 1 A at t = 0, at which time a voltage source E(t) = 5 cos 120t V is added. Determine the subsequent inductor current and voltage.

In Problems 23–27, assume that the rate of decay of a radioactive substance is proportional to the amount of the substance present. The half-life of a radioactive substance is the time it takes for one-half of the substance to disintegrate. If initially there are 50 g of a radioactive substance and after 3 days there are only 10 g remaining, what percentage of the original amount remains after 4 days?

In Problem 21 it is observed that when the velocity of the sailboat reaches 5 m/sec, the boat begins to rise out of the water and “plane.” When this happens, the proportionality constant for the water resistance drops to b0 = 60 N-sec/m. Now find the equation of motion of the sailboat. What is the limiting velocity of the sailboat under this wind as it is planning?

In Problem 9, suppose we have the additional information that the population of splake in 2004 was estimated to be 5000. Use a logistic model to estimate the population of splake in the year 2020. What is the predicted limiting population? [Hint: Use the formulas in Problem 12.

An object at rest on an inclined plane will not slide until the component of the gravitational force down the incline is sufficient to overcome the force due to static friction. Static friction is governed by an experimental law somewhat like that of kinetic friction (Problem 18); it has a magnitude of at most N, where m is the coefficient of static friction and N is, again, the magnitude of the normal force exerted by the surface on the object. If the plane is inclined at an angle a, determine the critical value 0 for which the object will slide if a>ao but will not move for a<ao.
Icon

Want to see more solutions like these?

Sign up for free to discover our expert answers
Get Started - It’s free

Recommended explanations on Math Textbooks

Decision Maths

Read explanation

Calculus

Read explanation

Probability and Statistics

Read explanation

Statistics

Read explanation

Mechanics Maths

Read explanation

Geometry

Read explanation

Pure Maths

Read explanation

94% of StudySmarter users get better grades.

Sign up for free
94% of StudySmarter users get better grades.