标题:关于用弦截法求根的问题
只看楼主
firewater
Rank: 1
等 级:新手上路
帖 子:47
专家分:4
注 册:2014-7-13
结帖率:87.5%
 问题点数:0 回复次数:11 
关于用弦截法求根的问题
第一个代码是我自己写的,是正确的,能得到答案,第二个贴的是书上的答案,但是得不到正确答案,代码都是一样的,只是几个函数定义的位置不同,我的函数定义在后面,而答案的函数定义在前面,为什么得不到结果
程序代码:
#include "stdio.h"
#include "math.h"

int main(void)
{
        float x1, x2, f1, f2, x;
        float f(float x);                  /*求函数值*/
        float xpoint(float x1, float x2);             /*求中间变量X*/
        float root(float x1, float 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);
}





程序代码:
#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);
}


[ 本帖最后由 firewater 于 2015-1-20 16:16 编辑 ]
搜索更多相关主题的帖子: 正确答案 
2015-01-20 11:44
wmf2014
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:216
帖 子:2039
专家分:11273
注 册:2014-12-6
得分:0 
我不相信你第一个代码能得到正确答案:很明显的语法错误“float f(float x); ”,调用函数是不需要注明函数类型,同时提交的实参x也不需要注明类型,另外这几个函数都是有返回值的,你根本就不应用,怎么可能得到正确答案。反正我修改了你的语法错误后调试,在调用root函数后死循环了,根本执行不到scanf函数,第二段代码无需修改即编译通过,有运行结果。

能编个毛线衣吗?
2015-01-20 12:34
诸葛欧阳
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
来 自:流年
等 级:贵宾
威 望:82
帖 子:2790
专家分:14619
注 册:2014-10-16
得分:0 
第二个程序能运行啊,只是有个警告就是第十二行结果可能变为double型,赋给y时最好进行强制类型转换

一片落叶掉进了回忆的流年。
2015-01-20 12:49
诸葛欧阳
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
来 自:流年
等 级:贵宾
威 望:82
帖 子:2790
专家分:14619
注 册:2014-10-16
得分:0 
我到觉得第一个有点问题,你在主函数中申明函数格式不对

一片落叶掉进了回忆的流年。
2015-01-20 12:51
firewater
Rank: 1
等 级:新手上路
帖 子:47
专家分:4
注 册:2014-7-13
得分:0 
回复 2楼 wmf2014
我用的是vs2013    第一个得到的答案是4.999999,这个是正确的     而第二个得到的是-1,这是错误的答案,
2015-01-20 16:09
firewater
Rank: 1
等 级:新手上路
帖 子:47
专家分:4
注 册:2014-7-13
得分:0 
回复 2楼 wmf2014
可能我写的那几个注释让你误解了,那三个是对函数的声明,并不是调用函数
2015-01-20 16:11
firewater
Rank: 1
等 级:新手上路
帖 子:47
专家分:4
注 册:2014-7-13
得分:0 
回复 4楼 诸葛欧阳
不好意思,没看清
2015-01-20 16:14
firewater
Rank: 1
等 级:新手上路
帖 子:47
专家分:4
注 册:2014-7-13
得分:0 
回复 2楼 wmf2014
函数声明是错了,改了之后,能的到正确答案,但是的二个代码依然得不到正确答案啊
2015-01-20 16:15
小军123
Rank: 1
等 级:新手上路
帖 子:10
专家分:2
注 册:2014-11-5
得分:0 
第二个出现了数值越界的错误,但不知道哪里出问题了?提示的是-1.#INF
2015-01-20 18:57
小军123
Rank: 1
等 级:新手上路
帖 子:10
专家分:2
注 册:2014-11-5
得分:0 
y = (x1*f(x2) - x2*f(x1) / (f(x2) - f(x1)));
括号放的不对
2015-01-20 18:59



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-441288-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.139259 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved