Como señala @ΜάρκοςΚαραμέρης, las soluciones aquí son coordenadas enteras (primas) en una curva elíptica, y la teoría general (Thm. de Siegel) nos dice que para un cúbico irreducible, hay como mucho muchos puntos enteros. Obtener el número real de éstos puede ser a veces un cálculo tedioso, pero en el presente caso podemos obtener límites estrechos con métodos elementales (por ejemplo, completando el cuadrado).
Como se ha señalado en los comentarios de Keith Backman y otros, la búsqueda de coordenadas primarias $(p,q)$ (que asumimos positivo por simplicidad) de:
$$ p(p^2-p-1)=q(2q+3) $$
se reduce a $p^2-p-1 = qk$ y $2q+3=pk$ (así que $k \gt 0$ ) donde:
$$ 2p^2 - (k^2 +2)p + (3k-2) = 0 $$
Una solución integral $p$ de esta ecuación cuadrática requiere $k$ que hace del discriminante un cuadrado perfecto, es decir:
$$ m^2 = (k^2 + 2)^2 - 8(3k-2) $$
Para limitar estas posibles $k$ basta con notar que para que sea suficientemente grande $k$ el discriminante se encuentra estrictamente entre los dos cuadrados consecutivos:
$$ (k^2 + 1)^2 \lt (k^2 + 2)^2 - 8(3k-2) \lt (k^2 + 2)^2 $$
La desigualdad de la mano derecha se mantiene para $k \ge 1$ y para probar la desigualdad de la mano izquierda la reescribimos como:
$$ 8(3k-2) \lt (k^2 + 2)^2 - (k^2 + 1)^2 = 2k^2 + 3 $$
$$ 0 \lt 2k^2 - 24k + 19 = 2(k-6)^2 - 53 $$
Por la inspección, esto último es cierto para todos $k \ge 12$ así que nuestra búsqueda de soluciones sólo necesita cubrir $k = 1, \ldots , 11$ . Por supuesto, este suelo ya está fregado, produciendo la única solución primaria (positiva) ( $k=5, p=13, q=31$ ).