更新以匹配已编辑问题的预期输出:
g=c(0.65,0.70,0.75,0.80,0.85)
VijkArray <- array(0,dim=c(3,2,3))
result.1 <- array(0, dim = c(3,2,3*length(g)))
for(t in 1:length(g)){
for (i in 1:3) {
for (j in 1:3) {
for (q in 1:2) {
x <-xnew1[i,q,j]
VijkArray[i,q,j]<-Vijk(x,q)
}
}
}
indices <- ((t-1)*3+1):(t*3)
result.1[,,indices] <- VijkArray
}
输出:
print(result.1)
>, , 1
>
> [,1] [,2]
>[1,] 0.000000 0.5285917
>[2,] -1.433881 -0.5823365
>[3,] -1.866261 -1.1893313
>
>, , 2
>
> [,1] [,2]
>[1,] -0.5823365 0.4061262
>[2,] -1.6576946 -0.9137839
>[3,] -2.0629455 -1.4338807
>
>, , 3
>
> [,1] [,2]
>[1,] -0.9137839 0.2588162
>[2,] -1.8662608 -1.1893313
>[3,] -2.2500000 -1.6576946
>
>, , 4
>
> [,1] [,2]
>[1,] 0.000000 0.5032942
>[2,] -1.433881 -0.5823365
>[3,] -1.866261 -1.1893313
>
>, , 5
>
> [,1] [,2]
>[1,] -0.5823365 0.3789291
>[2,] -1.6576946 -0.9137839
>[3,] -2.0629455 -1.4338807
>
>, , 6
>
> [,1] [,2]
>[1,] -0.9137839 0.2332582
>[2,] -1.8662608 -1.1893313
>[3,] -2.2500000 -1.6576946
>
>, , 7
>
> [,1] [,2]
>[1,] 0.000000 0.4792073
>[2,] -1.433881 -0.5823365
>[3,] -1.866261 -1.1893313
>
>, , 8
>
> [,1] [,2]
>[1,] -0.5823365 0.3535534
>[2,] -1.6576946 -0.9137839
>[3,] -2.0629455 -1.4338807
>
>, , 9
>
> [,1] [,2]
>[1,] -0.9137839 0.2102241
>[2,] -1.8662608 -1.1893313
>[3,] -2.2500000 -1.6576946
>
>, , 10
>
> [,1] [,2]
>[1,] 0.000000 0.4562733
>[2,] -1.433881 -0.5823365
>[3,] -1.866261 -1.1893313
>
>, , 11
>
> [,1] [,2]
>[1,] -0.5823365 0.3298770
>[2,] -1.6576946 -0.9137839
>[3,] -2.0629455 -1.4338807
>
>, , 12
>
> [,1] [,2]
>[1,] -0.9137839 0.1894646
>[2,] -1.8662608 -1.1893313
>[3,] -2.2500000 -1.6576946
>
>, , 13
>
> [,1] [,2]
>[1,] 0.000000 0.4344368
>[2,] -1.433881 -0.5823365
>[3,] -1.866261 -1.1893313
>
>, , 14
>
> [,1] [,2]
>[1,] -0.5823365 0.3077861
>[2,] -1.6576946 -0.9137839
>[3,] -2.0629455 -1.4338807
>
>, , 15
>
> [,1] [,2]
>[1,] -0.9137839 0.170755
>[2,] -1.8662608 -1.189331
>[3,] -2.2500000 -1.657695