For loop R中嵌套for循环的问题

我试图通过复制他们在课堂上给我们的定量系统性能分析算法来练习R。输出为系统中一定数量的项目（n）和一定数量的机器（k）的队列长度（Q）、吞吐量（X）和服务时间（R）。 我从一个简化的版本开始，当机器的数量=1并且代码正在运行时

N1 <-c(1,2,3)
K1 <- 1
Q <- 0 
R <- 0
D <- 3 # service rate of the machine
for(z in 1:length(N1))
{if(z==1){R[z] <-D} else{R[z] <- 3*(1+Q[z-1])}
X<- z/R[z];
Q[z] <- X*R[z]}

---

N1您必须分配二维矩阵R

使用：

R
N1 <-c(1,2)
K <- c(1,2,3,4)
D <- c(3,4,3,6) 
Q <- 0 
R <- 0
X <-0 
for(z in 1:length(N1)) 
  {for(k in 1:length(K))
  {if(z==1){R[k,z] <-D[k]} else{R[k,z] <- D[k]*(1+Q[k,z-1])}
  X[z]<- z/sum(R[z]); 
  if(z==1){Q[k,z] <- X[z]*R[k,z]} else {Q[k,z] <- X[z]*R[k,z]}
 }}

R <- matrix(nrow=length(K), ncol=length(N1))

R <- 0