Метод Рунге-Кутта 2 порядка на C++

дано задание реализовать данный метод, но возникла проблема с функцией. я приведу ниже код, он работает для диф уравнений первого порядка. но с новой функцией данные не совпадают с приведенными в таблице. Как возможно исправить или доработать программу?

#include "pch.h"
#include <iostream>
#include <cmath>
using namespace std;

double f(double x, double y) {
return 2*x - 3 + y;
//return 3 * sin(2 * y) + x;
//return -0.2*y - 10 * sin(x);
}
void enterData(double *a, double *b, double *h) {
cout << "Введите границы отрезка a и b: \n";
cin >> *a >> *b;
cout << "\nВведите шаг h:\n";
cin >> *h;
}

int main()
{
setlocale(LC_ALL, "rus");
double a, b, h;
int n;
double *x, *y, *k1, *k2, *k3, *k4;
enterData(&a, &b, &h);
//cout << endl << endl << "......." << a << "....." << b << "......" << h;

n = (a + b) / h;
//cout << endl << endl << "......." << a << "....." << b << "......" << 
h<<".... "<<n;
//system("pause");
x = new double[n+1];
y= new double[n+1];
k1= new double[n+1];
k2= new double[n+1];
k3= new double[n+1];
k4= new double[n+1];
cout << "\nВведите начальные условия.\nx = ";
cin >> x[0];
cout << "\ny = ";
cin >> y[0];
cout << endl;

//расчёт значений
for (int i = 1; i <= n; i++) {
    x[i] = x[0] + i * h;
    k1[i] = h * f(x[i - 1], y[i - 1]);
    k2[i] = h * f(x[i - 1] + h / 2, y[i - 1] + k1[i] / 2);
    k3[i] = h * f(x[i - 1] + h / 2, y[i - 1] + k2[i] / 2);
    k4[i] = h * f(x[i - 1] + h, y[i - 1] + k3[i]);
    y[i] = y[i - 1] + (k1[i] + 2 * k2[i] + 2 * k3[i] + k4[i])/6;
}
//печать результатов
for (int i = 0; i < n + 1; i++) {
    cout << "X[" << i << "] = " << x[i] << "    Y[" << i << "] = " << y[i] << 
";\n";
}

}