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

#include "pch.h"
#include <iostream>
using namespace std;
 
 
int fact(int i) // функция к заданию
{
    if (i < 0) return 0;
    if (i == 0) return 1;
    else return i * fact(i - 1);
 
}
 
int main()
{
    int step;
    double m, x, sign = 1, S = 0;
 
    cout << "m= "; cin >> m; cout << "x= "; cin >> x;
 
    for (int i = 1; i <= m; i++)
    {
        step = 2 * i - 1;
        S += sign * (pow(x, step) / (fact(2*i - 1)));
        sign *= -1;
    }
 
    cout << endl << "Sum = " << S << "\n";
 
    system("pause");
    return 0;
}

double mysin(double x,double eps)
{
    double s,an;
    int n;
    n=1;
    an=x;
    s=0;
    while(fabs(an)>eps)
    {
         s+=an;
         n++;
         an*=-x*x/(2.*n-1.0)/(2.0*n-2.0); 
    }
    return s;
}
int main()
{
    double x,eps;
 
    cout << "precision= "; cin >> eps; 
    cout << "x= "; cin >> x;
    cout << endl << "Sum = " << mysin(x,eps) <<' '<<sin(x)<< endl;
 
    system("pause");
    return 0;
}