Needs["NumericalDifferentialEquationAnalysis`"]
Module[{precision,nmax,d,GLQ},
precision=32;(*digits of precision*)
nmax=36;(*largest order*)
GLQ=Flatten@{nmax,
(d=GaussianQuadratureWeights[#,-1,1,precision];
{d,ConstantArray[{0.,0.},nmax-Length[d]]})
&/@Range[nmax]};

(*export*)
Export[StringJoin[
NotebookDirectory[],"gauss_legendre_quadrature_n",
ToString@nmax,"p",ToString@precision,".dat"],
GLQ];
]