Compute the electric field of a dipole from its potential and verify that it is a vacuum solution by computing the divergence:
Compute the equivalent fields in Cartesian coordinates for
p1 and visualize the equipotential surfaces and lines of force: