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)