add codes to compare calculated E0 with analytic soln