Example 5 (p. 734)
> F := [ x^2*y*z, 3*x*y*z^3, x^2-z^2 ];
> diverge( F, [x,y,z] );
> curl( F, [x,y,z] );
>