SET FORMAT=F8.3.
INPUT PROGRAM.
compute #rho=.50.
compute #n=25.
compute #mux=50.
compute #sigmax=10.
compute #muy=50.
compute #sigmay=10.
LOOP #I=1 TO #n.
compute x=RV.NORMAL(#mux, #sigmax).
compute #y=#rho*((x-#mux)/#sigmax)+RV.NORMAL(0, sqrt(1-#rho**2)).
compute y=#sigmay*#y+#muy.
END CASE.
END LOOP.
END FILE.
END INPUT PROGRAM.
GRAPH
  /SCATTERPLOT(BIVAR)=x WITH y
  /MISSING=LISTWISE .
CORRELATIONS
  /VARIABLES=x y
  /PRINT=TWOTAIL NOSIG
  /STATISTICS DESCRIPTIVES
  /MISSING=PAIRWISE .