4matlab正割线法求解非线性方程的数值解

正割线法求解非线性方程的数值解：M——文件（secant.m）functionxr=secant(fun,x0,x1,eqs)%secant.m函数为用于正割线法求解非线性方程的数值解法%fun是给定的非线性方程、%xox1是给定的初始值%eqs是精度ifnargin==3eqs=1e-6;endf0=feval(fun,x0);f1=feval(fun,x1);whileabs(x0-x1)eqsx2=x1-f1*(x0-x1)/(f0-f1);x0=x1;f0=f1;x1=x2;f1=feval(fun,x1);endxr=x2;命令窗口m文件：（secant1.m）clearallclcfplot('exp(x)-x-5',[-5,5]);gridonholdonfun=inline('exp(x)-x-5');xr=secant(fun,3.5,3);fxr=feval(fun,xr);xrfxr结果：xr=1.9368fxr=3.5527e-015-5-4-3-2-1012345-20020406080100120140