谁能用Matlab帮我算下13*x^2+1=y^2的整数解

syms y
x=0:1000;[m,n]=size(x);
y=sqrt(13.*x.^2+1);
for i=1:n
if max(1:y(i))==y(i)
X=x(i),Y=y(i)
end
end
运行结果：
X =

0

Y =

1

X =

180

Y =

649
x在0到1000内就两个正解
所以x在[-1000,1000],y在[-1000,1000]内
有解:
(0,1),(180,649),(-180,649),(-180,-649),(0,-1)

(0,±1),
(±180,±649),
共有6个解.对应6种组合.
其实不止-1000～1000内只有6个解,整个实数范围内也只有这6个解.