Q <- function(x,y) sqrt(x^2+y^2); x1 <- y1 <- -4; x2 <- y2 <- 4 x <- seq(x1,x2,len=21); y <- seq(y1,y2,len=21); z <- outer(x,y,Q) # Traccio il grafico e lo metto in fig # Il comando col=NULL serve a fare il grafico senza colore, trasparente fig <- persp(x,y,z,theta=35,phi=15,scale=FALSE,col=NULL,ticktype="detailed") # Il contorno del piano con cui seziono il cono lines(trans3d(c(-4,4,4,-4,-4),c(-4,-4,4,4,-4),4,fig), col="red",lwd=2) # Il cerchio base superiore phi <- seq(0,2*pi,len=201);xr <- cos(phi)*4;yr <- sin(phi)*4 lines(trans3d(xr,yr,4,fig), col="blue", lwd=3) # Gli assi lines(trans3d(c(-4,4),0,0,fig), col="red") lines(trans3d(0,c(-4,4),0,fig), col="red") lines(trans3d(0,0,c(0,4),fig), col="red")