Dr. Kevin G. TeBeest
As always, you are expected to rework my examples successfully on your own BEFORE attempting the homework.
- Answer each in a complete and coherent sentence. What are the global errors of:
(You should know these.)
- the trapezoidal rule,
- the rectangle method,
- Simpson's3/8 rule, and
- Simpson's1/3 rule?
- When using Simpson's3/8 rule, what must be true of the number of subintervals? Answer in a complete and coherent sentence.
- Do Problems 3–6 in this problem set.
- Modify the Maple code for Simpson's1/3 rule (see this previous assignment) to turn it into Simpson's3/8 rule.
Then use it to approximate the sample problem (integral of sin x on the interval [0, π]) using 6 subintervals.
You should get the answer we obtained in class. (Also see the pseudocode for Simpson's3/8 rule.)
- Now use your Simpson's3/8 code to re-do Problem 4 above using 12 subintervals. You should get the answer we obtained in class.
- Use the Simpson's3/8 rule code to approximate the integral of e–x2 on the interval [0,3] using:
- 9 subintervals Answer: 0.88619 38607 33682 32028
- 18 subintervals Answer: 0.88620 71620 01857 85809
- Use Richardson extrapolation to obtain an improved estimate.
Answer: +0.88620 80487 53069 56062- Without knowing the actual errors, determine how much more accurate we expect result (b) to be than result (a). Answer in a complete and coherent sentence.
- What is the error estimate of the result in (b)? Answer: +0.000000 886751 211702 520374
- What is the error estimate of the result in (b) in parts per million? Answer: +1.0006 ppm
- Play around with the code to approximate integrals of other functions on different intervals.
- In class I derived the formula for the trapezoidal rule. In a previous assignment, you were supposed to derive the formula for Simpson's1/3 rule. Now you will derive the composite formula for Simpson's3/8 rule:
- First construct the Newton-Gregory interpolating polynomial P[03](x) containing points indexed [0,1,2,3]. It approximates function f(x) on Section 1.
- Integrate the polynomial on Section 1: [x0, x3] and simplify the result.
You should obtain I1 = (3/8) · h · ( f0 + 3 f1 + 3 f2 + f3 ).
(You might want to review Problem 3 in this Assignment.)
- Use the formula for I1 on Section 1 as the template to construct the formulas for I2 on Section 2, I3 on Section 3, etc.
- Add I1 + I2 + I3 + · · · to obtain the composite Simpson's3/8 rule.
- Answer each in a complete and coherent sentence. How much do we expect the accuracy to increase if we double the number of subintervals using:
- Simpson's3/8 rule?
- the trapezoidal rule?
- Simpson's1/3 rule?
- the rectangle method learned in Calc-2?
Did you omit a VERY IMPORTANT WORD in your answers above?
- Answer each in a complete and coherent sentence. How much do we expect the accuracy to increase if we quadruple the number of subintervals using:
- Simpson's3/8 rule?
- the trapezoidal rule?
- Simpson's1/3 rule?
- the rectangle method learned in Calc-2?
Did you forget a VERY IMPORTANT WORD in your answers above?
- Is the rectangle method for approximating integrals a NewtonCotes method? Explain why or why not.
Newton-Cotes Methods: See a summary of some of the NewtonCotes integration rules.