[已解决]牛顿迭代法求根，代码运行错误，找不到原因求帮助呀！

798487012

用牛顿迭代法和二分法求根。方程为ax3+bx2+cx+d=0,系数a,b,c,d的值依次为2，-4，3，-6，由主函数输入。用牛顿迭代法求x在1附近的根；

#include<stdio.h>
#include<math.h>

double f(double x)
{        
        return 2*x*x*x+(-4)*x*x+3*x-6;
}

double f1(double x)
{        
        return 3*2*x*x+2*(-4)*x+3;
}


double main()
{

        double x=1,x1;
        x1=x-f(x)/f1(x);
        
        while (fabs(x1-x)>=1e-6)
        {
        
                x=x1;
                x1=x-f(x)/f1(x);
        }        
        return x1;
        printf("%f\n",x1);
        return 0;

最佳答案

月排行榜 / 总排行榜

claws0n

2018-9-28 22:19:55

798487012 发表于 2018-9-28 18:14
改成这样了，貌似可以了

#include <stdio.h>
#include <math.h>

double f(double x)
{        
    return 2*x*x*x+(-4)*x*x+3*x-6;
}

double f1(double x)
{        
    return 3*2*x*x+2*(-4)*x+3;
}

int main()
{
        double x,dx = 1;
        printf("输入一个估计值");
        scanf("%lf", &x);
        
        while (fabs(dx) >= 1e-7)
        {
                dx = f(x)/f1(x);
                x = x - dx;
        }
        
        printf("%f\n", x);
        return 0;

}

你不可以写 x1 = 6; 因为如果你输入 x = 6, fabs(x1-x) == 0，就不会进入循环了
没有加好友，呼叫是无效的 ^_^

跳转到最佳答案楼层

798487012

@claws0n

798487012

改成这样了，貌似可以了

#include<stdio.h>
#include<math.h>

double f(double x)
{        
        return 2*x*x*x+(-4)*x*x+3*x-6;
}

double f1(double x)
{        
        return 3*2*x*x+2*(-4)*x+3;
}


double main()
{        double x,x1;
     printf("输入一个数");
        scanf("%lf",&x);
        
        x1=6;
        
        while (fabs(x1-x)>=1e-6)
        {
        
                x=x1;
                x1=x-f(x)/f1(x);
        }        

        printf("%f\n",x1);
        return 0;

}

claws0n

798487012 发表于 2018-9-28 18:14
改成这样了，貌似可以了

晚一些才帮你改，应该还有问题

claws0n

798487012 发表于 2018-9-28 18:14
改成这样了，貌似可以了

#include <stdio.h>
#include <math.h>

double f(double x)
{        
    return 2*x*x*x+(-4)*x*x+3*x-6;
}

double f1(double x)
{        
    return 3*2*x*x+2*(-4)*x+3;
}

int main()
{
        double x,dx = 1;
        printf("输入一个估计值");
        scanf("%lf", &x);
        
        while (fabs(dx) >= 1e-7)
        {
                dx = f(x)/f1(x);
                x = x - dx;
        }
        
        printf("%f\n", x);
        return 0;

}

你不可以写 x1 = 6; 因为如果你输入 x = 6, fabs(x1-x) == 0，就不会进入循环了
没有加好友，呼叫是无效的 ^_^