Вычисление приближенного значения tan(x) через ряд Тейлора/ряд Маклорена

#include <iostream>
#include <math.h>
#include <iomanip>
#include <omp.h>
#define  PI          3.14159265358979323846
double tgx(double); //функция, высчитывающая тангенс
long long  fact(int); //факториал
double vetBern(double); // числа Бернулли
 
int main() {
    int n;
    std::cin >> n;
    double start = omp_get_wtime();
    double Result = tgx((n * PI) / 180);
    double end = omp_get_wtime() - start;
    std::cout << std::fixed << std::setprecision(7) << " result = " << Result << std::endl;
    std::cout << "elapsed time = " << end << std::endl;
 
    return 0;
}
 
double tgx(double x)
{
    double z = 0.0;
    int i;
    if (fabs(x) < PI / 2) {
        for (i = 1; i < 14; i++)
        {
            //это ряд Тейлора z += vetBern(2 * i) *  ((pow(2, 2 * i)) * (pow(2, (2 * i)) - 1)  * pow(x, (2 * i) - 1)) / fact(2 * i);
            z += (vetBern(2 * i) * pow(-4, i) * (1 - pow(4, i)) * pow(x, ((2 * i) - 1))) / fact(2 * i);// ряд Маклорена
            std::cout << z << std::endl;
        }
    }
    else std::cout << " Wrong variable" << std::endl;
    return z;
}
 
long long fact(int  x) { // факториал
    if (x == 0 && x == 1)
        return 1;
    else {
 
        long long result = 1;
        for (int i = 2; i <= x; i++) {
            result *= i;
        }
        return result;
    }
}
 
 
double vetBern(double n)// расчет чисел Бернулли
{
    static_cast<int>(n);
    double bern = 1.0;
    double sum;
    if (n == 0)
    {
        sum = 1.0;
    }
    else
    {
        sum = 0.0;
        for (int k = 1; k <= n;k++)
        {
            sum += ((fact(k + 1) / (fact(n + 1)*fact(k - n))) * vetBern(n - k));
        }
        bern = ((-1) / (n + 1)) * sum;
    }
    return fabs(bern);
}

double vetBern(double n)// расчет чисел Бернулли
{
    static_cast<int>(n);
   
    double *ber=new int n;
     double bern = 1.0;
     ber[0]=1;
    double sum;
    if (n == 0)
    {
        sum = 1.0;
    }
    else
    {
        sum = 0.0;
        for (int k = 1; k < n;k++)
        {
            sum += ((-1) / (n + 1)) * ((fact(k + 1) / (fact(n + 1)*fact(k - n))) * ber[n - k]);
            ber[k]=sum;
            
        }
       // bern = sum;
    }
    delete []ber;
    return sum;
}

p = 1;
for(i = 1; i <= 14; ++i) {
    p *= 29 - i;
    p /= i;
}