Firt cut work done:

Things still remaining:
1. Automate plotting with appropriate markers
2. For each of the 4 plots.
3. Think about the aberration near x = 0