add norm=... keyword argument to quadgk

As discussed in JuliaLang/ODE.jl#7, there is no reason not to let the user specify the desired error norm to use in the convergence criterion.   (e.g. there is potentially a big difference betwen converging in the L2 and Linf norms if you have a vector-valued integrand.)

cc: @tlycken