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

#include <iostream.h>
#include <math.h>
#include <iomanip.h>
 typedef double (*uf)(double, double, int &);
void tabl(double, double, double, double, uf);
double y(double, double, int &);
double s(double, double, int &);
int main()
{
cout << setw(8) <<"x"<< setw(15) <<"y(x)"<< setw(10) << "k" << endl;
tabl(0.1,0.8,0.1,0.001,y);
cout << endl;
cout << setw(8) <<"x"<< setw(15) <<"s(x)"<< setw(10) << "k" <<endl ;
tabl(0.1,0.8,0.1,0.001,s);
return 0;
}
void tabl(double a, double b, double h, double eps, uf fun)
{
33
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 sin(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(x,2)/(2*k*(2*k+1));
 a *= -c;
 sum += a;
 k++;
 }
return sum;
}

#include <iostream>
#define M_PI 3.141596
using namespace std;
 
void Y(double x)
{
    cout << "Y(" << x << ") = " << M_PI * (M_PI - 2 * abs(x)) / 8 << '\n';
}
 
void S(double x, double e, int &iterCount)
{
    double summ = 0, lVal = 1;
    cout << "S(" << x << ") = ";
    for (iterCount = 1; abs(lVal) > e; iterCount ++)
    {
        lVal = (cos(2.0 * iterCount - 1) * x) / (pow(2.0 * iterCount - 1, 2.0));
        summ += lVal;
        cout << lVal;
        if (abs(lVal) > e)
            cout << " + ";
        else
            cout << " ";
    }
    cout << "= " << summ << '\n';
}
 
int main(int argc, char *argv[])
{
    setlocale(0, "Russian");
    double A = -2, B = 2, e = 0.01;
    int iter;
    for(double x = A; x <= B; x += 0.1)
    {
        for(int i = 0; i < 60; i++) cout << '-'; cout << '\n';
        cout << "x = " << x << '\n';
        Y(x);
        S(x, e, iter);
        cout << "Количество итераций: " << iter << '\n';
    }
    system("pause");
    return 0;
}

cos((2.0 * iterCount - 1) * x)

#define M_PI 3.14159

const double M_PI = 3.141596

 int
            a=-2,//min
            b=2;//max
    double
            e=0.0001,//e
            step=1;//step;
 
    cout<<"results for e= "<<e<<endl;
    for(double x=a,i=1,iter=1;x<=b;x+=step,i++,iter=1){
        double
                yabs=(M_PI*(M_PI-2*std::abs(x)))/8,//abs value
                real=(cos((2*iter-1)*x))/pow((2*iter-1),2); //real value
        cout<<"abs value x"<<i<<"="<<x<<"   y"<<i<<"="<<setw(8)<<setprecision(6)<<yabs<<endl;
        for(;abs(yabs-real)>e;iter++,real+=(cos((2*iter-1)*x))/pow((2*iter-1),2));
        cout<<"count iteration "<<iter<<" value y"<<i<<"="<<real<<endl<<endl;
    }