*! version 0.1.2  24sep2026  Joris Pinkse
*! ghpos d [q] [, options]: positive-weight cubature rule for the Gaussian weight; same as quadriceps gh
program ghpos, rclass
    version 16
    quadriceps gh `0'
    return add
end