SalmonLA

For future reference, it helps to specify an absolute tolerance to ensure convergence: hquadrature(sin, 0.0, 2*pi,abstol = 1e-8) (Found on https://stackoverflow.com/questions/29292614/how-to-do-two-variable-numeric-integration-in-julia#comment46814673_29302840)