gethdot:

In musto101/wilcox_R: What the package does (short line)

Usage

gethdot(theta)

Arguments

theta

Examples

##---- 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 (theta) 
{
    p = 3
    ps = p * (p + 1)/2
    q = ps
    a = theta[1]
    b = theta[2]
    c = theta[3]
    vx = theta[4]
    vem = theta[5]
    vey = theta[6]
    sigmadot = matrix(NA, ps, q)
    sigmadot[1, ] = c(0, 0, 0, 1, 0, 0)
    sigmadot[2, ] = c(vx, 0, 0, a, 0, 0)
    sigmadot[3, ] = c(b * vx, a * vx, vx, a * b + c, 0, 0)
    sigmadot[4, ] = c(2 * a * vx, 0, 0, a^2, 1, 0)
    sigmadot[5, ] = c((2 * a * b + c) * vx, a^2 * vx + vem, a * 
        vx, a^2 * b + a * c, b, 0)
    sigmadot[6, ] = c((2 * b * c + 2 * a * b^2) * vx, (2 * c * 
        a + 2 * a^2 * b) * vx + 2 * b * vem, (2 * a * b + 2 * 
        c) * vx, c^2 + 2 * c * a * b + a^2 * b^2, b^2, 1)
    hdot = solve(sigmadot)
    return(hdot)
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.