Add codes to compare our numerical results with analytic solutions