弦截法的matlab实现(共1页).doc

弦截法的matlab实现编辑Secant.m文档：Funtionp1,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)for k=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