Consider scalar potential of electric field \phi (x,y)=\alpha xy\alpha is a constantx,y in limits -5,5, and assume constant \alpha =3. Add the following line
after surf =dots command to plot equipotential lines:
ax.contour(X, Y, Z, zdir='z', offset=-3, cmap='coolwarm').
C