Curve fitting in Stata