y := 100*exp(-x^2/2)/sqrt(2*Pi); plot(y,x=-4..4,labels=["Standard Units",""]); Example: Answer: Aaron_SUs := (48 - 53)/18.; Cecilia_SUs := (74-53)/18.; Question: Answer: Score := 53*pts + (-2.5)*18*pts; read `scores.mpl`; with(statplots): histogram(scores,area=count); read `ToSus.mpl`; with(stats): with(describe): AVE := evalf(mean(scores),3); SD := evalf(sqrt(variance(scores)),3); scrs_in_SUs := evalf(ToSUs(scores,AVE,SD),2); with(statplots): with(plots): Hist := histogram(scrs_in_SUs): NCurve := plot(stats[statevalf,pdf,normald], -3..3, color=red): display({Hist,NCurve}); sort(scores);