#---- 
pausa <- function() {
  cat ("\n >> Guarda il plot e dai enter per continuare\n")
  scan()
}

G <- 6.67e-11

#--- Giove  --------------------------------------------------------- 
nome <- c('Io', 'Europa', 'Ganimede', 'Callisto')
d <- c(422, 671, 1070,  1883) * 1e6 # metri  ('Semiasse maggiore')
T <- c(1.769, 3.551, 7.155, 16.689) * 86400 # s

plot(T^2, d^3, main="Satelliti di Giove", xlab='T^3 (s^2)', ylab='d^3 (m^3)')
pausa()

plot(T, d, log='xy',  main="Satelliti di Giove (plot 'log-log')", xlab='T (s)', ylab='d (m)')
pausa()

GM <- (4*pi^2) * d^3 / T^2
M.Giove <- GM/G
print(M.Giove)
media.M.Giove <- mean(M.Giove)
print("valore medio:")
print(media.M.Giove)

#Valore Wiki:  Mass: 1.898E27 kg (317.8 Earth mass)

#--- Sole  -----------------------------------------------------------
UA <- 149.6e9     # m
anno <- 31557600  # s
nome <- c('Mercurio', 'Venere', 'Terra', 'Marte', 'Giove', 'Saturno')
d <- c(0.387, 0.723, 1,  1.524, 5.203, 9.537) * UA
# http://it.wikipedia.org/wiki/Lista_dei_pianeti_del_sistema_solare
T <- c(0.241, 0.615, 1, 1.881, 11.86, 29.45) * anno # s

plot(T^2, d^3, main="Pianeti del sistema solare", xlab='T^3 (s^2)', ylab='d^3 (m^3)')
pausa()

plot(T, d, log='xy',  main="Pianeti del sistema solare (plot 'log-log')", xlab='T (s)', ylab='d (m)')
pausa()

GM <- (4*pi^2) * d^3 / T^2
M.Sole <- GM/G

print(M.Sole)
media.M.Sole <- mean(M.Sole)
print("valore medio:")
print(media.M.Sole)



# Wiki: 1.9891E30 kg