I want to explicitly determine the subfield of GF$(16)$=$\mathbb{F}_{16}$ with 4 elements. I know this to be $\mathbb{F}_4$ (or isomorphic to it) since a field of prime power (in this case $2^4$) has a unique subfield of order $2^m$ iff $m$ divides $4$. The method I used seemed quite cumbersome.
I determined the roots of $Y^4-Y$ in $\mathbb{F}_{16}$ since $\mathbb{F}_4$ is the decomposition field of the polynomial $Y^4-Y$. I found that $$Y^4-Y=Y(Y-1)(Y^2+Y+1),$$ and so $0,1$ are already in $\mathbb{F}_4$ as a subfield of $\mathbb{F}_{16}$. Also $$\mathbb{F}_4\cong\dfrac{\mathbb{Z}_2[X]}{\langle X^2+X+1\rangle},$$ which means that $X^2=X+1$ in $\mathbb{F}_4$ and $$\mathbb{F}_{16}\cong \dfrac{\mathbb{Z}_2[X]}{\langle X^4+X+1\rangle},$$ which means $X^4=X+1$ in $\mathbb{F}_{16}$ and so elements of $\mathbb{F}_{16}$ are of the form \begin{equation} a_3X^3+a_2X^2+a_1X+a_0 \end{equation} with $a_i \in \mathbb{Z}_2$ for $i=0,1,2,3$. Next, I substituted $a_3X^3+a_2X^2+a_1X+a_0$ into $Y^2+Y+1$ and used the properties $X^4=X+1$, $X^2=X+1$ and found, after a long calculation, something of the form $$b_1X+b_0,$$ with $b_1,b_0\in\mathbb{Z}_2$ some sum of $a_i$'s and powers of them. This means the two remaining elements of $\mathbb{F}_4$ are $X+1$ and $X$ and so I found that $\mathbb{F}_4=\{0,1,X,X+1\}$ as a subfield of $\mathbb{F}_{16}.$
First of all, is this method correct? I believe it is since I found a correct answer for $\mathbb{F}_4$.
Second, is there a more simple solution to this? I was hoping for a more generic way of finding such subfields.