Found out that wrong pressure has been used.

The right hand side of the differential equation returns the
gravitational pressure. How pressure is calculated therefore has to be
reimagined.