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

798487012

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

1. #include<stdio.h>
   
2. #include<math.h>
   
3. 
   
4. double f(double x)
   
5. {       
   
6.         return 2*x*x*x+(-4)*x*x+3*x-6;
   
7. }
   
8. 
   
9. double f1(double x)
   
10. {       
    
11.         return 3*2*x*x+2*(-4)*x+3;
    
12. }
    
13. 
    
14. 
    
15. double main()
    
16. {
    
17. 
    
18.         double x=1,x1;
    
19.         x1=x-f(x)/f1(x);
    
20.        
    
21.         while (fabs(x1-x)>=1e-6)
    
22.         {
    
23.        
    
24.                 x=x1;
    
25.                 x1=x-f(x)/f1(x);
    
26.         }       
    
27.         return x1;
    
28.         printf("%f\n",x1);
    
29.         return 0;

复制代码

最佳答案

月排行榜 / 总排行榜

claws0n

2018-9-28 22:19:55

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

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. double f(double x)
   
5. {        
   
6.     return 2*x*x*x+(-4)*x*x+3*x-6;
   
7. }
   
8. 
   
9. double f1(double x)
   
10. {        
    
11.     return 3*2*x*x+2*(-4)*x+3;
    
12. }
    
13. 
    
14. int main()
    
15. {
    
16.         double x,dx = 1;
    
17.         printf("输入一个估计值");
    
18.         scanf("%lf", &x);
    
19.        
    
20.         while (fabs(dx) >= 1e-7)
    
21.         {
    
22.                 dx = f(x)/f1(x);
    
23.                 x = x - dx;
    
24.         }
    
25.        
    
26.         printf("%f\n", x);
    
27.         return 0;
    
28. 
    
29. }

复制代码

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

跳转到最佳答案楼层

798487012

@claws0n

798487012

改成这样了，貌似可以了

1. #include<stdio.h>
   
2. #include<math.h>
   
3. 
   
4. double f(double x)
   
5. {       
   
6.         return 2*x*x*x+(-4)*x*x+3*x-6;
   
7. }
   
8. 
   
9. double f1(double x)
   
10. {       
    
11.         return 3*2*x*x+2*(-4)*x+3;
    
12. }
    
13. 
    
14. 
    
15. double main()
    
16. {        double x,x1;
    
17.      printf("输入一个数");
    
18.         scanf("%lf",&x);
    
19.        
    
20.         x1=6;
    
21.        
    
22.         while (fabs(x1-x)>=1e-6)
    
23.         {
    
24.        
    
25.                 x=x1;
    
26.                 x1=x-f(x)/f1(x);
    
27.         }       
    
28. 
    
29.         printf("%f\n",x1);
    
30.         return 0;
    
31. 
    
32. }
    

复制代码

claws0n

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

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

claws0n

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

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. double f(double x)
   
5. {        
   
6.     return 2*x*x*x+(-4)*x*x+3*x-6;
   
7. }
   
8. 
   
9. double f1(double x)
   
10. {        
    
11.     return 3*2*x*x+2*(-4)*x+3;
    
12. }
    
13. 
    
14. int main()
    
15. {
    
16.         double x,dx = 1;
    
17.         printf("输入一个估计值");
    
18.         scanf("%lf", &x);
    
19.        
    
20.         while (fabs(dx) >= 1e-7)
    
21.         {
    
22.                 dx = f(x)/f1(x);
    
23.                 x = x - dx;
    
24.         }
    
25.        
    
26.         printf("%f\n", x);
    
27.         return 0;
    
28. 
    
29. }

复制代码

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