#------------------------------------------------------
# Potenziale dentro e fuori la Terra
# in unità del modulo del potenziale sulla superficie 
# e con la distanza in unità di raggi terrestri
#
# GdA 12/11/2021
#-------------------------------------------------------


xM = 10  # x è il rapporto r/RT

curve(-1/x, 1, xM, xlim=c(0, xM), ylim=c(-3/2, 0), lwd=2, col='blue',
      xlab='r / RT', ylab=' VG / |VG(RT)| ')
curve(-1/2*(3-x^2), 0, 1, lty=2, lwd=2, col='blue', add=TRUE)
abline(h=0, col='gray')
abline(h=-1/2, col='gray')
abline(h=-1, lty=2, col='gray')
for(i in 0:xM) abline(v=i, lty=2, col='gray')