d. Run the Matlab code to perform FOUR iterations of the Crank-Nicolson scheme for this problem. For each time-step considered, print out a table showing the exact solution, the numerical
solution, and the difference between these values for each mesh point. These should show at least 4 iterations of the numerical method.