利用fsolve可以求解非线性方程的数值解,但是只有最后的结果,没有每步迭代的x取值。
在求解之前设置迭代显示的选项
options=optimset('Display','iter');
输出窗口中显示如下:
Norm of First-order Trust-region
Iteration Func-count f(x) step optimality radius
0 2 2.61482e-007 0.000171 1
1 3 2.61482e-007 0.00152725 0.000171 1
2 5 1.51338e-007 0.000381811 0.000109 0.000382
3 6 1.51338e-007 0.000954528 0.000109 0.000955
4 8 8.78862e-008 0.000238632 8.52e-005 0.000239
5 10 7.98377e-008 0.00059658 8.75e-005 0.000597
6 11 7.98377e-008 0.00059658