for y=0:1000
    x=sqrt(7*y^2+1)
    if x==floor(x) then
        afficher([x,y])
    end
end