gauss_quad_test:

(nq,nn) =  1 2   err =  2.6737861956088516
(nq,nn) =  2 2   err =  4.952372474195119
(nq,nn) =  4 2   err =  15.846545858225449
(nq,nn) =  8 2   err =  4.872175110465303
(nq,nn) =  16 2   err =  0.5935253374754232

(nq,nn) =  1 4   err =  36.60358319688305
(nq,nn) =  2 4   err =  11.223385608114558
(nq,nn) =  4 4   err =  2.6854943428877256
(nq,nn) =  8 4   err =  0.022686884229045745
(nq,nn) =  16 4   err =  0.0034647229847344363

(nq,nn) =  1 8   err =  6.093409776251892
(nq,nn) =  2 8   err =  1.0719751914246594
(nq,nn) =  4 8   err =  0.03481386342239645
(nq,nn) =  8 8   err =  3.897080620163251e-05
(nq,nn) =  16 8   err =  1.0270895245412248e-11

(nq,nn) =  1 16   err =  0.04051313353795294
(nq,nn) =  2 16   err =  0.004867168343409389
(nq,nn) =  4 16   err =  0.0006374271004752075
(nq,nn) =  8 16   err =  2.087421790974986e-09
(nq,nn) =  16 16   err =  3.552713678800501e-15

(nq,nn) =  1 32   err =  0.013209186871865342
(nq,nn) =  2 32   err =  0.0005010046178419714
(nq,nn) =  4 32   err =  2.315887904558167e-06
(nq,nn) =  8 32   err =  7.815970093361102e-14
(nq,nn) =  16 32   err =  0.0