library(npmlreg)

g.func<-function(xv,nv){
  mv<-floor(nv/2)
  gv<-nv*choose(nv-1,mv)*xv^2
        *(pnorm(xv)^mv)*((1-pnorm(xv))^mv)
  return(gv)
}

n<-11
nk<-20
gh<-gqz(nk,minweight=1e-14)
x<-gh$location
w<-gh$weight
g<-g.func(x,n)
print(c(n,ik,sum(w*g)))