[julia-users] Re: Comparison of FE codes; Julia versus commercial FE solver

[Steven G. Johnson]

In particular, it looks like you just need:

param_coords, weights = Base.QuadGK.gauss(Float64, order)

in your GaussRule function.