==> arithmetic/pell.s <== x=1 y=0 x=1151 y=120 x=2649601 y=276240 etc. Each successive solution is about 2300 times the previous solution; they are every 8th partial fraction (x=numerator, y=denominator) of the continued fraction for sqrt(92) = [9, 1,1,2,4,2,1,1,18, 1,1,2,4,2,1,1,18, 1,1,2,4,2,1,1,18, ...] Once you have the smallest positive solution (x1,y1) you don't need to "search" for the rest. You can obtain the nth positive solution (xn,yn) by the formula (x1 + y1 sqrt(92))^n = xn + yn sqrt(92). See Niven & Zuckerman's An Introduction to the Theory of Numbers. Look in the index under Pell's equation.