接上次问题，我想找线性模型的若干组整数解，我自己编了个循环语句，但是有错误，希望高手们指点下：
Y<-matrix(0,200,4)
for(x1 in 0:200){
    for(x2 in 0:200){
         for(x3 in 0:200){
             for(x4 in 0:200)
       {if (x3+x4==98&x2+x4==98&x1+x2+x3+x4==200 )Y<-[,x1:x4]}
                          }
                     }
                 }

思想，我想挨个试验，如果是解的话，就存储到Y这个矩阵中。