sin(x)1.c

这是我大一的时候做的sin(x)源程序,库函数的实现

#include <stdio.h>
#include <math.h>
#define k 1000
#define pi 3.141592653589793
long double f(int n)
{
 int i;
 long double j;
 j=1;
 for(i=1;i<=n;i++)j=j*i;
 return(j);
}
long double sin1(float x)
{
 int n;
 long double j;
 long double y=0,p,q;
 for(n=1;n<k;n++)
 {
 j=f(n*2-1);
 p=pow((-1),(n-1));
 q=pow(x,(n*2-1));
 y=y+(p*q/j);
 }
 return(y);
}
void main()
{
 float x;
 long double y,temp;
 printf("Enter x:");
 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=sin1(x);
 printf("The result is:\nsinx:%lf\nsin1x:%lf\n",sin(x),y,(1-sin(x)/y)*100);
 getch();
}