4. The Runge Example. Let f(x) = 1 / (1 + 25x^2), x in [-1, 1]. Using your Barycentric Formula code (Prob. 3) and (4) and (5) below, evaluate and plot the interpolating polynomial pn of f corresponding to (a) The equidistributed nodes xj = -1 + j(2/n), j = 0, ..., n for n = 4, 8, and 12. (b) the Chebyshev nodes xj = cos(j*pi/n), j = 0, ..., n for n = 4, 8, 12, and 100. As seen in class, for equidistributed nodes one can use the barycentric weights lambda_j = (-1)^j * (n choose j), j = 0, ..., n, and for the Chebyshev nodes we can use lambda_j = 1/2*(-1)^j for j = 0 or j = n, (-1)^j for j = 1, ..., n - 1. Make sure to employ (4) and (5) in your Barycentric Formula code for this problem. To plot the corresponding pn evaluate this at sufficiently large number of points ne as in Prob. 2. Note that your Barycentric Formula cannot be used to evaluate pn when x coincides with an interpolating node! Plot also f for comparison. (c) Plot the error en = f - pn for (a) and (b) and comment on the results. (d) Repeat (a) for f(x) = e^(-x^2) for x in [-1, 1] and comment on the result.