R编程:线性方程

时间:2014-08-15 05:44:38

标签: r

如何在R代码中说明给定方程式的某些可能值的列表?例如(这只是一个随机方程式,可随意使用任何合适的公式):

  

对于正整数a,b和c,其公式为x ^ 3 + y ^ 2 = z。

如何测试x和y小于或等于1000的所有可能组合,c如何满足公式并检查变量是否也是有效输入?

1 个答案:

答案 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