R编程：线性方程

Question

如何在R代码中说明给定方程式的某些可能值的列表？例如（这只是一个随机方程式，可随意使用任何合适的公式）：

  

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

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

Answer 1

您可以使用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