divergence_sym(V,X,​coordinate_system)

This Function calculates the divergence of the 3D symbolic vector in Cartesian, Cylindrical, and Spherical coordinate system.
function Div = divergence_sym(V,X,coordinate_system)
V is the 3D symbolic vector field
X is the parameter which the divergence will calculate with respect to.
coordinate_system is the kind of coordinate system at which the vector field is specified.
the divergence is calculated according to (Engineering Electromagnetics Sixth Edition William H. Hayt, Jr. . John A. Buck)
Example (1):
G = x*z*exp(2*y)*[z x*z x]
Div=divergence_sym(G, [x,y,z],'Cartesian')
Example (2):
G=0.2*r^3*phi*(sin(theta))^2*[1 1 1]
divergence_sym(G, [r theta phi],'Spherical')
Example (3):
G=2*rho^2*[cos(5*phi) -sin(5*phi) 1]
DivG=divergence_sym(G, [rho, phi, z],'cylindrical')