twoway function y=normal(x), range(-3 3)
— visualizing comparative features of different explicit functions
— visualizing theoretical vs. empirical results (e.g. goodness-of-fit)
— etc.
loc mn=(min x to plot)
loc mx=(max x to plot)
loc pts=`mx'-`mn'+1
tw fun y=f(x), ra(`mn' `mx') n(`pts') recast(scatter) [any options applicable to –twoway scatter– can go here]