非線形方程式  戻る

x-cos(x)=0 を解くには,初期値を定めて(ここでは0.7)

deff('y=f(x)','y=x-cos(x)');fsolve(0.7,f)
ans =

0.7390851

と求まる。 初期値はグラフを描いて,x軸との交点近くの x座標を採ればよい
x=-5:0.05:5;y=x-cos(x);plot(x,y);plot(x,0)

exp(x)-x=2 を解くとき,まずグラフを描いて
x=-5:0.05:2;y=exp(x)-x-2;plot(x,y);plot(x,0)
初期値を1とすれば
deff('y=f(x)','y=exp(x)-x-2');fsolve(1,f)
 1.1461932 が得られる。
また−1とすれば - 1.8414057が得られる。
初期値を 0 とすると解は 0 と表示されるが,この場合 f'(0)=0 となるから解なしのはず。

 整方程式の場合はもっと簡単にできる。係数が[ ]内の数で与えられた昇べきの整方程式(2−3x+x2=0)の解を求めるには

p=poly([2 -3 1],'x','coeff');x=roots(p)
この方法で虚の解も求まる。