function?root?=?Montecarlo（fBx0eps）
format?long;
if（nargin==3）
????eps=1.0e-4;
end
Fx?=?subs（sym（f）findsym（sym（f））x0）;

while?abs（Fx）>eps
????Fx?=?subs（sym（f）findsym（sym（f））x0）;
????Fx1?=?abs