Write a MATLAB script using the quiver and contour commands to visualize the field and its divergence. Assume the region of interest to be -20-20\leqslant x\leqslant20 ,-40\leqslant y\leqslant 40.

Hint: Mesh the 2D space as

gx=-20:1:20;

gy=-40:1:40