So this is
[tex]f_h(a+ bx+ cx^2+ dx^3)= \begin{bmatrix} a+ b+ c+ hd & b+ c \\ -b- c- hd & hb \end{bmatrix}[/tex]
from the space of polynomials of degree 3 (or less) to the space of 2 by 2 matrices. And the problem is to find all possible values of h such that what?
As says says, since this is from a 4 dimensional vector space to a four dimensional vector space, it would be written
as a 4 by 4 matrix. Taking the "obvious" bases for the two spaces, [itex]\{1, x, x^2, x^3\}[/itex] and [itex]\{\begin{bmatrix}1 & 0 \\ 0 & 0 \end{bmatrix}, \begin{bmatrix}0 & 1 \\ 0 & 0 \end{bmatrix}, \begin{bmatrix}0 & 0 \\ 1 & 0 \end{bmatrix}, \begin{bmatrix}0 & 0 \\ 0 & 1 \end{bmatrix}\}[/itex], this maps [itex]1= 1+ 0x+ 0x^2+ 0x^3[/itex] (a= 1, b= c= d= 0) to [itex]\begin{bmatrix}1 & 0 \\ 0 & 0 \end{bmatrix}[/itex], it maps [itex]x= 0+ 1x+ 0x^2+ 0x^3[/itex] (b= 1, a= c= d= 0) to [itex]\begin{bmatrix}1 & 1 \\ -1 & h\end{bmatrix}= \begin{bmatrix}1 & 0 \\ 0 & 0 \end{bmatrix}+ \begin{bmatrix}0 & 1 \\ 0 & 0 \end{bmatrix}- \begin{bmatrix}0 & 0 \\ 1 & 0 \end{bmatrix}+ h\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix}[/itex], it maps [itex]0+ 0x+ 1x^2+ 0x^3[/itex] (c= 1, a= b= d= 0) to [itex]\begin{bmatrix} 1 & 1 \\ -1 & 0\end{bmatrix}= \begin{bmatrix}1 & 0 \\ 0 & 0 \end{bmatrix}+ \begin{bmatrix}0 & 1 \\ 0 & 0 \end{bmatrix}- \begin{bmatrix}0 & 0 \\ 1 & 0 \end{bmatrix}[/itex], and maps [itex]0+ 0x+ 0x^2+ 0x^3[/itex] (d= 1, a= b= c= 0) to [itex]\begin{bmatrix} h & 0 \\ -h & 0\end{bmatrix}= h\begin{bmatrix} 1 & 0 \\ 0 & 0 \end{bmatrix}- h\begin{bmatrix}0 & 0 \\ 1 & 0 \end{bmatrix}[/itex]
So the matrix representing this linear transformation, in these bases, is
[tex]\begin{bmatrix} 1 & 1 & 1 & h \\ 0 & 1 & 1 & 0 \\ 0 & -1 & -1 & -h \\ 0 & h & 0 & 0 \end{bmatrix}[/tex].
Now, again, exactly what is it you want to do with, or know about, this matrix?