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 
  
       }

       }