model

           {
     
           for (i in 1:1) {       # one subject
           c~dnorm(mu, 100)       # c:drug concentration
           mu<-(0.2+10*abs(PTD-7)/((PTD+10)*60))*D*1000*0.3/(4*bw*(0.3-k))*((1/(1-exp(-k*tau)))*exp(-k*ts)-(1/(1-exp(-0.3*tau)))*exp(-0.3*ts))    # oral concentration equation at steady state
           k<-cl/(4*bw)

           cl~dnorm(theta1, 0.106797)           

           theta1<-0.281*bw*(0.629*Dia+1-Dia)   # population estimate of Cyclosporine clearance
  
           }

           }