[SOLVED] 5DV005-Lab 6

Description

Problem 1 Copy rintmwe1 into /lab6/work/l6p1.m and adapt it to the problem of computing the integral

numerically using the trapezoidal rule.

1. What evidence do you find to support support the conjecture that there exist an asymptotic error expansion of the form

T − A_h = αh^p + βh^q + O(h^r),              0 < p < q < r.

2. Based on the numerical evidence, what is a reasonable value of p?
3. Based on the numerical evidence, what is a reasonable value of q?
4. What is the smallest value of k for which the integral can be computed with a relative error less than τ = 10⁻⁶?
5. Compute the exact value of the integral and include this information in l6p2.
6. Is the behavior of Richardson’s fraction related to the quality of Richardson’s error estimate?

Problem 2         rintmwe1.m into /work/l6p2.m and adapt it to the problem of com-

puting the integral

using the trapezoidal rule as your approximation A_h.

1. What evidence do you find to support support the conjecture that there exist an asymptotic error expansion of the form

T − A_h = αh^p + βh^q + O(h^r),              0 < p < q < r.

2. Based on the numerical evidence, what is a reasonable value of p?
3. Based on the numerical evidence, what is a reasonable value of q?
4. What is the smallest value of k for which the integral can be computed with a relative error less than τ = 10⁻⁶?
5. Compute the exact value of the integral and include this information in l6p1. Hint: It is quite easy to compute the integral if you make a drawing of the graph first.
6. Is the behavior of Richardson’s fraction related to the quality of Richardson’s error estimate?

Problem 3         rintmwe1.m into /work/l6p3.m and adapt it to the problem of com-

puting the integral

using the trapezoidal rule as your approximation A_h.

1. What evidence do you find to support support the conjecture that there exist an asymptotic error expansion of the form

T − A_h = αh^p + βh^q + O(h^r),              0 < p < q < r.

2. Based on the numerical evidence, what is a reasonable value of p? 3. Based on the numerical evidence, what is a reasonable value of q?
3. Why is Richardson’s fraction not close to 2^p for small values of k?
4. Why is Richardson’s fraction not close to 2^p for very large values of k?.
5. What is the smallest value of k for which the integral can be computed with a relative error less than τ = 10⁻⁶?

Problem 4 rdifmwe1 into /work/l6p4.m and adapt it to the problem of computing the the target T = f⁰(x), where f is you favorite differentiable function and x is your favorite real number using the mysterious rule

where A_h is your favorite rule for computing f⁰(x) which obeys an asymptotic error expansion of the form

T − A_h = αh^p + βh^q + O(h^r),              0 < p < q < r.

1. What evidence can you uncover that suggests that M_h obeys an asymptotic error expansion of the form

T − M_h = ¯αh^q + βh^{¯ r} + O(h^s),            0 < q < r < s.                  (1)

2. Based on your numerical evidence, what is a reasonable value of q? 3. Based on your numerical evidence, what is a reasonable value of r?
3. Include the exact value of the derivative of f in the script.
4. Examine the relationship between Richardson’s fraction and the quality of the error estimate.