“…Below we have enclosed the Sage code that was used to determine the small solutions to T n1 −T n2 = b m1 −b m2 in Section 6, as well as the code that was used to determine the bounds for the reduction rounds in Table 3. if b^x * (b^y -1) == diff: m2 = x m1 = x + y c = t2 -b^x print("T_{",n1,"} -T_{",n2,"} =",b,"^{",m2,"} (",b,"^{",y,"}-1), \\quad c&=", c, "= T_{",n1,"} -",b,"^{",m1,"} = T_{",n2,"} -",b,"^{",m2,"}; \\\\") # reduction steps alpha = n(solve(x^3 -x^2 -x -1 == 0, x, solution_dict=True) [2][x], digits=2000) a = 1/(-alpha^2 + 4*alpha -1)…”