# 2004 年 11 月 11 日作成
rep.no<-100000
ell.1<-rep(0,rep.no)
ell.2<-rep(0,rep.no)
ell.1d<-ell.1;ell.2d<-ell.2
#ell.1; ell.2
for (i in 1:rep.no){
xw.1<-rnorm(10,0,1)
xw.2<-rnorm(10,0,1)
xw<-matrix(c(xw.1,xw.2),10,2)
#xw
w<-t(xw)%*%xw
#w
eigen.w<-eigen(w)
#eigen.w
lam<-c(eigen.w$values)
hmat<-matrix(eigen.w$vectors,2,2)
#lam; hmat
ell.1d[i]<-lam[1]/11; ell.2d[i]<-lam[2]/9
ell.1[i]<-lam[1]/10; ell.2[i]<-lam[2]/10
}
plot(density(ell.2d),xlim=c(0,4),ylim=c(0,2.0),xlab="",ylab="",main="density of eigenvalues")
par(new=T)
plot(density(ell.2),xlim=c(0,4),,xlab="",ylab="",ylim=c(0,2.0),main="")
mean(ell.1d);mean(ell.2d)
var(ell.1d); mean(ell.2d)
mean(ell.1);mean(ell.2)
var(ell.1); mean(ell.2)