###   Obiettivo 1   ###     
  Valore dell'eccentricita' e? (0 <= e < 1)
    0.7
  Valore dell'anomalia media M? (0 <= M <= pigreco)
    1.2
  Calcolo del corrispondente valore dell'anomalia vera v,
  applicando il metodo di Newton all'equazione di Keplero.
  Il corrispondente valore dell'anomalia vera v e':
    v = 2.539778
  Test di verifica del calcolo precedente, utilizzando le
  relazioni inverse; l'errore assoluto | M(E) - M | vale:
    3.124671e-16






     ###   Obiettivo 2   ###     
  Valore dell'indice n? (n >= 1)
    4
  I valori della funzione di Bessel J_n e della sua derivata in
  corrispondenza a x=ne e calcolati con 2 integrali sono uguali a:
    1.066687e-01    1.203148e-01



     ###   Obiettivo 3   ###     
  Test : la differenza tra il valore di J_n(n*e) calcolato grazie a
  una serie di Taylor e quello calcolato con un integrale e' uguale a
    8.326673e-17
  Test : la differenza tra la derivata di J_n calcolata in x=n*e grazie
  a una serie di Taylor e quella calcolata con un integrale e' uguale a
    9.159340e-16



     ###   Obiettivo 4   ###     
  Integrazione numerica dell'equazione differenziale che e'
  risolta proprio dalla funzione di Bessel di prima specie J_n
  Nuovo valore dell'indice n? (n >= 1)
    2
  Valore iniziale dell'intervallo [xin,xfin]? (xin > 0)
    0.05
  Valore finale dell'intervallo [xin,xfin]? (xin < xfin <= 40)
    25
  Numero di passi da effettuare con metodo di Heun? (10 <= npassi <= 1000000)
    100000
  Test : confronto tra il valore di Jn(xfin) (calcolato tramite
  un'espansione in serie) e l'altro suo valore che si ottiene
  integrando un'opportuna eq. diff. risolta proprio da y(x)=Jn(x):
    | Jn(xfin) - y(xfin) | = 5.921262e-07



     ###   Obiettivo 5   ###     
  Scrittura su file dei dati necessari a effettuare un grafico della
  funzione di Bessel di prima specie J_n.
  Si veda il file 'bessel_mec_cel.gnp' (scritto indipendentemente dal
  programma che produce questi messaggi); tale file contiene i comandi
  che consentono di visualizzare il grafico della funzione di Bessel
  di prima specie Jn(x) utilizzando opportunamente il programma gnuplot.