pyTOVpp Notebook by N. Stergioulas on Tolman-Oppenheimer-Volkoff solutions with piecewise polytropic EOS
pyTOVpp computes a nonrotating relativistic star obeying a piecewise polytropic equation of state
Data for a convergence study is included at the end, showing near 4th-order convergence to machine accuracy