Вывести на экран таблицу значений функции

#include <cmath>
#include <iomanip>
#include <iostream>
using namespace std;
 
double PI = 2*asin(1.0);
 
double Y(double x)
{
    return (3*pow(x, 2.0) - 6*PI*x + 2*pow(PI, 2.0)) / 12;
}
 
double S(double x, double e)
{
    double Sum = 0;
    double  k  = 1;
    double an  = cos(k*x) / pow(k, 2.0);
    for(k = k + 1; e < fabs(an); k = k + 1)
    {
        Sum += an;
        an  = cos(k*x) / pow(k, 2.0);
    }
    return Sum;
}
 
int main()
{
    double a;
    double b;
    double h;
    double e;
    double y;
    double s;
    cout<<"Enter a : ";cin>>a;
    cout<<"Enter b : ";cin>>b;
    cout<<"Enter h : ";cin>>h;
    cout<<"Enter e : ";cin>>e;
    cout<<"---------------------------"<<endl;
    cout<<"|   x   |  Y(x)  |  S(x)  |"<<endl;
    cout<<"---------------------------"<<endl;
    for(double x = a; x <= b; x = x + h)
    {
        y = Y(x);
        s = S(x, e);
        cout<<"|"<<setw(7)<<(x < 0 ? setprecision(4) : setprecision(5))<<x    
            <<"|"<<setw(8)<<(y < 0 ? setprecision(4) : setprecision(5))<<y  
            <<"|"<<setw(8)<<(s < 0 ? setprecision(4) : setprecision(5))<<s
            <<"|"<<endl;
    }
    return 0;
}