priormean <- 3
priorvar <- 4
S0 <- 2*priormean/priorvar
nu <- S0*priormean
x <- c(3,4,2,1,2,3)
n <- length(x)
T <- sum(x)
nuprime <- 2*T + nu
S1 <- 2*n + S0
postmean <- nuprime/S1
postvar <- 2*nuprime/S1^2
cat("Posterior has mean",postmean,"and variance",postvar,"\n")