model{
   for(i in 1:Nt){
      predy[i]<-M*exp(-c*days[i])+b
      recall[i]~dnorm(predy[i],lambda)
      predicts[i]~dnorm(predy[i],lambda)
   }
   # Prior
   lambda<-pow(sigma,-2)
   sigma~dunif(0,0.05)
   b~dunif(0,24)
   c~dunif(0,10)
}