Вычислить значение функции заданной разложением в ряд Тейлора

#include <iostream>
#include <string>
#include <math.h>
 
using namespace std;
 
double cos_series_sum(double x,int n){
    //1-x^2/2!+x^4/3!+...
    double res = 1;
    double sign = 1, fact = 1, power = 1;
    
    for (int i = 1; i < n; i++){
        sign = sign * -1;//знак
        fact = fact * (2 * i - 1)*(2 * i); //факториал
        power = power * x * x; //возведение x в степень
        res = res + sign * power / fact; //значение
    } 
 
    return res;
}
 
double calculate_series(double x,int n){
    //sum((-1)^n*x^(2n)/(2n)!)
    double res = 1;//значение при n=0
    double sign = 1, fact = 1, power = 1;
    int k = 1;
    double x1 = 0;
 
    for (int i = 1; i < n; i++) {
        sign = pow((-1),i);//знак
        fact = fact * (2 * i - 1)*(2 * i); //факториал
        x1=pow(x, 2*i); //возведение x в степень
        res = res+sign*x1/fact; //значение
    }
 
    return res;
}
 
 
void main(){
    setlocale(LC_ALL, "Russian");
    double x = 0;
    double PI = 3.142;
    int n = 0;
    cout << "Введите x" << endl;
    cin >> x;
 
    cout << "Введите максимальное n" << endl;
    cin >> n;
 
    x = x * (PI / 180.0); //x переводим в радианы
    cout.precision(10); //количество знаков
    
    cout << "cos(x)=" << cos(x) << endl;
    cout << "sum((-1)^n*x^(2n)/(2n)!)=" << calculate_series(x, n) << endl;
    cout << "1-x^2/2!+x^4/3!+...=" << cos_series_sum(x, n) << endl;
    system("pause");
}

double mycos(double x,double eps)
{
    double s,an;
    int n;
    n=0;
    an=1;
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n); 
    }
    return s;
}