Denotiamo con η(t) il numero di attualmente positivi nell'istante t. Per quanto precede, la variabile tempo è in realtà campionata con intervallo di campionamento Δ=1d , a partire da t0 = 0=24/02. Plottando i dati ηk e interpolando otteniamo il grafico di figura:
È interessante studiare l'evoluzione dinamica di questo sistema fisico quando è sottoposto a una sollecitazione esterna. Il codice non differisce molto: basta aggiungere la legge della forza (o meglio, di singola componente cartesiana) a secondo membro delle equazioni differenziali del moto. Tuttavia, il caso più interessante è quello in cui tale sollecitazione è una variabile aleatoria, ovvero un noise. (altro…)