h = 0.0001;
x = x1;
i = 0;
while(x)
{
        y[i]=y0;
        y[i]+=y1(t)*h;
        y[i]+=y2(t)*h*h/2.0;
        y[i]+=y3(t)*h*h*h/6.0;
        y[i]+=y4(t)*h*h*h*h/24.0;
        y0=y[i]; x+=h; i++;
}

// Написано под "Turbo C++ v3.0", запускать в эмуляции MS-DOS
#include <bios.h>
#include <bios.h>
#include <stdio.h>

// Вывод точки
#define SetPixel(X,Y,C) *((unsigned char far *)(0xA0000000L+0x0140*(Y)+(X)))=((unsigned char)(C))

// Шаг интегрирования
double h=0.0001;

// Наше диф. уравнение
double f(double X, double Y)
{
        double F;
        F=2.70*Y*Y*sin(M_PI*(0.50+8.00*X))/(X+0.10)-X*Y;
        return(F);
}

// факториал
double NF(unsigned int N)
{
        unsigned long int P=1.00;
        if(!N) return P;
        while(N) { P*=N; N--; }
        return (double)P;
}

// коэффициенты бинома Ньютона
double Cnm(unsigned int N, unsigned int M) { return NF(N)/(NF(M)*NF(N-M)); }

// дополнительная функция
double DF(double X, double *Y, unsigned int N)
{
        unsigned int I;
        double       S=0.00,Sign=1.00;
        for(I=0x00;I<=N;I++)
        {
                S+=Sign*Cnm(N,I)*f(X-h*(double)I,*(Y+0x03-I));
                Sign=-Sign;
        }
        return S;
}

// ну, если вы не знаете, что такое main(), то зачем это читаете?
void main()
{
        double X1=0.00,X2=1.00,Y0=1.00;
        double X,DY;
        double Y[0x04];
        unsigned int x,y,i;
        double MaxX,MaxY;

        // коэффициенты, см. Конечно разносные методы
        double G[0x04]={ 1.00, 1.00/2.00, 5.00/12.00, 3.00/8.00    };
        double g[0x04]={ 1.00,-1.00/2.00,-1.00/12.00,-1.00/24.00   };

        // точка максимального значения
        MaxX=X1; MaxY=Y0;

        // Включение графического режима 320*200*256
        asm MOV AX,0x0013;
        asm INT 0x10;

        // Подставляем начальные условия и производим расчет
        // первых 4-ёх значений по неявной формуле Адамса 2-го порядка точности
        X=X1;
        Y[0x00]=Y0;
        Y[0x01]=Y[0x00]+h*(f(X,Y[0x00])+f(X+h,Y[0x00]+h*f(X,Y[0x00])))/2.00; X+=h;
        Y[0x02]=Y[0x01]+h*(f(X,Y[0x01])+f(X+h,Y[0x01]+h*f(X,Y[0x01])))/2.00; X+=h;
        Y[0x03]=Y[0x02]+h*(f(X,Y[0x02])+f(X+h,Y[0x02]+h*f(X,Y[0x02])))/2.00; X+=h;

        // основной расчетный цикл
        while(X < X2)
        {
                // рассчитываем по явной формуле
                DY=0.00;  for(i=0x00;i<0x04;i++) DY+=G[i]*DF(X,Y,i);  DY*=h;
                // перемещаемся на шаг вперед
                for(i=0x00;i<0x03;i++) Y[i]=Y[i+0x01];  Y[0x03]+=DY;  X+=h;
                // две итерации на основе неявной формулы
                DY=0.00;  for(i=0x00;i<0x04;i++) DY+=g[i]*DF(X,Y,i);   DY*=h;
                Y[0x03]=Y[0x02]+DY;
                DY=0.00;  for(i=0x00;i<0x04;i++) DY+=g[i]*DF(X,Y,i);   DY*=h;
                Y[0x03]=Y[0x02]+DY;

                // проверяем максимальное значение
                if(MaxY < Y[0x03]) { MaxY=Y[0x03]; MaxX=X; }

                // вывод на экран
                x=0x0140*(X/(X2-X1));
                if(x<0x140)
                {
                        DY=f(X,Y[0x03]); y=0xC8*(0.50-0.0014*DY); if(y<0xC8) SetPixel(x,y,0x0D);
                        y=0xC8*(1.00-0.1*Y[0x03]); if(y<0xC8) SetPixel(x,y,0x0E);
                }
        }

   // пауза
   bioskey(0x00);

   // выкл. графического режима
   asm MOV AX,0x0003;
   asm INT 0x10;

   // Вывод данных
   printf("Max - %.14g,%.14gn",MaxX,MaxY);

   // пауза
   bioskey(0x00);
  }