function f=density(theta)

global tau k y_low

%
f = 
