According to Example 3.39, the $n \times n$ Hilbert matrix $H_n$ is positive definite, and hence we can apply the Conjugate Gradient Method to solve the linear system $H_n \mathbf{u}=\mathbf{f}$. For the values $n=5,10,30$, let $\mathbf{u}^* \in \mathbb{R}^n$ be the vector with all entries equal to 1 .
(a) Compute $\mathbf{f}=H_n \mathbf{u}^*$. (b) Use Gaussian Elimination to solve $H_n \mathbf{u}=\mathbf{f}$. How close is your solution to $\mathbf{u}^*$ ? ${ }^n$ (c) Does pivoting improve the solution in part (b)?
(d) Does the conjugate gradient algorithm do any better?