转载：R做空间三维图问题

在R中可使用rgl这个包来绘制三维图像，比R中自带的函数要好。

http://cos.name/cn/topic/12269

亮点：

广告语：rgl是您居家旅行画三维图形必备之良药。

fn = function(x, y) {
x^2 + y^2 - 10 * (cos(2 * pi * x) + cos(2 * pi * y)) + 20
}
a = 5
x = seq(-a, a, 0.05)
y = seq(-a, a, 0.05)
z = outer(x, y, fn)
persp(x, y, z)
library(rgl)
zorder = rank(z)
persp3d(x, y, z, col = rainbow(as.integer(max(zorder)))[zorder])