qq_45687141 2019-11-04 23:04 采纳率: 0%
浏览 1317

用牛顿迭代法求方程的根c++,麻烦注释一下代码,搞不懂啥意思?

自定义一个函数,实现用牛顿迭代法求方程的根。
方程为ax^3+bx^2+cx+d=0。系数a, b, c, d的值依次为1, 2, 3, 4,由主函数输入。求x在1附近的一个实根。求出根后由主函数输出。
牛顿迭代法的公式是x = x0 - f(x0) / f '(x),精度为10^-5。

#include
#include
#include
using namespace std;
float fun1(float x)
{
return x*x*x + 2*x*x + 3*x + 4;
}
float fun2(float x)
{
return 3*x*x + 4*x + 3;
}

int main()
{
    float x , f1 ,f2 , d;
    x = 1;
    do
    {
        f1 = fun1( x);
        f2 = fun2( x);
        d = f1 / f2;
        x -= d;
    }
        while(fabs(d) > 1e-5);//这一句fabs没学过不懂,能换别的方式表达吗 

        printf("%.5f\n\n",x);

        return 0;
 }

    }
  • 写回答

2条回答 默认 最新

  • 豆丷 2019-11-05 09:15
    关注
    int main()
    {
        float x , f1 ,f2 , d;
        x = 1;//初值
        do
        {
            f1 = fun1( x);//f(x0)的值
            f2 = fun2( x);//f‘(x0)的值
            d = f1 / f2;//x的偏移值
            x -= d;//下一个x值
        }
            while(fabs(d) > 1e-5);//判断精度,不满足则继续循环,满足的跳出,fabs是取绝对值
    
            printf("%.5f\n\n",x);
    
            return 0;
     }
    
    评论

报告相同问题?

悬赏问题

  • ¥15 素材场景中光线烘焙后灯光失效
  • ¥15 请教一下各位,为什么我这个没有实现模拟点击
  • ¥15 执行 virtuoso 命令后,界面没有,cadence 启动不起来
  • ¥50 comfyui下连接animatediff节点生成视频质量非常差的原因
  • ¥20 有关区间dp的问题求解
  • ¥15 多电路系统共用电源的串扰问题
  • ¥15 slam rangenet++配置
  • ¥15 有没有研究水声通信方面的帮我改俩matlab代码
  • ¥15 ubuntu子系统密码忘记
  • ¥15 保护模式-系统加载-段寄存器