Homework
7 (due 3/11)
CSC 431
This week we finished chapter 6 (6.5.3), covered parts of chapter 7
(7.1, 7.3.1, 7.3.2, parts of 7.3.4/7.3.5) and chapter 8 (8.1, 8.3). Next week we
will cover some of the remaining topics in chapter 8 and go on to chapter 9.
Also, next week I will hand out the take-home final (due one week later).
During class we saw some of Heath's simulations:
and we discussed some of my solutions to homework 6:
Bernoulli and the accuracy
problem.
Submission: you can submit the homework by hardcopy
in class or by sending it to me as an email.
Note: for all of the problems on this homework you can use
the computer, but still explain your steps.
1. [Horner's Scheme, 5pt] Rewrite f(x) 5x^4 - 10x^3 + 5x^2 - x + 7 using
Horner's scheme and evaluate f(3) using Horner's scheme.
2. [Interpolation, 20pt] You have observed the following data:
t |
y = f(t) |
-1 |
10 |
0 |
4 |
1 |
6 |
2 |
2 |
- Using polynomial interpolation (Vandermonde matrix), find an
interpolant in the monomoial basis.
- Using Lagrange interpolation, find an interpolant in the Lagrange
basis.
- How many additions/multiplications are needed to evaluate each of
the two interpolants (assuming that you pre-compute divisors in the
Lagrange interpolant and evaluate the polynomial interpolant using
Horner's scheme)?
3. [Integration, 25pt] Newton derived the following formula to compute
pi:
-
- Use Midpoint, trapezoid and Simpson's quadrature rules to
approximate the integral and find an approximation to pi. Compute errors
and compare.
- Using Gaussian quadrature with two and three points to approximate
the integral and obtain an approximation to pi. We derived the locations
of two points for the interval [-1,1]. See
Gaussian
quadrature for the locations of three points and how to change the
interval from [-1,1] to [0,1/4]. Again, compute errors and compare (to
the methods in a as well).
- Review composite quadrature rules (Example 8.6 of the book) and
approximate the integral using composite midpoint and trapezoid rules by
breaking the interval into n subintervals. Try various values of n.
Again compute errors; at what value of n, roughly, are the composite
rules comparable to Gaussian quadrature with three points?
- (Extra Credit) Derive the composite Simpson's rule for arbitrary n
and apply it to approximate pi using the above formula. How does its
accuracy compare to the composite midpoint and trapezoid rules?
Marcus Schaefer
Last updated: March 5th, 2010.