>> syms x;f=inline('sin(x)-x*(1/3)')
f =
Inline function:
f(x) = sin(x)-x*(1/3)
>> fplot(f,[-10,10])
>> grid
由圖像可以看出在[-4,-2]和[2,4]之間有兩個(gè)非零解,所以使用下面兩條命令
>> [xroot,y]=fsolve(f,[2,4],1e-5)
Optimization terminated:first-order optimality is less than options.TolFun.
xroot =
2.278862660075828 2.278862660076503
y =
1.0e-012 *
0 -0.663136212608606
>> [xroot,y]=fsolve(f,[-4,-2],1e-5)
Optimization terminated:first-order optimality is less than options.TolFun.
xroot =
-2.278862660076468 -2.278862660075828
y =
1.0e-012 *
0.629607477264926 0
注意這里最終結(jié)果xroot是最終的區(qū)間（精度為1e-5）,也就是說(shuō)方程的根在[2.278862660075828 2.278862660076503]之間,y是把該值代入f得到的函數(shù)值.