关于用弦截法求根的问题

firewater

1. /*第二个代码是我自己写的，是正确的，能得到答案，第一个贴的是书上的答案，但是得不到正确答案，代码都是一样的，只是几个函数定义的位置不同，我的定义在后面，而答案的定义在前面，为什么得不到结果</span><span style="line-height: 1.5;">*/</span>

复制代码

1. #include <stdio.h>
   
2. #include <math.h>
   
3. float f(float x)               /*定义f函数，以实现f（x）=x*x*x-5*x*x+16*x-80*/
   
4. {
   
5.         float y;
   
6.         y = ((x - 5.0)*x + 16.0)*x - 80.0;
   
7.         return (y);
   
8. }
   
9. 
   
10. float xpoint(float x1, float x2)        /*定义xpoint函数，求出弦与x轴的交点*/
    
11. {
    
12.         float y;
    
13.         y = (x1*f(x2) - x2*f(x1) / (f(x2) - f(x1)));
    
14.         return (y);
    
15. }
    
16. 
    
17. float root(float x1, float x2)          /*定义root函数，求近似根*/
    
18. {
    
19.         float x, y, y1;
    
20.         y1 = f(x1);
    
21.         do
    
22.         {
    
23.                 x = xpoint(x1,x2);
    
24.                 y = f(x);
    
25.                 if (y*y1 > 0)
    
26.                 {
    
27.                         y1 = y;
    
28.                         x1 = x;
    
29.                 }
    
30.                 else
    
31.                 {
    
32.                         x2 = x;
    
33.                 }
    
34.         } while (fabs(y) >= 0.0001);
    
35.                 return (x);
    
36. }
    
37. 
    
38. void main()
    
39. {
    
40.         float x1, x2, f1, f2, x;
    
41.         do
    
42.         {
    
43.                 printf("input x1, x2:\n");
    
44.                 scanf("%f%f", &x1, &x2);
    
45.                 f1 = f(x1);
    
46.                 f2 = f(x2);
    
47.         } while (f1 * f2 >= 0);
    
48.         x = root(x1, x2);
    
49.         printf("A root of equation is%8.4f\n", x);
    
50. }

复制代码

以下代码是我自己写的：

1. #include "stdio.h"
   
2. #include "math.h"
   
3. 
   
4. int main(void)
   
5. {
   
6.         float x1, x2, f1, f2, x;
   
7.         float f(float x);                  /*求函数值*/
   
8.         float xpoint(x1, x2);             /*求中间变量X*/
   
9.         float root(x1, x2);              /*求根*/
   
10.         do
    
11.         {
    
12.                 printf("Please input two numbers :\n");
    
13.                 scanf("%f%f", &x1, &x2);
    
14.                 f1 = f(x1);
    
15.                 f2 = f(x2);
    
16.         } while (f1*f2 >= 0);
    
17.         x = root(x1, x2);
    
18.         printf("%f", x);
    
19.         return 0;
    
20. }
    
21. float f(float x)
    
22. {
    
23.         float y;
    
24.         y = x*x*x - 5.0 * x*x + 16.0 * x - 80.0;
    
25.         return (y);
    
26. }
    
27. float xpoint(float x1, float x2)
    
28. {
    
29.         float y;
    
30.         y = (x1*f(x2) - x2*f(x1)) / (f(x2) - f(x1));
    
31.         return (y);
    
32. }
    
33. float root(float x1, float x2)
    
34. {
    
35.         float x, y, y1;
    
36.         y1 = f(x1);
    
37.         do
    
38.         {
    
39.                 x = xpoint(x1, x2);
    
40.                 y = f(x);
    
41.                 if (y*y1>0)
    
42.                 {
    
43.                         x1 = x;
    
44.                         y1 = y;
    
45.                 }
    
46.                 else
    
47.                 {
    
48.                         x2 = x;
    
49.                 }
    
50.         } while (fabs(y) >= 0.0001);
    
51.         return (x);
    
52. }
    

复制代码

firewater

没人回吗？？？

雪是梅之香

如图所示，那句代码打错了，括号加错了地方应该是：

1. y = (x1*f(x2) - x2*f(x1)) / (f(x2) - f(x1));

复制代码

这样就可以

firewater

雪是梅之香 发表于 2015-1-20 19:04
如图所示，那句代码打错了，括号加错了地方应该是：
这样就可以

没有，括号没有打错。公式是这个，正确答案是5.000000

雪是梅之香

firewater 发表于 2015-1-20 21:25
没有，括号没有打错。公式是这个，正确答案是5.000000

但是改了以后答案就是5.000000
而且如实没有打错的话那就是和你的代码就不一样了

firewater

雪是梅之香 发表于 2015-1-20 21:28
但是改了以后答案就是5.000000
而且如实没有打错的话那就是和你的代码就不一样了

好吧，确实是我打错了，多谢你这么仔细的看啊，这个错误真是太不应该了。我被这个括号弄晕了.......谢谢了