OldClaudio” post=79588In ogni caso la cllonna è la sesta, da quel che dici; è ho il sospetto che quei numeri derivino da una differenza fra l’energia attesa e quella simulata dal calcolo. La differenza da 1 è minima ma si tratta sostanzialmente di runore dovuto agli errori di scalamento e troncamento dell’algoritmo che ti porta a quel risultato. Hai studato un po’ di clacolo numerico? non occore andare molto a fondo, ma l’errore che si ottiene da differenza fra numeri relativamante grandi ma molto vicini, è fortemente “rumoroso” e quindi non è assolutamente attendibile.
Se invece quella sesta colonna è il rapporto fra l’energia attesa e quella simulata, la retta orizzontale all’ordinata 1 va più che bene, benissimo! e le oscillazioni che risultano invisibili nel grafico devono restare invisibili perché sono rumore.
Hai ragione, la colonna è la sesta. Avevo in mente il numero di colonna in pgfplots, che parte a contare da zero. =).
Comunque non è nessuna delle due alternative di cui hai parlato. È semplicemente l’energia meccanica (nessuna differenza, nessun rapporto). Fisicamente dovrebbe essere costante ma numericamente voglio (e devo) mostrare che oscilla (anche se poi negli schemi di Euler-Cromer e Stormer-Verlet è conservata come media temporale). Questo per gli errori introdotti dall’algoritmo: ma sono proprio questi errori che voglio studiare.
Spero che ora sia più chiaro il mio scopo.
Grazie,
R.