Wyjaśnienie pisania funkcji sin i cos

Witam!
Mógłby mi ktoś wyjaśnić na czym polega zasada aproksymacji funkcji sinus i cosinus? Wzory widziałem- mam umieć wyliczyć przyblizona wartość funkcji sin, cos, sinh, cosh szeregiem taylora. Tylko... nie rozumiem co mam pisac w ogóle.

double sinus(double x, double eps)
{
	int i;
	double dd, s, licznik, mianownik;

	s = 0;
	i = 0;
	licznik = x;
	mianownik = 1.0;
	dd = x;

	while (dd >= eps)
	{
		i += 2;
		
	}
	sh += x;
	return s;
} 

rozumiem że while jest robiony względem precyzji eps. Jeżeli stosunek licznika do mianownika będzie mniejszy od żądanej dokładności to wylatujemy i mamy wynik. Ale co w środku napisać? Jak zrobić sinusa bez potęg i silni?

https://sourceware.org/git/?p=glibc.git;a=blob;f=sysdeps/ieee754/dbl-64/s_sin.c;hb=HEAD#l281

Bez łaski. Sam se zrobiłem....

double sinus(double x, double eps)
{
	x = x*M_PI / 180;
	int i;
	double dd, s, licznik, mianownik;

	s = x;
	i = 0;
	licznik = x;
	mianownik = 1.0;
	dd = x;

	while (abs(dd) >= eps)
	{
		i += 1;
		licznik = (-1)*licznik*(x*x);
		mianownik = mianownik* (2 * i) * (2 * i + 1);
		dd = licznik / mianownik;
		s += dd;
	}
	return s;
}