数值分析程序纪要(MATLAB)

一、求方程的根

    (1)二分法

图片发自简书App

function xc=bisect(f,a,b,tol)

if sign(f(a))*sign(f(b)) ＞＝0

    error('f(a)f(b)<0 not satisfied')

end

fa=f(a);

fb=f(b);

while (b-a)/2＞tol

    c=(a+b)/2;

    fc=f(c);

    if fc==0

    break

end

if sign(fc)*sign(fa)<0

  b=c;fb=fc;

else

  a=c;fa=fc;

    end

end

xc=(a+b)/2;

定义函数f=@(x) 表达式，

xc=bisect(f，a，b，tool)


＊注记:

 

最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容