数值计算方法作业-二分法

编程实现，用二分法求$f(x)=x^3+4x^2-10$在区间$[1,2]$内的一个实根，要得到精确到小数点后三位位精确值的近似
注意：编程的循环次数

易证该方程在区间内单调递增，所以以此来确定二分方向

C++实现：

#include <bits/stdc++.h>
/****************
author:WANG_zibi
time:2020/2/25
title:二分法求方程近似值
course:数值计算方法
****************/
using namespace std;
const int N = 1e6+5;
typedef long long ll;
typedef unsigned long long ull;
const int INF= 0x3f3f3f3f;
#define f(i, a, b) for (int i = (a); i <= (b); ++i)
const double eps = 1e-4;//精度
double func(double x) //函数
{
    return x*x*x+4*x*x-10.0;
}
int main()
{
    double left,right,mid,answer,sum;//左右区间，答案等变量生命
    left=1,right=2；
    mid=(left+right)/2.0;
    sum=func(mid);
    answer=mid;
    while(abs(sum)>eps)
    {
        if(sum>0)
        {
            right=mid;
            mid=(right+left)/2.0;
            sum=func(mid);
            answer=mid;
        }
        else{
            left=mid;
            mid=(right+left)/2.0;
            sum=func(mid);
            answer=mid;
        }
    }
    
    cout<<"最终的答案为："<<answer<<'\n';
    return 0;
}

MATLAB实现：

%author:WANG_zibi
%date:2020/2/25
%course:数值计算方法
%name:二分法求方程近似解
%方程
l=1.0;
r=2.0;
x=(l+r)/2.0;
eps=0.0001;
f=x^3+4*(x^2)-10;
while abs(f)>eps
    if(f>0)
        r=x;
        x=(l+r)/2.0;
        f=x^3+4*(x^2)-10;
    end
    if(f<0)
        l=x;
        x=(l+r)/2.0;
        f=x^3+4*(x^2)-10;
    end
end
ans=x