谭浩强c第四版218页12题~

用牛顿迭代法求根。方程ax^3+bx^2+cx+d=0,系数a,b,c,d的值依次为1,2,3,4,由主函数输入。求x在1附近的一个实根。求出根后由主函数输出。

不懂牛顿迭代法~求解~

#include<stdio.h>
#include<math.h>
int main()
{
    double solut(double ,double ,double ,double );
    double a,b,c,d;
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
    printf("%.2f",solut(a,b,c,d));;
    return 0;
}
double solut(double a,double b,double c,double d)
{
    double x=1,x1=2,f,f1;
    while(fabs(x1-x)>=0.00000001)
    {
        f=((a*x+b)*x+c)*x+d;
        f1=(3*a*x+2*b)*x+c;
        x=x1;
        x1=x-f/f1;
    }
    return x1;
}

啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦啦代码不是我滴，百度滴，，，，，

[此贴子已经被作者于2017-1-20 22:02编辑过]

满分~结~

好像结贴过早了~

突然发现结贴过早了还有些细节没有消化~虽然编程还可以可惜高数不过关~

好像和求导有关

收到的鲜花

• 九转星河 于 2017-01-20 23:46 送鲜花  10朵   附言：看得出来~

#include <stdio.h>
int main(void)
{
    double a = 1 , b = 2 , c = 3 , d = 4;
    double jump = 0.000001;    //这个是精度，可自行修改
    double x = 1+jump;
    double y = 0;

    do {
        x = x - jump;    //x从1开始以jump的精度缩小
        y = a*x*x*x+b*x*x+c*x+d;    //公式结果赋予y
        printf(" x的值等于：%lf，结果就等于：%lf\n",x,y);
    }while(y>jump || x>jump);    //如果y小于精度，则足够少，可以结束并打印结果，如果x小于精度，说明求不出结果，结束。
    
    if (x < jump){    //判断一下，怕数据越界
        printf("\n 最后需要设置x为：%lf，公式的结果就等于：%lf。\n",x,y);    //最后输出结果。
    }else{
        printf("\n 最后x超出精度也求不出目标精度的结果。\n");
    };
    return 0;
}

不知有否理解错误题目，请大家检阅。

收到的鲜花

• 九转星河 于 2017-01-21 10:44 送鲜花  5朵   附言：……

函数：f(x)=x*x*x+2x*x+3x+4
导函数：f'(x)=3x*x+4x+3
X(n+1) = Xn - f(Xn)/f'(Xn)

收到的鲜花

• 九转星河 于 2017-01-21 10:46 送鲜花  15朵   附言：谢谢吹版解答~

这个看上虽然很好理解~但看上去用二分法更快~牛顿迭代法应该不会运算那么多步骤~不知我有没有理解错题目~
        还是谢谢你的解答~

谢啦~我还是去看看牛顿迭代法怎么用~