I am running this code, unfortunately I have no idea that why it is asking me a yes or no question.
Normu(D) = my(k = bnfinit(y^2-D, 1), u = bnfunits(k)[1][1]); nfeltnorm(k, u);
{
L = List();
forprime (q = 5, 100,
if (q%8==5,
N2 = Normu(q);
forprime(p = q+1, 100,
if(p%8==5,
N1 = Normu(2*p);
N3 = Normu(2*p*q);
istput(L, [[q,p], [N1, N2, N3]])))));
}
Unfortunately I can not undrestand what PARI means by 'showing all possibilities'and asks me to type(y or n).
I tried to change commands, and change the default PARI size, but not successful. I know the problem is not in the function defined above. I appreciate your help.
Update: I wrote all the code in one line and it is working but it doesn't consider internal 'if's . It is not clear to me why it should work with this change.
Normu(D) = my(k = bnfinit(y^2-D, 1), u = bnfunits(k)[1][1]); nfeltnorm(k, u);
{
L = List();
forprime (q = 5, 100, if (q%8==5, N2 = Normu(q); forprime(p = q+1, 100, if(p%8==5, N1 = Normu(2*p); N3 = Normu(2*p*q); listput(L,[[q,p], [N1, N2, N3]])))));
}