r<-c(-10:10)
r<-r/10
sigma<-matrix(0,2,2)
sigma[1,1]<-1
sigma[2,2]<-1
x<-matrix(0,100,2)
dados<-matrix(0,100,42)
for (i in 1:21){
sigma[1,2]<-r[i]
sigma[2,1]<-r[i]
B<-mroot(sigma)
x[,1]<-rnorm(100)
x[,2]<-rnorm(100)
j<-1+(i-1)*2
dados[,j:(j+1)]<-x%*%t(B)
}
win.graph()
par(mfrow=c(3,2))
plot(dados[,1],-dados[,2],xlab="x",ylab="y",main="r=-1")
plot(dados[,3],dados[,4],xlab="x",ylab="y",main="r=-0.9")
plot(dados[,5],dados[,6],xlab="x",ylab="y",main="r=-0.8")
plot(dados[,7],dados[,8],xlab="x",ylab="y",main="r=-0.7")
plot(dados[,9],dados[,10],xlab="x",ylab="y",main="r=-0.6")
plot(dados[,11],dados[,12],xlab="x",ylab="y",main="r=-0.5")
win.graph()
par(mfrow=c(3,2))
plot(dados[,13],dados[,14],xlab="x",ylab="y",main="r=-0.4")
plot(dados[,15],dados[,16],xlab="x",ylab="y",main="r=-0.3")
plot(dados[,17],dados[,18],xlab="x",ylab="y",main="r=-0.2")
plot(dados[,19],dados[,20],xlab="x",ylab="y",main="r=-0.1")
plot(dados[,21],dados[,22],xlab="x",ylab="y",main="r=0")
plot(dados[,23],dados[,24],xlab="x",ylab="y",main="r=+0.1")
win.graph()
par(mfrow=c(3,3))
plot(dados[,25],dados[,26],xlab="x",ylab="y",main="r=+0.2")
plot(dados[,27],dados[,28],xlab="x",ylab="y",main="r=+0.3")
plot(dados[,29],dados[,30],xlab="x",ylab="y",main="r=+0.4")
plot(dados[,31],dados[,32],xlab="x",ylab="y",main="r=+0.5")
plot(dados[,33],dados[,34],xlab="x",ylab="y",main="r=+0.6")
plot(dados[,35],dados[,36],xlab="x",ylab="y",main="r=+0.7")
plot(dados[,37],dados[,38],xlab="x",ylab="y",main="r=+0.8")
plot(dados[,39],dados[,40],xlab="x",ylab="y",main="r=+0.9")
plot(dados[,41],dados[,42],xlab="x",ylab="y",main="r=+1")