how to compute curl for velocity components expressed interms of spherical coordinate​s(r,alpha,​beta).

How mupad interface will do it and how can i plot its resulting curl in matlab. Thanks