The \(n \times n\) Hilbert matrix \(H^{(n)}\) (see page 512) defined by
$$
H_{i j}^{(n)}=\frac{1}{i+j-1}, \quad 1 \leq i, j \leq n
$$
is an ill-conditioned matrix that arises in solving the normal equations for
the coefficients of the least-squares polynomial (see Example 1 of Section
8.2).
a. Show that
$$
\left[H^{(4)}\right]^{-1}=\left[\begin{array}{rrrr}
16 & -120 & 240 & -140 \\
-120 & 1200 & -2700 & 1680 \\
240 & -2700 & 6480 & -4200 \\
-140 & 1680 & -4200 & 2800
\end{array}\right],
$$
and compute \(K_{\infty}\left(H^{(4)}\right)\).
b. Show that
$$
\left[H^{(5)}\right]^{-1}=\left[\begin{array}{rrrrr}
25 & -300 & 1050 & -1400 & 630 \\
-300 & 4800 & -18900 & 26880 & -12600 \\
1050 & -18900 & 79380 & -117600 & 56700 \\
-1400 & 26880 & -117600 & 179200 & -88200 \\
630 & -12600 & 56700 & -88200 & 44100
\end{array}\right],
$$
and compute \(K_{\infty}\left(H^{(5)}\right)\).
c. Solve the linear system
$$
H^{(4)}\left[\begin{array}{l}
x_{1} \\
x_{2} \\
x_{3} \\
x_{4}
\end{array}\right]=\left[\begin{array}{l}
1 \\
0 \\
0 \\
1
\end{array}\right]
$$
using five-digit rounding arithmetic, and compare the actual error to that
estimated in (7.25).