函数求零点(二分法,牛顿法)

二分

找两个值:l(low)与h(high)。使得f(l)<0且f(h)>0。
函数连续性,x0在l与h之间。
设m=l+h2
,判断f(m)的正负性:
f(m)>0⇔x0在l与m之间
f(m)=0⇔x0=m
f(m)<0⇔x0在m与h之间

牛顿法

迭代公式
这里写图片描述

#include<bits/stdc++.h>
using namespace std;
double f(double x,int c)//函数f(x) = x^3-C
{
  return x*x*x-c;
}
double f_d(double x) //函数的导数
{
  return 3*x*x;
}
double mySqrt(int x)//二分法求零点
{
     double low = 0;         //f(low)<0
     double high = INT_MAX;  //f(high)>0
     double mid;
     while(high-low>0.000001) 
     {
        mid = low + (high-low)/2;
        if(f(mid,x)>0) high  = mid;
        else
           low = mid;  
     }
     return mid;
}

double mySqrt2(int x)
{
  double xnew,xnow;
  xnow = x;
  while(1)
  {
    xnew = xnow - f(xnow,x)/f_d(xnow);
    if(xnow-xnew<0.000001)break;
    xnow =xnew;
  }
  return xnow;
}
int main()
{
    while(1)
    {
        int  x;
        cin>>x;
      cout<<mySqrt(x)<<endl;
      cout<<mySqrt2(x)<<endl;
    }
    return 0;
}
  • 1
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值