*! NJC 1.0.0 7 Dec 2006
program invgaussfit_lf
	version 8.1
	args lnf mu lambda
	qui replace `lnf' = ///
	0.5 * ln(`lambda') - ln(2 * _pi) - 3 * ln($S_MLy) ///
	- (`lambda'/(2 * `mu'^2 * $S_MLy) * ($S_MLy - `mu')^2)
end