Почему не заполняется массив vector?

#include <iostream>
#include <cmath>
#include <vector>
using namespace std;
 
class RunKutMer{
private:
    float a,y0,InEps,h;//x=a, f(x)=f(a)=y0, h - шаг
    float k1,k2,k3,k4,k5; //коэффициенты для рассчетов
    int i=0;
public:
vector <float> H; //создали массив для h
vector <float> Y; // создали массив значений для у
vector <float> X; // создали массив значений для x
 
vector <float> dY;
 
vector <float> Eps; //массив значений для точности эпсилон
 
 
    float f(float x, float y){
       return x*x+y*y-4;
      }
//нужно по-видимому запихивать получившиеся значения в массив
    //мы используем vector
 
    void EnterData(){
        cout<<"Enter a= ";cin>>a;
        X.push_back(a); //теперь в массиве X есть первый элемент
        cout<<"Enter the amount of function f(a)= ";cin>>y0;
        Y.push_back(y0);//теперь в массиве Y есть первый элемент
        cout<<"Enter Eps= ";cin>>InEps;//точность введенная пользователем, ОНА НЕ МЕНЯЕТСЯ
        //Eps.push_back(InEps);
        cout<<"Enter step h= ";cin>>h;
        H.push_back (h);  //записал в массив H первый элемент h
    }
    void OutPut(){
        cout<<"a= "<<a<<endl;
        cout<<"f(a)= "<<y0<<endl;
    }
    void Body(){
        i=0;
X[1]=2;X[2]=0;
    //   for(;;){//бесконечный цикл
//Почему не заполняется массивы X, Y, H и Eps???
           X[i+1]=X[i]+H[i];
 
           k1=(H[i]/3)*(f(X[i],Y[i]));
           k2=(H[i]/3)*(f(X[i]+H[i]/3,Y[i]+k1));
           //cout<<endl<<"k1= "<<k1<<endl<<"k2= "<<k2<<endl;
           k3=(H[i]/3)*(f(X[i]+H[i]/3,Y[i]+k1/2+k2/2));
           k4=(H[i]/3)*(f(X[i]+H[i]/2,Y[i]+(3.0/8)*k1+(9.0/8)*k3 ));
           k5=(H[i]/3)*(f(X[i]+H[i],Y[i]+(3.0/2)*k1-(9.0/2)*k3+6*k4 ));
 
           dY[i]=(1.0/2)*(k1+4*k4+k5);
           Y[i+1]=Y[i]+dY[i];
 
           i++;
           //сделали выше первую итерацию, теперь уточняем эпсилон и соответственно шаг h
           Eps[i]=(1.0/5)*(k1+4*k4-(9.0/2)*k3-(1.0/2)*k5);
 
           //условие проверки для эпсилон
           cout<<Eps[i];
      // }
    }
 
 
};
 
int main()
{
      RunKutMer W;
      W.EnterData();
      W.OutPut(); //так как переменные типа private, то обращаться к ним можно только через функции
      cout<<endl<<W.H[0]<<endl;cout<<W.X[0]<<endl;cout<<W.Y[0]<<endl;//проверил, что действительно занес все в vectora
      W.Body();
 
 
}

void Body(){
 
//X[1]=2;X[2]=0;
    //   for(;;){//бесконечный цикл
 
           X[i+1]=X[i]+H[i];
           X.push_back(X[i+1]);           // записали в массив Х
 
           k1=(H[i]/3)*(f(X[i],Y[i]));
           k2=(H[i]/3)*(f(X[i]+H[i]/3,Y[i]+k1));
 
           //cout<<endl<<"k1= "<<k1<<endl<<"k2= "<<k2<<endl;
 
           k3=(H[i]/3)*(f(X[i]+H[i]/3,Y[i]+k1/2+k2/2));
           k4=(H[i]/3)*(f(X[i]+H[i]/2,Y[i]+(3.0/8)*k1+(9.0/8)*k3 ));
           k5=(H[i]/3)*(f(X[i]+H[i],Y[i]+(3.0/2)*k1-(9.0/2)*k3+6*k4 ));
 
           dY[i]=(1.0/2)*(k1+4*k4+k5);
           dY.push_back(dY[i]);            // записали значение в массив dY
           Y[i+1]=Y[i]+dY[i];
                      Y.push_back(Y[i+1]); // записали значение в массив Y
 
           i++;
           //сделали выше первую итерацию, теперь уточняем эпсилон и соответственно шаг h
           Eps[i]=(1.0/5)*(k1+4*k4-(9.0/2)*k3-(1.0/2)*k5);
                      Eps.push_back(Eps[i]);
 
           //условие проверки для эпсилон
 
      // }
    }