Вычисление тангенса С++

[firerap]

Требуется вычислить значения тангенса, используя разложение в ряд (где-то в пределах от -9 до 9), но например при x=6 выпадает floating point invalid, я так понимаю частично проблема в факториале, но я думаю еще где-то накосячил. Подскажите, можно ли поправить код, что бы он корректно работал без длинной арифметики и как это можно реализовать?

PHP код:

#include <stdio.h>
#include <conio.h>
#include <math.h>

#include <windows>


int fact(int n) //Max n=12, overflow
{
if (n>12) puts("ERROR OVERFLOW");
    if(n==0) {return 1;}
       else return n*fact(n-1);
}

double DzetaRiman(int s)
{
 double t = 0.000001;
 double dzeta = 0;
 int i = 1;
 double mem = 1;
 while (mem>t) {
 mem = 1/ pow(i,s);
 i++;
 dzeta +=mem;
 }
 return dzeta;
}

double EvenBernulli(int n)
{ //Only for even n
int fsign=-1;
if((n/2+1)%2 == 0) fsign=1;
double r = 2*fsign*DzetaRiman(n)*fact(n)/pow(2*M_PI,n);
return  r;
}

double _tg(double x)
{
 double mem = 1;
 int n = 1;
 double rtg = 0;
 double t = 0.001;
 while (mem>t){
  mem = EvenBernulli(2*n)*pow(-4,n)*(1-pow(4,n))*pow(x,2*n-1)/fact(2*n);
  rtg+=mem;
  n++;
  }
  return rtg;
}


void main(void){

 printf("%f", _tg(6));
 system("PAUSE");

} 


[type_Oleg]

Вот тут написано, что формула для |x|< пи/2
http://ru.wikipedia.org/wiki/%D7%E8%...ED%F3%EB%EB%E8

Может в этом дело, надо приводить x к интервалу -пи/2 ..пи/2

[firerap]

Проверил при х = 1.5, не работает. Корректно находит тангенс только начиная с где-то х=0.8