在matlab中寻找曲线上离原点距离最近的点的坐标的思考

syms y
x=linspace(-2,2,200)
n=length(x)
g=[];        
for i=1:n
    a=x(i);
    t=vpa(solve(y/(1+a^2+y^2)-sin(a+cos(y)),y),3);
    g=[g t];
end

对于这个代码主要是想算出对应x的y，因为直接用get函数得到的x和y是按照顺序先后排列而不是一一对应。下面那个代码得出的结果是错误的。

syms x y
f=y/(x^2+y^