Runge Kutta implemented. Runs:

- no pressure
- peak in momentum rhov
- constant density
- plots derivative operator