heqb.SGB <-
function(x,d,u,bound,shape1,index,...){
   h <- x[index]
   return(h)
}