Il grafico di: 2 2 x +y - ————— 1 2 f(x,y) = ——— e 2π ottenuto con i comandi seguenti:

f <- function(x,y) { q <- x^2+y^2; exp(-q/2)/(2*pi) }
x1 <- -5; x2 <- 5; y1 <- -5; y2 <- 5; z12 <- c(-0.02,0.18)
x <- seq(x1,x2,length=40); y <- seq(y1,y2,length=40); z <- outer(x, y, f)
persp(x,y,z,zlim=z12,theta=30,phi=20,expand=0.5,col="lightblue",ticktype="detailed")
# expand indica il fattore di espansione applicato all'altezza del box; con 1
# si ha una forma più "cubica". Se aggiungo d=... con un numero > 1 vedo
# da un punto di vista più lontano (le facce tendono ad diventare parallele),
# con uno < 1 lo avvicino.
persp(x,y,z,zlim=z12,theta=30,phi=20,expand=1,col="lightblue",ticktype="detailed",d=10)
persp(x,y,z,zlim=z12,theta=30,phi=20,expand=0.3,col="lightblue",ticktype="detailed",d=0.1)