798487012 发表于 2018-9-28 17:26:07

牛顿迭代法求根,代码运行错误,找不到原因求帮助呀!

用牛顿迭代法和二分法求根。方程为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;

798487012 发表于 2018-9-28 17:29:15

@claws0n

798487012 发表于 2018-9-28 18:14:06

改成这样了,貌似可以了{:10_279:}
#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 发表于 2018-9-28 20:00:03

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

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

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,就不会进入循环了
没有加好友,呼叫是无效的 ^_^
页: [1]
查看完整版本: 牛顿迭代法求根,代码运行错误,找不到原因求帮助呀!