Omitting solutions that are obtained from smaller solutions by multiplying by powers of 10, the squares of these numbers satisfy the

- condition
1. (105,145), (3375,4625), (14025,17225), (326625,454625),

(10846875,14753125), (43708125,53948125), ...

2. (45,55), (144375,175625), (463171875,560828125), ...

7. (2824483699753370361328125,2996282391593370361328125), ...

Here is how to find them. We have (y+x)*(y-x) = 10^n, and so we must have
{y+x, y-x} as {5^m*10^a, 2^m*10^b} in some order. It is also necessary (and
sufficient) that y/x lies in the interval [sqrt(3__/2),sqrt(2)__?, or equivalently
that (y+x)/(y-x) lies in __[3+sqrt(8),5+sqrt(24)__?] = __[5.82842...,9.89897...__?].
Thus we need to make (5/2)^m*10^(a-b), or its reciprocal, in this range.
For each m there is clearly at most one power of 10 that will do. m=2,a=b
gives (105,145); m=3,b=a+2 gives (3375,4625), and so on.

There are infinitely many non-equivalent solutions, because log(5/2) / log(10) is irrational.

One can use exactly the same argument to find squares whose initial 2 can be replaced by a 3, of course, except that the range of (y+x)/(y-x) changes.