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)
}