# http://e...content-available-to-author-only...a.org/wiki/Pell%27s_equation
http://e...content-available-to-author-only...a.org/wiki/Continued_fraction#Best_rational_approximations
 
from fractions import *
from math import *
 
 
def PellsEquation(convergent, n):
 return (convergent.numerator**2) - n*(convergent.denominator**2)
 
# The sequence of convergents for the square root of seven are
L = [Fraction(2,1), Fraction(3,1), Fraction(5,2), Fraction(8,3)]
 
for f in L:
 print(PellsEquation(f,7))
 
print(Fraction(str(sqrt(7))).limit_denominator())

IyBodHRwOi8vZS4uLmNvbnRlbnQtYXZhaWxhYmxlLXRvLWF1dGhvci1vbmx5Li4uYS5vcmcvd2lraS9QZWxsJTI3c19lcXVhdGlvbgpodHRwOi8vZS4uLmNvbnRlbnQtYXZhaWxhYmxlLXRvLWF1dGhvci1vbmx5Li4uYS5vcmcvd2lraS9Db250aW51ZWRfZnJhY3Rpb24jQmVzdF9yYXRpb25hbF9hcHByb3hpbWF0aW9ucwoKZnJvbSBmcmFjdGlvbnMgaW1wb3J0ICoKZnJvbSBtYXRoIGltcG9ydCAqCgoKZGVmIFBlbGxzRXF1YXRpb24oY29udmVyZ2VudCwgbik6CiByZXR1cm4gKGNvbnZlcmdlbnQubnVtZXJhdG9yKioyKSAtIG4qKGNvbnZlcmdlbnQuZGVub21pbmF0b3IqKjIpCgojIFRoZSBzZXF1ZW5jZSBvZiBjb252ZXJnZW50cyBmb3IgdGhlIHNxdWFyZSByb290IG9mIHNldmVuIGFyZQpMID0gW0ZyYWN0aW9uKDIsMSksIEZyYWN0aW9uKDMsMSksIEZyYWN0aW9uKDUsMiksIEZyYWN0aW9uKDgsMyldCgpmb3IgZiBpbiBMOgogcHJpbnQoUGVsbHNFcXVhdGlvbihmLDcpKQoKcHJpbnQoRnJhY3Rpb24oc3RyKHNxcnQoNykpKS5saW1pdF9kZW5vbWluYXRvcigpKQ==