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