Hi all,
I'm wanting to write a small program simulating a 1D lattice with some motion. I have the equation:
m_{n}\frac{d^{2}u_{n}}{dt^{2}}= k_{n,n+1}(u_{n+1}-u_{n})+k_{n-1,n}(u_{n-1}-u_{n})
Then using a simple trial plane wave (u_{n}=Ae^{-i\omega t}). It boils down to:
- \omega^{2}...