function x=normar1(xmean,xrho,xsd,sampsiz) ;
cm=xmean*(1-xrho) ;
xb=zeros(2*sampsiz,1) ;
z=randn(1) ;
xb(1,1)=xmean+xsd*z/sqrt(1-xrho*xrho) ;
for i=2:(2*sampsiz) ;
 xb(i,1)=cm+xrho*xb(i-1,1)+xsd*rand ;
end ;
x=xb(sampsiz+1:2*sampsiz,1) ;
