Счёт с первого элемента массива

#include<iostream>
#include<math.h>
#include<iomanip>
using namespace std;
int main()
{
    double a, b, h, x, y, s, p;
    int n, i = 1;
 
    cout << "input a,b,n" << endl;
    cin >> a >> b >> n;
    h = (b - a) / 10;     
    x = a;
    do
    {
        s= p = x/6;
        for (i = 2; i <= n; i++)
        {
            p *=x/(2*i*(2*i+1));    
            s +=i*i*p;
        }
        y = (1. / 4)*(((x + 1) / sqrt(x))*sinh(sqrt(x)) - cosh(sqrt(x)));       
        cout << setw(15) << x << setw(15) << y << setw(15) << s << endl;
        x += h;
    } while (x <= b + h / 2);
    cout << endl;
 
    return 0;
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{
    int n=120;
    double a=0.2, b=0.8, h, Y, S, an, eps=1e-9;     
    h = (b-a)/10; 
    
    cout <<"  x       Y         S\n";
    for (double x=a; x<=b+eps; x+=h) 
    {
       p=1.0;
       S=0.0;
       for (int i = 1; i <= n; i++)
       {
       p*=x/((2*i+1)*2*i); 
       S+=i*i*p;
       }
       Y = 1./4*((x+1)*sinh(sqrt(x))/sqrt(x)-cosh(sqrt(x))); 
       cout <<setw(4)<<x<<"  "<<Y<<"  "<<S<<endl;
    }     
system("pause");
return 0;
}