Кому не сложно, помогите решить

Сильно модифицированный метод Эйлера (точность ~h⁴):
#include <cstdio>
#include <iostream>
double f(double x, double y)
{
return 2. * x + y * y;
}
int main()
{
double x, h = 0.1, hp = h * 0.5,
hq = h * 0.25, hphq = hp + hq,
y, u, v, w, Y[3];
int i;
Y[0] = 0.5;
for (i = 0; i < 10; i++)
{
x = i * h; y = Y[i];
u = y + hp * f(x + hq, y + hq * f(x, y));
v = u + hp * f(x + hphq, u + hq * f(x + hp, u));
w = y + h * f(x + hp, y + hp * f(x, y));
Y[i + 1] = (4. * v - w) / 3.;
}
for (i = 0; i < 11; i++)
{
x = i * h;
printf("\033[36m%3.1f\033[35m%8.4f\n", x, Y[i]);
}
getchar();
}
Слегка модифицированный метод Эйлера (точность ~h²):
#include <cstdlib>
#include <iostream>
double f(double x, double y)
{
return 2. * x + y * y;
}
int main()
{
double x, h = 0.1, hp = h * 0.5, y, Y[11];
int i; Y[0] = 0.5;
for (i = 0; i < 10; i++)
{
x = i * h; y = Y[i];
Y[i + 1] = y + h * f(x + hp, y + hp * f(x, y));
}
for (i = 0; i < 11; i++)
{
x = i * h;
printf("\033[36m%3.1f\033[35m%8.4f\n", x, Y[i]);
}
getchar();
}
Сравнение результатов (слева -более точный результат):