cwelłer

#include <iostream>
#include <iomanip>
using namespace std;

float pochodna(float x,float y)
{
    return x+x*y+y+1;
}
float yn1(float yn,float h,float xn)
{
    return yn+h*pochodna(xn,yn);
}
int main()
{
    float h=0.001;
    float b=0;
    float yn=1;
    float xn=-1;
    float k=(b-xn)/h;
    float x,y;

    for(int i=0;i<4;i++)
    {
        y=yn1(yn,h,xn);
        x=xn+h;

        yn=y;
        xn=x;
    }
    cout<<setw(10)<<xn;
    cout<<setw(10)<<yn<<endl;
    yn=1;
    xn=-1;
     for(int i = 0; i < k; i++)
    {

        y = yn + (h/2)*(pochodna(xn,yn) + pochodna(xn+h, yn+h*pochodna(xn,yn)));
        x = xn + h;
        yn=y;
        xn=x;
    }
    cout<<setw(10)<<xn;
    cout<<setw(10)<<yn<<endl;
    //cout<<"x"<<k<<"= "<<xn<<"  y"<<k<<"= "<<yn<<endl;
    return 0;
}