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

#include <stdio.h>
#include <math.h>
void main(void)
{double x,y,e,a,b,c;
 int k,s;
 printf(“Введите аргумент и точность: “);
 scanf(“%lf%lf”,&x,&e);
 a=x; y=0;
 b=x; s=1; c=1;
 k=1;
 while (fabs(a)>e)
 {
 y=y+a;
 k=k+1;
 b=b*x*x;
 c=c*(2*k-2)*(2*k-1);
 s=s*(-1);
 a=s*b/c;
 }
 printf(“Для x=%lf и e=%lf y=%lf”,x,e,y)
}

// x - аргумент, значение функции которого надо вычислить
// epsilon - точность вычисления суммы
double MySum(double x,double epsilon)
{
    double s=0;// Тут копим сумму
    double an; // Тут храним N-ый член ряда
    int n;
    n=1; // начальное значение n
    an=x*x*x/3.0; // начальное значение an
    while(fabs(an)>epsilon) 
// Суммировать будем пока член ряда an не станет достаточно маленьким по модулю
    {
          s+=an; // суммируем очередной член ряда
          n++; // переходим к следующему члену 
          an*=-x*x*(4.*n*n-2.*n-5)/(4.*n*n-1); // пересчитываем a(n) через a(n-1)
    }
    return s; // получившаяся сумма
}
//---------------------------------------------
int main()
{
     double x,y,eps;
     cout<<"Введите аргумент и точность: ");
     cin>>x>>eps;
     cout<< MySum(x,eps);
     return 0;
}

#include <iostream>
#include <cmath>
int main()
{
    double x, eps;
    std::cout << "Аргумент x = "; std::cin >> x;
    std::cout << "Точность eps = "; std::cin >> eps;
    double t = x * x, u = t * x / 3, s = u;
    int n = 1;
    while (fabs(u) > eps)
    {
        n++;
        u *= -(t * (4 * n * n - 8 * n + 3)) / (4 * n * n - 1);
        s += u;
    }
    std::cout << "s = " << s;
    return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{
    int n=1;
    double x, q, p, an=1., S=0., eps;
    cout << "x = "; cin >> x;
    cout << "eps = "; cin >> eps;    
    
    p=-x;
    q=x*x;
    
    while (fabs(an) > eps)
    {        
        p*=-q;        
        an=p/(4.*n*n-1.);        
        S+=an;
        n++;
    }
    
    cout << "S=" << S << "\n";
    
system("pause");
return 0;
}