Laboratorio di Calcolo (G. D'Agostini)

Soluzione numerica del moto di un oggetto sospeso a una molla

[Nota: nella terza equazione 'k/m' è preso unitario! (*) Un risultato è riportato da Feynman in Table 9-1]

Soluzione numerica del moto di un oggetto sospeso a una molla

[Posizione, velocità e accelerazione in funzione del tempo
(normalizzate ai valori massimi, riportati nella didascalia)]

• studiare bene lo script cercando di capire il significato dei vari step durante le iterazioni (si raccomanda la lettura dei paragrafi 9-5 e 9-6 del primo libro de La Fisica di Feynman);
• provare a cambiare k e m per vedere come cambiano periodo e velocità massima (potrebbe essere necessario riaggiustare 'ylim');
• provare anche a cambiare la granularità temporale, controllata dalla variabile dt;
• cominciare a pensare a un programma in C (e possibilmente scriverlo) che riproduca questi risultati (da controllare mediante stampe e con pochi punti, ovvero tMax pari a una decina di dt), seppur senza la parte grafica.
  [Nota: a lezione non abbiamo ancora affrontato i vettori nel linguaggio C, quindi si può usare per ora il trucco di memorizzare i valori vecchi in apposite variabili, come xOld e yOld nel lancio della moneta. Ma, ovviamente, se si dà un'occhiata ai tutorial raccomandati non sarà difficile usarli.]