function f=density(theta)

global tau k y_low

f = 2*k*(1-tau)^(-k)*theta^(-2*k-1)*y_low^k;