The function wilkpoly in the software distribution computes the MATLAB representation of the Wilkinson polynomial  Compute the roots of the polynomial after a perturbation of −23 in the coefficient of...

The function wilkpoly in the software distribution computes the MATLAB representation of the Wilkinson polynomial

 Compute the roots of the polynomial after a perturbation of

⁻²³
in the coefficient of

¹⁹.

 Using the definition of

<sub>ai

k</sub>
, compute the sensitivity of each root

_k
 relative to a change in the coefficient of

¹⁹. Use the function polyder discussed in Problem 10.33. Do your results agree with the errors in the roots you computed in part