如何在R代码中说明给定方程式的某些可能值的列表?例如(这只是一个随机方程式,可随意使用任何合适的公式):
对于正整数a,b和c,其公式为x ^ 3 + y ^ 2 = z。
如何测试x和y小于或等于1000的所有可能组合,c如何满足公式并检查变量是否也是有效输入?
答案 0 :(得分:4)
您可以使用expand.grid
生成所有可能的值,然后将其分配到符合条件的值:
vals <- expand.grid(x=seq(1000), y=seq(1000))
subset(vals, x^3 + y^2 == 108)
# x y
# 8003 3 9
# 9002 2 10