弦截法的matlab实现.doc

《弦截法的matlab实现.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、弦截法的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

2、);p0=p1;p1=p2;k,p1,err,y=feval(‘f’,p1)if(err

3、(y==0),break,endend