solve_Pell(): get the NO-th solution of Pell's equation: x^2 + d y^2 = n (n=+1 || -1).