We apply the Midpoint method and the Modified Euler (Heun's) method to find approximate solutions a first-order initial value problem. We use Mathematica to generate the approximate solution. After finding the approximate solution, we find a curve, in this video, the Lagrange polynomial through all the points and check how well the Lagrange polynomial satisfies the differential equation.