第一步求解出的y有多组解，你可以运行一下看结果：
y=solve('(10*y+2.4448e-04*y)^2+(-y*(x-101.4)/0.1+3*(6.2723e+07)*(y^3)/(8*101.4)+y*(101.4^2)/(2*(101.4^2+25)))^2-(1/1.2168)^2','y')
我取的是第一组解（fy=matlabFunction(y(1));），因为我们版本不一样，所以你的排在第一组的解与我的不同，于是出现这样的结果，如果我换第三组解，图形就和你一样了。但是这组解是复数解，画图时只取了实部，应该是不对的。
不知你对这个y有没有要求，比如正负？
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024