Add analytic solution for comparison with Euler solution