#-------------------------------------------------------------------------
# Dati delle orbite dei pianeti del sistema solare
# https://nssdc.gsfc.nasa.gov/planetary/factsheet/planet_table_ratio.html
# semiasse maggiore in unità di quello terrestre
# 
#   Variante di pianeti_log.R, con 
# - plot 'normale'   
# - plot dei i VALORI DEI LOGARITMI decimali
#    
#  GdA marzo 2025
#--------------------------------------------------------------------------

pausa <- function() { cat ("\n >> press Enter to continue\n"); scan() }

# periodo di rivoluzione in unità dell'anno terrestre 
T <- c(0.241, 0.615, 1, 1.88, 11.9, 29.4, 83.7, 163.7)
# 'distanza' dal Sole (-> semiasse maggiore) in funzione dell'U.A.
d <- c(0.387, 0.723,  1 , 1.52, 5.20, 9.57, 19.17, 30.18)
n = length(T)

# limitiamo il nr di pianeti da analizzare (min 2, max 8)
n1 = as.integer(readline("Numero di pianeti (2-8; default 8): "))
if ( is.na(n1) ) n1 = 8
if (n1 > 1 && n1 <= 8) {
    n <- n1
    T <- T[1:n]
    d <- d[1:n]
    cat(sprintf(" => Consideriamo soli i primi %d pianeti \n", n))
} else {
    cat(sprintf(" %d non è valido \n", n1))
    stop()
}

# primo plot, d vs T
plot(T, d, pch=19, xlab='T (y)', ylab='d (R_T)',
     col='blue', main="Pianeti del sistema solare: 'distanza' vs periodo")
grid()
pausa()

# linearizzazione, plottando i logaritmi decimali delle grandezze 
# ATT! sono 'adimensionali' in quanto sono state divise per
# unità opportune (vedi sopra)
plot(log10(T),log10(d), pch=19, xlab='log10( T / y )', ylab='log10( d / R_T )', col='blue',
     main='Pianeti del sistema solare -- linearizzazione mediante i logaritmi decimali')
grid()