$$\begin{aligned} 0 &= -8x^2 + 8xy + 17y^2 - 34xz - 2yz - 28z^2 - 10xw - 9yw - 18zw + 2w^2, \\ 0 &= 4x^3 - 6x^2 y - 6x y^2 + 12x^2 z + 6xyz + 24y^2 z - 12x z^2 - 24z^3 + 2x^2 w + 7xyw \\ &\qquad + 4y^2 w + 4xzw - 13yzw - 8z^2 w - 20x w^2 - 3z w^2 - 12w^3 \end{aligned}$$
has real multiplication by the maximal order of $\mathbf{Q}(x)/(x^4 - x^3 - 3x^2 + x + 1)$.The first step to show that, under Langlands, it corresponds to a specific Hilbert modular form $f$, i.e., $J_{\mathbf{Q}(\sqrt{3})} \sim A_f$. We used this in a recent project, where we show that the 2-isogeny field of $A_f$ solves the inverse Galois problem for $\operatorname{PSL}_2(\mathbf{F}_{16}) \rtimes C_2 \simeq \texttt{17T7}$.