model
   {
      for( i in 1 : N )
      {
      x[i] ~ dflex.weib(alpha, beta)
      }   
      # Prior distributions of the model parameters      
         
         alpha ~ dgamma(0.001, 0.001)
         beta ~ dgamma(0.001, 0.001)            
   }