###   Obiettivo 1   ###     
  E e J sono, rispettivamente, l'energia e il momento angolare
  corrispondenti a un fissato valore del parametro eps e alle
  condizioni iniziali ( x(0), y(0), vx(0), vy(0) ).
  In particolare, le condizioni iniziali sono tali che: x(0) = 1,
  y(0) = 0, vx(0) = 0 e il valore di vy(0) verra' introdotto da
  tastiera tra poco; invece, il valore del parametro eps viene
  automaticamente calcolato a partire dal rapporto Jbarra/J
  che verra' anch'esso introdotto da tastiera tra poco.
  L'utente di questo programma DEVE introdurre dei valori di
  vy(0) e Jbarra/J tali che i corrispondenti valori di J e di
  E siano rispettivamente DIVERSO da zero e MINORE di zero,
  altrimenti non esistono le orbite 'a rosetta' che vogliamo
  studiare.  Per ottenere che l'energia totale E sia NEGATIVA
  e' bene introdurre dei valori non troppo grandi sia per quanto
  riguarda la velocita' vy(0) che il rapporto Jbarra/J

  Valore iniziale della velocita' verticale vy(0)? (vy(0) > 0)
      4
  Valore del rapporto tra Jbarra e J? (Jbarra/J >= 1 , E < 0)
      1.8
  Il corrispondente valore dell'energia e':
      E = -13.558418
  e il corrispondente valore del momento angolare e':
      J = 4.000000



     ###   Obiettivo 2   ###     
  Il valore del periodo di oscillazione tra pericentro e apocentro e':
    T_rho = 1.756633



     ###   Obiettivo 3   ###     
  Numero nlobi delle escursioni complete (pericentro-apocentro
  e ritorno apocentro-pericentro) da calcolare? (nlobi >= 1)
      9
  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. Si
  ricordi che il tempo finale T e' dato da T = T_rho * nlobi
     T / h     Delta E    Delta J
  ---------------------------------
    1.0e+03    3.5e-03    1.1e-03
    1.0e+04    3.5e-06    1.1e-06
    1.0e+05    3.5e-09    1.1e-09
    1.0e+06    3.6e-12    1.2e-12
    1.0e+07    4.7e-13    1.3e-13
    1.0e+08    2.9e-13    1.4e-14



     ###   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)
      1000000
  Differenze tra posizione e velocita' al tempo T e quelle
  stesse quantita' corrispondenti alle condizioni iniziali:
      Delta x  =  1.1e-13
      Delta y  =  8.1e-09
      Delta vx =  2.5e-08
      Delta vy =  8.1e-13



     ###   Obiettivo 5   ###     
  Si veda il file 'heun_calogero.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 di Calogero,
  utilizzando opportunamente il programma gnuplot.