hammset: function to do ...
In quaddisc: Quadratic discrepancies

Description Usage Arguments Details Value Note Author(s) References See Also Examples

~~ A concise (1-5 lines) description of what the function does. ~~

1	hammset(iN, vP = 0)

`iN`	~~Describe `iN` here~~
`vP`	~~Describe `vP` here~~

~~ If necessary, more details than the description above ~~

~Describe the value returned If it is a LIST, use

`comp1`	Description of 'comp1'
`comp2`	Description of 'comp2'

...

~~further notes~~

~Make other sections like Warning with Warning .... ~

~~who you are~~

~put references to the literature/web site here ~

~~objects to See Also as help, ~~~

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (iN, vP = 0) 
{
    if (vP == 0) {
        mX <- matrix(nrow = iN, ncol = 1)
        for (k in 1:iN) mX[k, 1] <- (2 * k - 1)/(2 * iN)
        return(mX)
    }
    else {
        iD <- length(vP) + 1
        mX <- matrix(nrow = iN, ncol = iD)
        for (k in 1:iN) {
            mX[k, 1] <- (2 * k - 1)/(2 * iN)
            for (s in 2:iD) {
                dY <- radinv(k, vP[s - 1])
                mX[k, s] <- dY
            }
        }
        return(mX)
    }
  }