美文网首页
数值分析程序纪要(MATLAB)

数值分析程序纪要(MATLAB)

作者: 叫酿么 | 来源:发表于2018-03-29 16:17 被阅读0次

一、求方程的根

    (1)二分法

数值分析程序纪要(MATLAB)

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)

*注记:

 

相关文章

网友评论

      本文标题:数值分析程序纪要(MATLAB)

      本文链接:https://www.haomeiwen.com/subject/wvamcftx.html