Derivation of TOV equations, now with units (to add: numerical solution, M(R))