Edição das 16h45min de 4 de maio de 2022 (ver código-fonte)
a_{m} = \vec{x}\cdot{\vec{n}}
Isso se segue imediatamente da equação (13). Com essas amplitudes e as frequências de cada modo dadas por (15), é simples resolver a parte da energia potencial (<math>\frac{1}{2}\omega_{m}^2a_{m}^2</math>) da equação (14). A parte da energia cinética pode ser resolvida se guardarmos as amplitudes calculadas no tempo anterior e fizermos a seguinte discretização:
\frac{da_{m}}{dt} \approx \frac{a_{i}-a_{i-1}}{\Delta t}
Edição das 16h53min de 4 de maio de 2022
O trabalho anterior contém apenas um pequeno erro no cálculo da parte potencial da energia. O código está usando a soma da amplitude no tempo anterior com a amplitude no tempo atual, o que está evidentemente incorreto. O cálculo deve usar apenas a amplitude no tempo atual.
Disponível em "http://fiscomp.if.ufrgs.br/index.php?title=Oscilações_Acopladas/Problema_de_Fermi-Pasta-Ulam-Tsingou&oldid=7630"