(methods)=
# Methods
<!-- talk about why FD is used against FV and what are the basics of FD (it approximates derivatives at grid points: 1. a differential equation is discretized, 2. stores values at grid points, works nicely on regular grids where the domain is simple,  ) -->