A explicação do script a seguir está no artigo R na Prática – Parte 31

#fonte rdocumentation.org
z <- 2 * volcano #exagero vertical   
x <- 10 * (1:nrow(z))   
y <- 10 * (1:ncol(z))
def.par=par(no.readonly=TRUE)
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_22.jpeg",width=5,height=5,
#units="in",res=300)
par(bg = "white")
persp(x, y, z, theta = 135, phi = 30, 
col = "green3", scale = FALSE, 
ltheta = -120, shade = 0.75, 
border = NA, box = FALSE)
par(def.par)
#dev.off()