generating normal distribution for confidence bands

[vogelh]

Using version7

Sorry to bother, but I'm not a full-time expert programmer.

Have a bivariate regression line plotted thru scatter points.
Still having problem coding upper and lower confidence band (say 95%) above regression line.

As in elementary textbooks, want the band to be closest to mean and then flare out or curve as
it goes further out from mean in either direction.

Have tried to use @dnorm(0.025) but haven't figured proper way to code this.
I can generate the regression line and outcome matrix without any problem
and place any part of the outcome into any matrix.

Is there a suggestion for coding this easily.
thanks.