curl_sym(V,X,coordi​nate_system)

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