搜索
弦截法的matlab实现编辑Secant.m文档：Funtion[p1,err,k,y]=secant(f,p0,p1,delta,max1)%f是非线性函数%p0,p1是初始值%delta是给定允许误差%max1是迭代次数的上限%p1是所求得的方程的近似解%err是p1-p0的误差估计%k是所需要的迭代次数%y=f(p1)K=0,p0，p1,feval(‘f’,p0),feval(‘f’,p1)fork=1:max1p2=p1-feval(‘f’,p1)(p1-p0)/(feval(‘f’,p1)-feval(‘f’,p0));err=abs(p2-p1);p0=p1;p1=p2;k,p1,err,y=feval(‘f’,p1)if(errdelta)|(y==0),break,endend