Consider a uniform, infinitely thin disc of surface mass density $\sigma$ and radius $R$ placed in the $xy$-plane with its center as the origin. The gravitational potential at a point on the axis of the disc ($z$-axis) at a distance $z$ from the center of the disc is given by: $$\phi(z)= -2\pi G \sigma \left(\sqrt{R^2 +z^2}-z\right). $$
The $z$-component of the field at this point is given by $$f_z(z)=-2\pi G \sigma\left(\frac{z}{\sqrt{R^2+z^2}}-1\right).$$
Given the symmetry of the problem, shouldn't $\phi$ be symmetric functions of z? Here, we clearly have $\phi (-z) \neq \phi(z)$.
Also, the field should be antisymmetric, i.e $f_z(-z) = -f_z(z)$, which is not satisfied.