fsolve函数用于数值求解非线性方程。
coef=[1 2 2 1];
option=optimset('TolX',1e-8,'TolFun',1e-8);
root=fsolve(@(x) polyval(coef,x),0,option);
这段程序求解coef数组定义的代数方程 x^3+2x^2+2x+1=0 在 x=0 附近的解。optimset函数设置收敛条件。