This file conclude of five codes , four of them in Mathematica program and one in C++.

1.　　Erk4.nb. this code represent the explicit Runge Kutta method of order four for solving first order ODE.

2.　　RK45.nb This code represent the Embedded Runge Kutta Felhberg method for solving first order ODE.

3.　　EularStability.nb. This code represent how to find the stability region for the Eular method in 2D and 3D. its allow to modify the code to find the stability for any method.

4.　　RadauI.c . This code represent the fifth order implicit RadauI for solving first order ODE.

5.　　HEUN.nb. This code represent the third order Runge kutta method (Hune) Method for systems of ODE.

-This file conclude of five codes , four of them in Mathematica program and one in C++.

1.　　Erk4.nb.　this code represent the explicit Runge Kutta method of order four for solving first order ODE.

2.　　RK45.nb　This code represent the Embedded Runge Kutta Felhberg　method for solving first order ODE.

3.　　EularStability.nb. This code represent how to find the stability region for the Eular method in 2D and 3D.　its allow to modify the code to find the stability for any method.

4.　　RadauI.c . This code represent the fifth order implicit RadauI for solving first order ODE.

5.　　HEUN.nb. This code represent the third order Runge kutta method (Hune)　Method for systems of ODE.