R绘图隐式函数外部命令

我想绘制 x 和 y 的隐式函数：1 - 0.125 * y ^ 2 - x ^ 2 = 0.005我知道它可以绘制为等值线图，但使用“外部”命令时遇到问题 在下面的：

x<-seq(0.4,1.01,length=1000)
y<-seq(0,3,length=1000)
z<-outer(x,y,FUN="1-0.125*y^2-x^2=0.005")
contour(x,y,z,levels=0,drawpoints=FALSE)

我阅读了有关“外部”命令和矢量化函数的需要的常见问题解答（7.17），但我仍然陷入困境。

我认为您对“功能”的含义有点困惑。 函数中的所有运算 (+、-、​​^) 均已向量化，因此一切正常。

x <- seq(0.4,1.01,length=1000)
y <- seq(0,3,length=1000)
z <- outer(x,y,function(x,y) 1-0.125*y^2-x^2-0.005)
contour(x,y,z,levels=0,drawlabels=FALSE)

或者如果你想要一个小捷径：

library(emdbook)
curve3d(1-0.125*y^2-x^2-0.005,
        xlim=c(0.4,1.01),
        ylim=c(0,3),
        n=c(100,100),
        sys3d="contour",drawlabels=FALSE,levels=0)

这实际上比较慢，因为它使用了for内部循环而不是outer()，所以我将其设置为 100x100 而不是 1000x1000（无论如何，这对于这个示例来说是多余的），但它适用于无法轻松矢量化的更复杂的示例......