5DV005- Lab 6 Solved

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 − Ah = αhp + βhq + O(hr),              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−6?

You must explain why your error estimate is reliable!

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 Ah.

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

                                                        T − Ah = αhp + βhq + O(hr),              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−6?

You must explain why your error estimate is reliable

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 Ah.

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

                                                        T − Ah = αhp + βhq + O(hr),              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.    Why is Richardson’s fraction not close to 2p for small values of k?

5.    Why is Richardson’s fraction not close to 2p for very large values of k?.

6.    What is the smallest value of k for which the integral can be computed with a relative error less than τ = 10−6?

You must explain why your error estimate is reliable!


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

 

where Ah is your favorite rule for computing f0(x) which obeys an asymptotic error expansion of the form

                                                  T − Ah = αhp + βhq + O(hr),              0 < p < q < r.

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

                                                      T − Mh = ¯αhq + βh¯ r + O(hs),            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?

4.    Include the exact value of the derivative of f in the script.

5.    Examine the relationship between Richardson’s fraction and the quality of the error estimate.