Вывести на экран таблицу значений функции Y(x) и ее разложения в ряд S(x) с точностью ε

#include <iostream>
#include <math.h>
#include <iomanip>
using namespace std;
 
typedef double (*uf)(double, double, int &);
 
void tabl(double, double, double, double, uf);
double y(double, double, int &);
double s(double, double, int &);
 
int factorial(int );
 
 
void tabl(double a, double b, double h, double eps, uf fun) {
    int k=0;
    double sum;
    for (double x=a; x<b+h/2; x+=h) {
        sum=fun(x,eps,k);
        cout << setw(8) << x << setw(15) << sum << setw(10) << k << endl;
    } }
 
//////////// обозначим функцию у(х)//////////////////////////////////////
double y(double x, double eps, int &k) {
    return sqrt(x+1)-1-x/2;
}
 
/////////// обозначим разложение ряда s(x)///////////////////////////////
 
 
 
double s(double x, double eps, int &k) {
    double a,c,sum; sum=a=c=x;
    k=1;
    while (fabs(c)>eps) {
        c = pow(-1,k-1)*factorial(factorial(2*k-3))*pow(x,k)/factorial(2*k); a *= -c;
        sum += a;
        k++;
    }
     return sum; }
 
int factorial(int &k)
{
    int res=1;
    int i;
    for (i=1; i<=k;i++)
    {
        res *= i;
    }
    return res;
}
 
//////////////// основная функция/////////////////////////////////////////
int main()
{
    cout << setw(8) <<"x"<< setw(15) <<"y(x)"<< setw(10) << "k" << endl; tabl(-0.9,0.9,0.1,0.001,y);
    cout << endl;
    cout << setw(8) <<"x"<< setw(15) <<"s(x)"<< setw(10) << "k" <<endl ; tabl(-0.9,0.9,0.1,0.001,s);
    return 0;
}

#define PRECISION 1E-9   
 
double S(double x, unsigned &iterations)
{
  unsigned k = 2;
  double sum = 0.0, m = -(x * x) * (2 * k - 3) / (2 * 4);
  iterations = 1;
  
  while (fabs(m) > PRECISION)
  {
    sum += m;
    k++;
    m *= -x * (1 - 3.0 / (2.0 * k));  
    iterations++;        
  }
  
  return sum;
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{    
    int n;
    double a=-0.9, b=0.9, h=0.1, an, eps=1e-3, S, Y, delta=1e-9; 
    
    cout <<  "  x     S(x)      Y(x)\n";
    for (double x=a; x<=b+delta; x+=h)
    {
        if(fabs(x)<delta) {x=0.0;S=0.0;Y=0.0;}
        else 
        {
        n=2;
        Y=sqrt(1.+x)-1.-x/2.;
        
        an=-x*x/8.;
        S=0.0;        
        while(fabs(an)>eps)    
        {           
           S+=an;           
           n++;           
           an*=-x*(2*n-3)/(2*n);            
        }
        }
    cout << setprecision(1) << fixed << x << setprecision(5) << "  " << S << "  " << Y << "\n";  
    }
    
system("pause");
return 0;
}