R编程：线性方程

Question

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

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

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

Answer 1

您可以使用expand.grid生成所有可能的值，然后子集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