列出解决r中多项式的所有解

时间:2016-07-24 22:34:18

标签: r

我有一个50x + 50y = 200的等式,我在R中编程为:

a = array(c(50,50),c(1,2))
b = c(200)
qr.solve(a,b)

但事实证明只显示4和0。所以,我的问题是有没有办法列出所有解决方案?

提前致谢

1 个答案:

答案 0 :(得分:3)

如果你需要非负整数解,你可以通过强力来轻松地做到这一点,因为0< = x< = 4和0< = y< = 4:

dd <- expand.grid(x=0:4,y=0:4)
with(dd,dd[50*x+50*y==200,])
##    x y
## 5  4 0
## 9  3 1
## 13 2 2
## 17 1 3
## 21 0 4

(忽略第一列,它只是原始数据框中的行号)