Численные методы (лаба номер 4)

#include <iostream>
#include <cmath>
#include <iomanip>

#define _USE_MATH_DEFINES
using namespace std;
const double e = 2.718281828459045;
//фунkция для расчета производной
double func(double x,double y) {
    double yp = 0.218 * (x * x + sin(1.6 * x)) + 0.718*y;
    return yp;
}

//фунkция полченная аналитичесkи
double afunc(double x) {
    double f = -0.3036 * x * x - 0.8457 * x + 1.501 * pow(e,(0.718 * x)) - 0.0508 * sin(1.6 * x) - 0.1134 * cos(1.6 * x) - 1.1779;
    return f;
}

int main()
{
    //начальные условия
    double x = 0.2;
    double y = 0.25;

    cout << setw(8) << "lom otr";
    cout << "    ";
    cout << setw(8) << "analit";
    cout << "    ";
    cout << setw(8) << "difference" << endl;

    for (int i=0; i<10; i++) {
       // вычисление значения аналитичеkсого решения
        double f = afunc(x);
        // вывод
        cout<< setw(8)<<y;
        cout<<"    ";
        cout<<setw(8)<< f;
        cout << "    ";
        cout <<setprecision(4)<< abs(y-f) << endl;

        //вычисление значения y i+1(ломанные отрезkи)
        y = 0.1*func(x,y) + y;
        x += 0.1;
    }

}