According to MathWorld,
Pentagonal Triangular Number: A number which is simultaneously a pentagonal number $P_n$ and triangular number $T_m$. Such numbers exist when $$\frac{1}{2}n(3n-1)=\frac{1}{2}m(m+1).$$ Completing the square gives $$(6n-1)^2-3(2m+1)^2=-2.$$ Substituting $x=6n-1$ and $y=2m+1$ gives the Pell-like quadratic Diophantine equation $$x^2-3y^2=-2,$$ which has solutions $(x,y)=(5,3),(19,11),(71,41),(265,153), \ldots$.
However, it does not state how these solutions for $(x,y)$ were obtained.
I know that the solution $(5,3)$ can be obtained by observing that $1$ is both a pentagonal and a triangular number.
Does obtaining the other solutions simply involve trial-and-error? Or is there a way to obtain these solutions?

If you're a bit familiar with algebraic number theory:
$x^2 - 3y^2$ is the norm of the element $x + y\sqrt{3}$ in $\mathbb{Q}(\sqrt{3})$. Given the obvious element $1 + \sqrt{3}$ with norm $-2$, every other possibility differs by multiplication with an element of norm $1$. Dirichlet's unit theorem characterizes them: all powers of $2 + \sqrt{3}$ (up to $\pm 1$).
So the solutions are given by $\pm x \pm y\sqrt{3} = (1 + \sqrt{3})(2 + \sqrt{3})^n$ for $n \in \mathbb{Z}$.