Циклические алгоритмы

[DashaShirshova]

Для каждого x, изменяющегося от a до b с шагом h, найти значения функции Y(x), суммы S(x) и |Y(x)–S(x)| и вывести в виде таблицы. Значения a, b, h и n вводятся с клавиатуры.
Работу программы проверить для a = 0,1; b = 1,0; h = 0,1; значение параметра n=10.

#include <math.h>
#include <stdio.h>
#include <iostream.h>

int main()
{
   double a,b,h,S=0,Y,x,p=1,r;
   int i,n;
   
	cout<<"Vvedite a:";
	cin>>a;
	cout<<"Vvedite b:";
	cin>>b;
	cout<<"Vvedite h:";
	cin>>h;
	cout<<"Vvedite n:";
	cin>>n;
   
    for (x=a; x<=b; x+=h)
	{
		Y=(exp(x)+exp(-x))/2;
		{
			r=S=1;
			for(i=1; i<=n; i++)
			{
                       r=r*(x*i)/(2*i);
                       S+=r;
             }

		}

		printf("\n x= %8.2lf   S(x)= %8.5lf    Y(x)= %8.5lf    |Y(x)-S(x)|= %8.5lf", x,S,Y,fabs(Y-S));	
		i++;
	}
return 0;
}

Помогите правильно организовать S(x)
Подсказка: S(x) и Y(x) равны