## fixed a bug below: must be '%*% theta' instead of '* theta' (2010-03-15)
M <- rowSums((exp(-alpha * Wend) - exp(-alpha * Wstart)) %*% theta)
V <- exp(-alpha * W) * (1 - exp(-2 * alpha * (Tmax - W/2)))
## fixed a bug below: must be '%*% theta' instead of '* theta' (2010-03-15)
M <- rowSums((exp(-alpha * Wend) - exp(-alpha * Wstart)) %*% theta)
V <- exp(-alpha * W) * (1 - exp(-2 * alpha * (Tmax - W/2)))