Evaluating a double integral using Euler's method and Richardson extrapolation