y=sin(x)的源程序
请大家帮忙写一下y=sin(x)的源程序，不能用库函数注意精度问题，终止循环不是用#define N XXXX i<=N
谢谢~~~
离散公式是
sin x=x-x^3/3!+x^5/5!-...+(-1)^(n-1)*x^(2n-1)/(2n-1)!+...
（这里x是弧度）
2006-03-09 12:50
2006-03-09 13:08
2006-03-09 18:31
#define k 50
#define pi 3.1415926
double factorial(int z) /*阶乘，递归算法实现*/
{if(z==0||z==1)
return (1);
else
return (factorial(z-1)*z);
}
double mysin(double x)
{
int n;
long double v;
double y=0,o,p;
for(n=1;n<k;n++)
{
v=factorial(n*2-1);
o=pow((-1),(n-1));
p=pow(x,(n*2-1));
y=y+(o*p/v);
/*
printf("%f %f %f %f\n",v,o,p,y);
*/
}
return(y);
}
void main()
{
float x;
double y,temp;
printf("Enter the randian:");
scanf("%f",&x);
if (x>pi)
{
temp=-2*pi;
while(x>pi)
{
x+=temp;
}
}
if (x<-pi)
{
temp=2*pi;
while(x<-pi)
{
x+=temp;
}
}
y=mysin(x);
printf("sinx:%f\nmysinx:%f\ndisparity:%f%",sin(x),y,(1-sin(x)/y)*100);
getch();
}
怎样把#define k 50
去掉，然后设立另一个方法跳出循环~~~
2006-03-09 18:49