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

#include <iostream>
#include <cmath>
using namespace std;
 
double f(double d)
{
    double Rez;
    
    Rez = (exp(d) - exp(-d))/2;
    return Rez;
}
 
int main()
{
    double x, s, an;
    int n, i;
    
    x = 0.2;
    an = x;
    for(i = 1; i < n; i++)
    {
        s += an;
        an *= pow(x,2)/(2*n*(2*n + 1));
    }
    cout << "s   = " << s << '\n';
    cout << "Rez = " << f(x) << '\n';
    
    return 0;
}

#include <iostream>
double fexp(double x)
{
    const double eps = 1.0e-4; // Точность !
    int n = 0;
    double u = 1.0, res = 1.0;
    while (u > eps || u < -eps)
    {
        n++;
        u *= x / n;
        res += u;
    }
    return res;
}
 
double Y(double x)
{
    return (fexp(x) - fexp(-x)) / 2;
}
 
double S(double x, int n)
{
    const double eps = 1.0e-5; // Точность !
    int i = 0;
    double u = x, res = x;
    while (i <= n)
    {
        i++;
        u *= (x * x) / (2 * i * (2 * i + 1));
        res += u;
    }
    return res;
}
 
int main()
{
    int n;
    std::cout << "n = ";    std::cin >> n;
    double x;
    do  // x = 0.1 ... 1
    {
        std::cout << "x = ";
        std::cin >> x;
    } while (x < 0.1 || x>1.0);
    std::cout << "Result: x = " << x << ";\tY = " << Y(x) << ";\tS = " << S(x, n);
    return 0;
}