More specifically, 25 + 36(R^2-N) is the problem; my current algorithm

is as follows:

xSet = Set new.

1 to: 100 do: [:x | xSet add:(x*x)\\100]. {gives all square endings}

R = n isqrt2 + 1. {makes R greater than N}

while false:[

[t = R^2 -N.

q := 25 + 36*t. xSet includes:q\\100] whilefalse:[R := R+1].

p := q isqrt2.

P*p = q ifTrue:[^true]

ifFalse:[R:=R+1]]

thanks in advance

Ron