#include <stdio.h>
#include <math.h>
 
double f(double x, double y) {
 return y * (y - 1.0)/x;
} 
const int n = 50;
 
int main()
{
    int i;
    double x, y;
    double a, b, h, eps, k1, k2, k3, k4;
 
 a=1.0;b=2.0;h=(b-a)/n;
    x=a;y=0.5;
 
    for (i = 1; i <= n; i++)
    {
        k1=f(x,y);
        y=y+h*f(x+h/2.0, y+(h/20)*k1);
        x=x+h;
 
    }
    eps=1.0/3.0-y;
    printf("x = %.16lf\n", x);
 printf("eps = %.16lf\n", eps); 
    return 0;
}

I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+Cgpkb3VibGUgZihkb3VibGUgeCwgZG91YmxlIHkpIHsKIHJldHVybiB5ICogKHkgLSAxLjApL3g7Cn0gCmNvbnN0IGludCBuID0gNTA7CgppbnQgbWFpbigpCnsKICAgIGludCBpOwogICAgZG91YmxlIHgsIHk7CiAgICBkb3VibGUgYSwgYiwgaCwgZXBzLCBrMSwgazIsIGszLCBrNDsKCiBhPTEuMDtiPTIuMDtoPShiLWEpL247CiAgICB4PWE7eT0wLjU7CiAgICAKICAgIGZvciAoaSA9IDE7IGkgPD0gbjsgaSsrKQogICAgewogICAgICAgIGsxPWYoeCx5KTsKICAgICAgICB5PXkraCpmKHgraC8yLjAsIHkrKGgvMjApKmsxKTsKICAgICAgICB4PXgraDsKICAgICAgICAKICAgIH0KICAgIGVwcz0xLjAvMy4wLXk7CiAgICBwcmludGYoInggPSAlLjE2bGZcbiIsIHgpOwogcHJpbnRmKCJlcHMgPSAlLjE2bGZcbiIsIGVwcyk7IAogICAgcmV0dXJuIDA7Cn0=