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