In “Gravity anomalies of 2D bodies with variable density contrast,” Zhang et al. derive a set of equations (9 and 10 in their paper) for a gravity anomaly at the origin of the coordinate system for a polygon mass body when the density contrast is a polynomial function of x and z. For convenience of discussion, I reproduce and renumber their equations 9 and 10 as, respectively