I'll denote vectors in boldface. So the given vector field is
<em>B</em> = <em>x</em>² <em>x</em> - <em>y</em>² <em>y</em> + 3<em>z</em> <em>z</em>
<em></em>
Compute the divergence:
∇ • <em>B</em> = (∂/∂<em>x</em> <em>x</em> + ∂/∂<em>y</em> <em>y</em> + ∂/∂<em>z</em> <em>z</em>) • (<em>x</em>² <em>x</em> - <em>y</em>² <em>y</em> + 3<em>z</em> <em>z</em>)
∇ • <em>B</em> = ∂/∂<em>x </em>[<em>x</em>²] + ∂/∂<em>y</em> [-<em>y</em>²] + ∂/∂<em>z</em> [3<em>z</em>]
∇ • <em>B</em> = 2<em>x</em> - 2<em>y</em> + 3
Compute the curl:
∇ × <em>B</em> = (∂/∂<em>x</em> <em>x</em> + ∂/∂<em>y</em> <em>y</em> + ∂/∂<em>z</em> <em>z</em>) × (<em>x</em>² <em>x</em> - <em>y</em>² <em>y</em> + 3<em>z</em> <em>z</em>)
∇ × <em>B</em> = ∂/∂<em>x</em> [<em>x</em>²] (<em>x</em> × <em>x</em>) + ∂/∂<em>x</em> [-<em>y</em>²] (<em>x</em> × <em>y</em>) + ∂/∂<em>x</em> [3<em>z</em>] (<em>x</em> × <em>z</em>)
… … … + ∂/∂<em>y</em> [<em>x</em>²] (<em>y</em> × <em>x</em>) + ∂/∂<em>y</em> [-<em>y</em>²] (<em>y</em> × <em>y</em>) + ∂/∂<em>y</em> [3<em>z</em>] (<em>y</em> × <em>z</em>)
… … … + ∂/∂<em>z</em> [<em>x</em>²] (<em>z</em> × <em>x</em>) + ∂/∂<em>z</em> [-<em>y</em>²] (<em>z</em> × <em>y</em>) + ∂/∂<em>z</em> [3<em>z</em>] (<em>z</em> × <em>z</em>)
∇ × <em>B</em> = 0
(since each partial derivative not along the main diagonal vanishes, and for any vector <em>a</em> we have <em>a</em> × <em>a</em> = 0)