cap prog drop myreg1

program define myreg1

	version 8

	args lnf theta1 theta2

	quietly replace `lnf'= ln(normden($ML_y1,`theta1',`theta2'))

end

ml model lf myreg1 (reg:y=x1 x2) /sigma

ml maximize