Euler

#include<iostream>
using namespace std;
#define F(x,y) (2*x)/(3*y*y)
int main()
{
 
    float a,b,y,y0,h,x,S[3];
    int N;
    a=0;
    b=1;
    y0=1;        
    for(int i=0;i<3;++i)
    {
    
        y=y0;
        x=a;
        switch (i){
        case 0:
            h=0.01f;
            break;
        case 1:
            h=0.05f;
            break;
        case 2:
            h=0.1f;
            break;
        }
        N=(b-a)/h;
        for(int j=1;j<=N;++j)
        {
            y=y+h*F(x,y);
            x=x+h;
        }
        S[i]=y;
    }
    printf("The three result is :\n");
    printf("H1(0.01)   H2(0.05)   H3(0.1)\n");
    printf("%f   %f   %f\n",S[0],S[1],S[2]);
    return 0;
}