Make an ODESolver subclass for the iterated Midpoint method. Implement the numerical method (E.50)–(E.51) in a subclass of ODESolver. The code should reside in a separate file where the ODESolver...

Make an ODESolver subclass for the iterated Midpoint method.

Implement the numerical method (E.50)–(E.51) in a subclass of ODESolver. The code should reside in a separate file where the ODESolver class is imported. One can either fix N or introduce an ǫ and iterate until the change in |v q − vq−1| is less than ǫ. Allow the constructor to take both N and ǫ as arguments. Compute a new v q as long as q ≤ N

or |vq − vq−1| > ǫ. Let N = 20 and ǫ = 10−6 by default. Name of program file: MidpointIter.py.