标题:泰勒公式求sin(x)的近似值
只看楼主
嘻哈哈嘻嘻
Rank: 1
等 级:新手上路
帖 子:12
专家分:0
注 册:2017-3-1
结帖率:0
已结贴  问题点数:10 回复次数:2 
泰勒公式求sin(x)的近似值


代码:
程序代码:
#include <stdio.h>
#include <math.h>
int main()
{
    int n, i;
    double x, result, num, fenmu, fenzi;

    while(scanf("%lf %d",&x,&n)!=EOF)
    {
        i=1;
        fenzi=1;
        fenmu=result=x;
        num=pow(x,i)/fenzi;

        do
        {
            fenzi *= (i+1)*(i+2);
            fenmu *= -x*x;
            num = fenmu/fenzi;
            result += num;
            i += 2;
        }while(fabs(num)>=pow(10,-n));

        printf("%.*lf\n",n,result);
    }
    return 0;
}

结果在31后就开始出错了。求助大佬们。


[此贴子已经被作者于2017-3-28 12:49编辑过]

搜索更多相关主题的帖子: result 
2017-03-28 10:34
rjsp
Rank: 20Rank: 20Rank: 20Rank: 20Rank: 20
等 级:版主
威 望:507
帖 子:8890
专家分:53117
注 册:2011-1-18
得分:10 
不知道哪一步出错了
你都不知道哪一步出错了,别人怎么帮你?
最起码,你得知道是编译出错,还是运行结果出错。

若是运行结果不对,我猜是因为 double 存储 fenzi 和 fenmu 时,因为数值过大而失真。活该,滥定义没用的变量就该受惩罚
程序代码:
#include <stdio.h>
#include <math.h>

int main( void )
{
    double x; unsigned n;
    while( scanf("%lf%u",&x,&n) == 2 )
    {
        double result = x;

        double item = x;
        for( unsigned i=0; fabs(item)>=pow(10.0,0.0-n); i+=2 )
        {
            item *= -x*x/(i+2)/(i+3);
            result += item;
        }

        printf( "%.*lf\n", n, result );
    }

    return 0;
}

2017-03-28 12:38
嘻哈哈嘻嘻
Rank: 1
等 级:新手上路
帖 子:12
专家分:0
注 册:2017-3-1
得分:0 
回复 2楼 rjsp
了解了。是结果不对,我也猜测是fenzi和fenmu过大,但改不对。
2017-03-28 12:58



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




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

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