from sympy.vector import * 
from sympy import 
symbols N=CoordSys3D('N') 
x,y,z=symbols('x y z')
A=N.x**2*N.y*N.z*N.i+N.y**2*N.z*N.x*N.j+
N.z**2*N.x*N.y*N.k 
delop=Del() 
curlA=delop.cross(A) 
display(curlA)
print(f"\n Curl of {A} is \n") 
display(curl(A))