8月29日（金）午前

行列化

午前中全部使ってやっと成功した。終わってみるとなんていうことはない話だが、かなり苦労した。

やはりまだまだ習熟が足りない。

f1<-function(i){
y<-c()
for (j in 0:(9-i)){
x<-c(i,j,9-i-j)
y<-c(y,x)}
;y
}

f2<-function(){
z<-c()
for(i in 0:9)
z<-c(z,f1(i));z
}

matrix(f2(),ncol=3,byrow=TRUE)

> matrix(f2(),ncol=3,byrow=TRUE)
      [,1] [,2] [,3]
 [1,]    0    0    9
 [2,]    0    1    8
 [3,]    0    2    7
 [4,]    0    3    6
 [5,]    0    4    5
 [6,]    0    5    4
 [7,]    0    6    3
 [8,]    0    7    2
 [9,]    0    8    1
[10,]    0    9    0
[11,]    1    0    8
[12,]    1    1    7
[13,]    1    2    6
[14,]    1    3    5
[15,]    1    4    4
[16,]    1    5    3
[17,]    1    6    2
[18,]    1    7    1
[19,]    1    8    0
[20,]    2    0    7
[21,]    2    1    6
[22,]    2    2    5
[23,]    2    3    4
[24,]    2    4    3
[25,]    2    5    2
[26,]    2    6    1
…以下略

午後はこれを使って確率を計算してみる。また苦労しそうな予感がする。

以上