n = 1;
digits=100;
xh = NSolve[HermiteH[n,x]==0,x,digits];
xvec = x /. xh;
wvec = {};
s[t_] := HermiteH[n,t];
Do [{xx = xvec[[i]],w=Sqrt[Pi] * Gamma[n] * 2^n / s'[xx] / HermiteH[n-1,xx],wvec = Append[wvec,w]},{i,1,n,1}];
Do [ 
      Print["    x(",i,") = ",N[xvec[[i]],30] ],{i,1,n,1}];
Do [   
      Print["    w(",i,") = ",FortranForm[N[wvec[[i]],30]] ],{i,1,n,1}];

n = 3;
digits=100;
xh = NSolve[HermiteH[n,x]==0,x,digits];
xvec = x /. xh;
wvec = {};
s[t_] := HermiteH[n,t];
Do [{xx = xvec[[i]],w=Sqrt[Pi] * Gamma[n] * 2^n / s'[xx] / HermiteH[n-1,xx],wvec = Append[wvec,w]},{i,1,n,1}];
Do [ 
      Print["    x(",i,") = ",N[xvec[[i]],30] ],{i,1,n,1}];
Do [   
      Print["    w(",i,") = ",FortranForm[N[wvec[[i]],30]] ],{i,1,n,1}];

n = 7;
digits=100;
xh = NSolve[HermiteH[n,x]==0,x,digits];
xvec = x /. xh;
wvec = {};
s[t_] := HermiteH[n,t];
Do [{xx = xvec[[i]],w=Sqrt[Pi] * Gamma[n] * 2^n / s'[xx] / HermiteH[n-1,xx],wvec = Append[wvec,w]},{i,1,n,1}];
Do [ 
      Print["    x(",i,") = ",N[xvec[[i]],30] ],{i,1,n,1}];
Do [   
      Print["    w(",i,") = ",FortranForm[N[wvec[[i]],30]] ],{i,1,n,1}];

n = 15;
digits=100;
xh = NSolve[HermiteH[n,x]==0,x,digits];
xvec = x /. xh;
wvec = {};
s[t_] := HermiteH[n,t];
Do [{xx = xvec[[i]],w=Sqrt[Pi] * Gamma[n] * 2^n / s'[xx] / HermiteH[n-1,xx],wvec = Append[wvec,w]},{i,1,n,1}];
Do [ 
      Print["    x(",i,") = ",N[xvec[[i]],30] ],{i,1,n,1}];
Do [   
      Print["    w(",i,") = ",FortranForm[N[wvec[[i]],30]] ],{i,1,n,1}];

n = 31;
digits=100;
xh = NSolve[HermiteH[n,x]==0,x,digits];
xvec = x /. xh;
wvec = {};
s[t_] := HermiteH[n,t];
Do [{xx = xvec[[i]],w=Sqrt[Pi] * Gamma[n] * 2^n / s'[xx] / HermiteH[n-1,xx],wvec = Append[wvec,w]},{i,1,n,1}];
Do [ 
      Print["    x(",i,") = ",N[xvec[[i]],30] ],{i,1,n,1}];
Do [   
      Print["    w(",i,") = ",FortranForm[N[wvec[[i]],30]] ],{i,1,n,1}];

n = 63;
digits=100;
xh = NSolve[HermiteH[n,x]==0,x,digits];
xvec = x /. xh;
wvec = {};
s[t_] := HermiteH[n,t];
Do [{xx = xvec[[i]],w=Sqrt[Pi] * Gamma[n] * 2^n / s'[xx] / HermiteH[n-1,xx],wvec = Append[wvec,w]},{i,1,n,1}];
Do [ 
      Print["    x(",i,") = ",N[xvec[[i]],30] ],{i,1,n,1}];
Do [   
      Print["    w(",i,") = ",FortranForm[N[wvec[[i]],30]] ],{i,1,n,1}];

n = 127;
digits=100;
xh = NSolve[HermiteH[n,x]==0,x,digits];
xvec = x /. xh;
wvec = {};
s[t_] := HermiteH[n,t];
Do [{xx = xvec[[i]],w=Sqrt[Pi] * Gamma[n] * 2^n / s'[xx] / HermiteH[n-1,xx],wvec = Append[wvec,w]},{i,1,n,1}];
Do [ 
      Print["    x(",i,") = ",N[xvec[[i]],30] ],{i,1,n,1}];
Do [   
      Print["    w(",i,") = ",FortranForm[N[wvec[[i]],30]] ],{i,1,n,1}];