function?root=SimpleNewton（fabeps）
if（nargin==3）
????eps=1.0e-4;
end

f1=subs（sym（f）findsym（sym（f））a）;
f2=subs（sym（f）findsym（sym（f））b）;
if（f1==0）
????root=a;
end
if（f2==0）
????root=b;
end

if（f1*f2>0）
????disp（‘兩端點(diǎn)函數(shù)值乘積大于0!‘）;
????return;
else
????tol=1;
????fun=diff（sym（f））;
????fa=subs（sym（f）findsym（sym（f））a）;
????fb=subs（sym（f）findsym（sym（f））