关于用弦截法求根的问题
本帖最后由 firewater 于 2015-1-19 22:53 编辑/*第二个代码是我自己写的，是正确的，能得到答案，第一个贴的是书上的答案，但是得不到正确答案，代码都是一样的，只是几个函数定义的位置不同，我的定义在后面，而答案的定义在前面，为什么得不到结果</span><span style="line-height: 1.5;">*/</span>#include <stdio.h>
#include <math.h>
float f(float x) /*定义f函数，以实现f（x）=x*x*x-5*x*x+16*x-80*/
{
float y;
y = ((x - 5.0)*x + 16.0)*x - 80.0;
return (y);
}
float xpoint(float x1, float x2) /*定义xpoint函数，求出弦与x轴的交点*/
{
float y;
y = (x1*f(x2) - x2*f(x1) / (f(x2) - f(x1)));
return (y);
}
float root(float x1, float x2) /*定义root函数，求近似根*/
{
float x, y, y1;
y1 = f(x1);
do
{
x = xpoint(x1,x2);
y = f(x);
if (y*y1 > 0)
{
y1 = y;
x1 = x;
}
else
{
x2 = x;
}
} while (fabs(y) >= 0.0001);
return (x);
}
void main()
{
float x1, x2, f1, f2, x;
do
{
printf("input x1, x2:\n");
scanf("%f%f", &x1, &x2);
f1 = f(x1);
f2 = f(x2);
} while (f1 * f2 >= 0);
x = root(x1, x2);
printf("A root of equation is%8.4f\n", x);
}
以下代码是我自己写的：#include "stdio.h"
#include "math.h"
int main(void)
{
float x1, x2, f1, f2, x;
float f(float x); /*求函数值*/
float xpoint(x1, x2); /*求中间变量X*/
float root(x1, x2); /*求根*/
do
{
printf("Please input two numbers :\n");
scanf("%f%f", &x1, &x2);
f1 = f(x1);
f2 = f(x2);
} while (f1*f2 >= 0);
x = root(x1, x2);
printf("%f", x);
return 0;
}
float f(float x)
{
float y;
y = x*x*x - 5.0 * x*x + 16.0 * x - 80.0;
return (y);
}
float xpoint(float x1, float x2)
{
float y;
y = (x1*f(x2) - x2*f(x1)) / (f(x2) - f(x1));
return (y);
}
float root(float x1, float x2)
{
float x, y, y1;
y1 = f(x1);
do
{
x = xpoint(x1, x2);
y = f(x);
if (y*y1>0)
{
x1 = x;
y1 = y;
}
else
{
x2 = x;
}
} while (fabs(y) >= 0.0001);
return (x);
}
没人回吗？？？ 如图所示，那句代码打错了，括号加错了地方应该是：
y = (x1*f(x2) - x2*f(x1)) / (f(x2) - f(x1));这样就可以
雪是梅之香 发表于 2015-1-20 19:04
如图所示，那句代码打错了，括号加错了地方应该是：
这样就可以
没有，括号没有打错。公式是这个，正确答案是5.000000 firewater 发表于 2015-1-20 21:25
没有，括号没有打错。公式是这个，正确答案是5.000000
但是改了以后答案就是5.000000
而且如实没有打错的话那就是和你的代码就不一样了 雪是梅之香 发表于 2015-1-20 21:28
但是改了以后答案就是5.000000
而且如实没有打错的话那就是和你的代码就不一样了
好吧，确实是我打错了，多谢你这么仔细的看啊，这个错误真是太不应该了。我被这个括号弄晕了.......谢谢了
页:
[1]