Create plot.py

Python script which compiles and runs the algorithm and plots against the closed-form solution for some values of n.
It also extracts the relative error and plots log10(rel_err) vs log10(h).