#--------------------------------------------------------
# Lente convergente con oggetto in diversi punti
# (o anche specchio concavo)
#
# GdA gennaio 2022
#--------------------------------------------------------

f = 1  # lunghezza focale unitaria 
       # -> p e q saranno in unità della lunghezza focale

par(mfrow=c(2,1)) # 2 plot insieme
old.mar = par("mar")
par(mar=c(4.0,4.2,0.4,0.8))   # margini ottimizzati
curve(p*f/(p-f), 0.01*f, 0.95*f, xname='p', xlim=c(0,3), ylim=c(-10,10),
      xlab='p / f', ylab='q / f', col='blue', lwd=2, xaxs="i")
curve(p*f/(p-f), 1.1*f, 4*f, xname='p', col='blue', lwd=2, add=TRUE)
abline(h=0, lty=2, col='gray')
abline(v=1, lty=2, col='gray')
abline(h=1, lty=2, col='blue')
abline(h=-1, lty=2, col='blue')
par(mar=c(4.3,4.2,1.0,0.8))   # margini ottimizzati
curve(-f/(p-f), 0.01*f, 0.95*f, xname='p', xlim=c(0,3), ylim=c(-10,10),
      xlab='p / f', ylab='M', col='red', lwd=2, xaxs="i")
curve(-f/(p-f), 1.01*f, 4*f, xname='p', col='red', lwd=2, add=TRUE)
abline(h=0, lty=2, col='gray')
abline(v=1, lty=2, col='gray')
abline(h=1, lty=2, col='red')
abline(h=-1, lty=2, col='red')
par(mar=old.mar)
par(mfrow=c(1,1)) # resettiamo il parametro mfrow