###   Obiettivo 1   ###     
  Valore dell'estremo superiore di oscillazione? (z+ > 0)
    2

  Integrazione con il metodo del punto medio:
    numero dei sotto-intervalli = 10^3  ;  errore = 3.465078e-01
    numero dei sotto-intervalli = 10^4  ;  errore = 1.095113e-01
    numero dei sotto-intervalli = 10^5  ;  errore = 3.462849e-02
    numero dei sotto-intervalli = 10^6  ;  errore = 1.095043e-02

  Stampa dell'ultimo valore calcolato:
    periodo = 13.944391



     ###   Obiettivo 2   ###     

  Integrazione con il metodo del punto medio:
    numero dei sotto-intervalli = 10^3  ;  errore = 3.487083e-04
    numero dei sotto-intervalli = 10^4  ;  errore = 1.101996e-05
    numero dei sotto-intervalli = 10^5  ;  errore = 3.484593e-07
    numero dei sotto-intervalli = 10^6  ;  errore = 1.102070e-08
    periodo = 13.949456



     ###   Obiettivo 3   ###     
  Valore del tempo finale Tf della soluz. num. dell'eq. diff.? (Tf > 0)
    10
  Studio degli errori indotti dall'applicazione del metodo
  di Heun (per l'integrazione numerica delle equazioni diff.
  di Newton) in funzione del passo di integrazione h.
    Tf / h     Delta E
  ----------------------
    1.0e+02    8.2e-03
    1.0e+03    1.7e-05
    1.0e+04    8.9e-08
    1.0e+05    8.2e-10
    1.0e+06    8.1e-12



     ###   Obiettivo 4   ###     
  Valore del rapporto tra T e h? (T/h deve essere uguale a
  una potenza intera di 10 tale che 100 <= T/h <= 1milione)
      10000
  Differenze tra posizione e velocita' al tempo T e quelle
  stesse quantita' corrispondenti alle condizioni iniziali:
      Delta z  =  2.1e-07
      Delta v  =  5.3e-07



     ###   Obiettivo 5   ###     
  Si veda il file 'periodo_sitn_circ.gnp' (scritto indipendentemente dal
  programma che produce questi messaggi); tale file contiene i comandi che
  consentono di visualizzare il grafico di un'orbita del problema circolare di
  Sitnikov utilizzando opportunamente il programma gnuplot.