The matrix trace and determinant of the matrices in part (i) are
Since and , the extra conditions force and . Hence
Both choices are real, commute with , and satisfy the required trace and determinant conditions.
Let . The commutativity condition gives
Equality is equivalent to and . Therefore
Each displayed matrix commutes with , since it is a polynomial in . The centralizer of a nonscalar two by two matrix explains why this two-parameter form is typical.
The two orthogonal directions and are preserved by this symmetric matrix. Direct multiplication, using the hyperbolic function identities , gives the eigenvalue–eigenvector pairs
Each eigenvector can be multiplied by any nonzero scalar. For these give the two one-dimensional eigenspaces. When , the matrix is the identity: the only eigenvalue is , with multiplicity two, and every nonzero vector is an eigenvector.
The vector is an eigenvector with eigenvalue . The complementary subspace of vectors with first two components equal is preserved by : a vector maps to . There, the eigenvalue equations reduce to the matrix , whose characteristic polynomial is . For , ; for , .
Normalising by the Euclidean norm, a complete orthonormal basis of eigenvectors is
Sign changes of any of these eigenvectors are equally valid. The eigenvalues also check part (a): their sum is and the sum of their squares is .
Three vectors form a basis of three-dimensional space exactly when their scalar triple product is nonzero:
It is the determinant of the matrix with these columns, so this condition is equivalent to linear independence and spanning the whole space.
Take the dot product of the expansion of with . The terms containing and vanish, leaving . Cycling the same argument yields
The three cross-product numerators use the reciprocal basis.
In the nonsingular case the unique solution is
Geometrically this is a single selected point of the straight line through parallel to .
When and , the solution set is that entire straight line, . The remaining case has empty solution set. In matrix terms the equation is , a rank-one matrix update of the identity matrix.
The displayed matrix is not orthogonal: its first column has squared Euclidean norm three, whereas every column of an orthogonal matrix must have squared length one.
For the general assertion, put . In the Einstein summation convention, orthogonality gives and . The entries of satisfy
This is exactly , proving that the matrix product is orthogonal. The Kronecker delta contracts the two middle indices without changing the multiplication order of the matrices.
For , the matrix is invertible, so there is exactly one solution. To exhibit it, subtract the first linear equation from the third to obtain , then subtract twice the second from the first to obtain . Back-substitution gives
At , the first and third linear equations have identical left-hand sides but right-hand sides one and three, so there is no solution.
At , subtracting the first from the third forces , while the first minus twice the second forces . These contradict each other, so there is again no solution. Thus there is one solution for , no solution for either exceptional value, and no parameter value with infinitely many solutions.
The coefficient matrix has determinant , so it is an invertible matrix. Its homogeneous linear system therefore has only the zero solution:
One can also subtract the first equation from the second to get , then subtract the first from the third to get , giving and finally . Thus there are no nonzero solutions.
Substituting and into the matrix gives
Its row sums are , checking that for every real .
Differentiate the first ordinary differential equation and substitute the second:
The original first equation also gives , so
Therefore the requested coefficients are
They are the trace and determinant of the coefficient matrix. This derivation never divides by and is valid even if . In that degenerate case the scalar second-order equation is a necessary consequence, while the original constraint must also be retained when reconstructing the coupled system.
Suppose . The given quadratic form bound then gives
so . Thus the kernel of is trivial. A square matrix with trivial kernel is an invertible matrix, proving the assertion. No symmetry assumption on is needed.
The first coordinate direction is already an eigenvector of with eigenvalue . On the last two coordinates, the matrix is , with characteristic polynomial
Thus the three eigenvalues are . For , the last two components are equal; for , they are opposite. Normalising by the Euclidean norm gives the complete list
Each sign-reversed vector is an equally valid choice. The three displayed eigenvectors are mutually orthogonal, as expected for a real symmetric matrix with distinct eigenvalues.
The matrix is
so
Here , so with . This is a Householder transformation: reflection in the plane . The component parallel to the normal vector changes sign, while the perpendicular component is fixed.
The second matrix acts by
Applying from part (e) then gives
The Householder transformation reflects in the plane . Its normal vector is perpendicular to , so the two reflecting planes meet at a right angle along the axis. Their composition is a rotation through about the axis, as is also immediate from the displayed rotation matrix.