Question

C语言 sinx，cosx的近似值，求改错！急！

从键盘读入x值，按如下公式计算并输出sin x，cos x的值（精确度为0.000001，即按如下公式计算的值与实际sin x和cos
x的值相差不超过0.000001）。

(1) sin x = x-x3/3!+x5/5!-x7/7!+......

(2) cos x = 1- x2/2!+x4/4!-x6/6!+......
我写的
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int main()
{
int i=1,j=0,s=1,k,b;
double x,sinx=0.0,cosx=0.0;
printf("Please input x:");
scanf("%lf",&x);
for(;;i+=2,j++)
{
for(k=1;k<=i;k++)
{
s=s*k;
if(k==i)
{
b=s/k;
}
}
sinx=sinx+pow(-1,j)*pow(x,i)/s;
cosx=cosx+pow(-1,j)*pow(x,i-1)/b;
if(fabs(sinx-sin(x))<=0.000001&&fabs(cosx-cos(x))<=0.000001)
{
break;
}
}
printf("the sinx = %f \nthe cosx = %f",sinx,cosx);
return 0;
}
不知道为什么结果出不来

Answer 1

1、溢出问题，你这种方法，求sinx，当x比较大的时候，观察sinx的表达式第n项:x^(2n-1)/(2n-1)!，请注意，n!，当n>200，double型数据溢出;至于x，假设就等于10，显然double型数据范围为10^308，2n-1>308，也溢出…这是x=10的时候，但是x可以比10大的多，那么n可以变得很小也溢出
2、造成上述溢出原因主要是x太大…而sinx有周期性,sin(x)=sin(x-2pi)，那么x总可以减少到-派到+派之间…派是3.14，已经足够小，x^n/n!收敛速度较快，因为x^n较小，而n!增长很快，同时，n很小

Answer 2

int a=1;
int b=2;
int c=a/b;
你会发现c==0 懂了吧

Answer 3

s应该在for(；；i+2,j++)里面定义
没看出来表达式错了还是其他的。。。。