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

#include <iostream>
#include <cmath>
 
using namespace std;
 
double my_ch(double x, double eps)
{
    double s,a,n;
    s=0;
    a=1;
    n=0;
    while(abs(a)>eps)
    {
        s+=a;
        a=a*x*x/((2*n+1)*(2*n+2));
        n++;
    }
    return s;
}
 
int main()
{
    double x=0.123;
    cout << my_ch(x,1.0e-7) << endl;
    cout << 0.5*(exp(x)+exp(-x)) << endl;
    return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{   
    int n=1, N;
    double x, eps, S=1., Y, an=1.;    
        
    cout << "eps="; cin >> eps;
    cout << "x="; cin >> x;
    cout << "N="; cin >> N;
   
    while(fabs(an)>eps) // пока модуль очередного слагаемого больше eps цикл работает   
    {     
        an*=x*x/(2.*n*(2.*n-1.)); // вычисление очередного слагаемого ряда
        S+=an; // суммирование
        n++; // счётчик слагаемых
        if(n==N) break; // если количество слагаемых равно N выход из цикла      
    }
    
    Y=cosh(x); // вычисление ch(x) стандартной функцией
    
    cout << "S=" << S << "  Y=" << Y << "\n";  
    
system("pause");
return 0;
}