---

#include <iostream.h>

#include <math.h>



double euler (double h, double y, double n)



int main (void)

    {

      double a, b, c;

      int n;



      cout<<"Value of h: "; //Step size

      cin>>a;

      cout<<"Initial value y(0): "; //Initial value

      cin>>b;

      cout<<"Value to be estimated y(x): "; //ex: If estimating y(1), write 1

      cin>>c;

      

      n=c/a; //Number of necessary iterations (int)



      printf("%3.8f\n", euler(a,c,n))



      system("pause") //Only there so Windows doesn't retardedly close the window immediately after execution

      return 0;

    }



double euler (double h, double y, double n) //Does the actual work

	{

      int i=0;

      double x=0;

      

      while (i <= c)

      	{

          y=y+h*(2*x*y) //Whatever's in parentheses is F(x,y), change at will

          x+=h;

          i++;

        }



      return y;

    }



Example: A step size of 0.2, an initial value of 1 and estimating y(1) returns 2.05058.

---