program mynormal1_lf
	version 8.0
	args lnf mu sigma
	quietly replace `lnf' = ln(normalden($ML_y1,`mu',`sigma'))
	end