Это следует из нулевых рангов соответствующих эллиптических кривых в диапазоне

.
Здесь два семейства эллиптических кривых

при

и

при

В обоих случаях кривые приводятся к форме Вейерштрасса с помощью Maple, а затем с помощью Pari
считаются ранги кривых.
Вот коды Pari для определения нулевых рангов этих кривых.
Код:
(20:18) gp > {for(N=2,50,
E=ellinit([0,-N^4 - 6*N^2 - 1,0,4*N^2*(N^2 + 1)^2,0]);
r=ellanalyticrank(E)[1];
if(r<1,print1(N,", ")))}
y=2, 3, 5, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 19, 23, 25, 26, 28, 29, 31, 37, 39, 42, 43, 48, 50
(20:38) gp > {for(N=2,50,
E=ellinit([0,N^4 - 6*N^2 + 1,0,-4*N^2*(N^2 - 1)^2,0]);
r=ellanalyticrank(E)[1];
if(r<1,print1(N,", ")))}
x=2, 3, 4, 5, 6, 7, 8, 12, 20, 26, 27, 28, 30, 32, 35, 36, 38, 39, 42, 44, 47
Теперь о целых значениях

. На самом деле их бесконечно много.
Процедуру вычисления сообщу позже. Она интересная и связана с уравнением Пелля.