phL <- function(I,L) { t <- -10; par( mai = c(0.2,0.2,0.5,0.3) )
x <- c(0,L); y <- c(0,L); z <- array(c(0,0,0,0),dim=c(2,2)); z0 <- L
pr<- persp(x,y,z,zlim=c(0,z0),xlim=x,ylim=y,xlab="x",ylab="y",zlab="z",box=FALSE,border="white",phi=I,theta=t)
x1<- seq(-20,30,len=200); lines(trans3d(x1, z=0, y=(x1-5)^2, pmat=pr),col="blue")
  lines(trans3d(c(0,0),c(0,0),c(0,z0),pmat=pr),col="red",lty=3)
  lines(trans3d(c(0,z0),c(0,0),c(0,0),pmat=pr),col="black",lty=3)
  lines(trans3d(c(0,0),c(0,z0),c(0,0),pmat=pr),col="green4",lty=3)
  title(paste("x-bl y-gr z-re  I, L =",I,",",L)) }
phL(1,10)
# Una figura vista da un'inclinazione di 1° dal piano x,y.
# Assi lunghi 10, colorati nero, verde e rosso.  Che figura è?
# Prova a cambiare l'inclinazione I (ed eventualmente la
# lunghezza degli assi L):  batti phL(I,L) con i nuovi valori