方程的形式是
f=(y-y*x^2+x*(sqrt(y^2*x^2-y^2+b^2*m^2)))/(b*(y-y*x^2+x*(sqrt(y^2*x^2-y^2+n^2))))-a
已知参数a,b,m,n的值以及x的取值范围为（-1,0），现在想利用数值解析的方法求解y的值，不会写matlab程序，哪位大神帮忙看一下