Example 5 (p. 734)

> F := [ x^2*y*z, 3*x*y*z^3, x^2-z^2 ];

F := [x^2*y*z, 3*x*y*z^3, x^2-z^2]

> diverge( F, [x,y,z] );

2*x*y*z+3*x*z^3-2*z

> curl( F, [x,y,z] );

VECTOR([-9*x*y*z^2, x^2*y-2*x, 3*y*z^3-x^2*z])

>