model {
  A[1,1] <- -exp(logK1)
  A[1,2] <-  exp(logK2)
  A[2,1] <-  exp(logK1)
  A[2,2] <- -exp(logK2)-exp(logK3)

  tau <- 1/(sigma*sigma)
  x0[1] <- 0
  x0[2] <- 100
  for (i in 1:noObs){
    M[i] ~ dnorm (pred[i], tau)
    pred[i] <- 100-sum(mexp(A*time[i])%*%x0)
  }
  sigma ~ dunif (0, 1000)
  logK1 ~ dunif(-10,10)
  logK2 ~ dunif(-10,10)
  logK3 ~ dunif(-10,10)
}