sin(x)=x - x^3/3! + x^5/5! - ... + (-1^n)*x^(2n+1)/(2n+1)!

Yêu cầu bài toán: Tính sin(x) theo CT:
sin(x)=x - x^3/3! + x^5/5! - ... + (-1^n)*x^(2n+1)/(2n+1)! với độ chính xác 0.0001

Mã:

#include<stdio.h>
#include<conio.h>
#include<math.h>
double giaithua(double n)
{
	if(n==0||n==1)	return 1;
	return n*giaithua(n-1);
}
void main()
{
	clrscr();
	double tg,x,a,n;
	printf("Nhap goc can tinh (radian): ");
	scanf("%lf",&a);
	n=1;x=a;
	do
	 {
		tg=pow(a,2*n+1)/giaithua(2*n+1);
		x+=pow(-1,n)*tg;
		n++;
	 }
	while(fabs(tg)>0.0001);
	printf("Ket qua tinh theo CT:  %10.8lf",x);
	printf("\nKet qua tinh theo ham: %10.8lf",sin(a));
	getch();
}