another simple approach is to fit an intercept model in PROC MODEL to test bivariate (or multivariate) normality between X and Y or.....more variables.

proc model;
y=a;
x=b;
fit y x /normal;