MATLAB学习:方程求根,二分法

 在期间[0,1]内用二分法求方程e^x+10*x-2=0的近似根,要求误差不超过0.5*10^(-3)。

实验结果:x=0.09033

代码示例:

function x = agui_bisect(fname,a,b,e)
%  fname为函数名,a,b为区间端点,e为精度
%  二分法
fa=feval(fname,a);
fb=feval(fname,b);
if fa*fb>0
    error('两端函数值为同号');end
k=0;
x=(a+b)/2;
while(b-a)>(2*e)
    fx=feval(fname,x);
    if fa*fx<0
        b=x;
    else
        a=x;
        fa=fx;
    end
    k=k+1;
    x=(a+b)/2;
end

输入内联函数:

 fun=inline('exp(x) +10*x - 2')

调用agui_bisect函数:

x = agui_bisect(fun,0,1,0.0005)

测试结果:

f8cb03da230540bfa32c5dd5ad96ae3e.png

 

 

  • 1
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

晞微

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值