	model
	{
		for( i in 1 : N ) 
		   {		    
				x[i] ~ dbs(alpha, beta)
		   }
		
	# Prior distributions of the model parameters	
			alpha ~ dunif(0.01, 5.0)
			beta~ dunif(0.01, 300.0)		
	}