Для каждого x, изменяющегося от a до b с шагом h, найти значения функции Y(x), суммы S(x) и |Y(x)–S(x)|

#include <math.h>
#include <iostream>
#include <conio.h>
 
using namespace std;
 
int main()
{
    double x, a, b, h, n, k, sum, r, y;
 
    cout << "Enter a,b,h,n:" << endl;
 
    cout << "a=", cin >> a;
    cout << "b=", cin >> b;
    cout << "h=", cin >> h;
    cout << "n=", cin >> n;
 
    for (x=a; x <= b; x += h)
    {
        for (k=0; k <= n; k++)
        {
            r = -x * x / (2 * k + 1) * (2 * k + 3);
            }
        y = sin(x);
        cout<<"sum=" << sum << '\t'<<"y=" << y << '\t'<<"y-sum=" << fabs(y - sum)<<'\n';
    }
 
}

#include <iostream>
#include <cmath> 
using namespace std;
 
int main()
{
    int n;
    double a, b, h, ak, S, Y;
 
    cout << "Enter a,b,h,n:\n";
 
    cout << "a=", cin >> a;
    cout << "b=", cin >> b;
    cout << "h=", cin >> h;
    cout << "n=", cin >> n;
    
    cout <<" x     Y(x)       S(x)       |Y(x)-S(x)|\n";
    for (double x=a; x<=b+h/2.; x+=h)
    {
        ak=x; S=x; Y=sin(x);
        for (int k=1; k<=n; k++)
        {
            ak*=-x*x/(2*k*(2*k+1));
            S+=ak;            
        }        
    cout << x << "   " << Y << "   " << S << "   " << fabs(Y-S) << "\n";
    }
    
system("pause");
return 0;
}