The divergence is correct

Not sure what you have in mind for rotF

but it is 5z^2 i +3xz^2 j +4 y^4x k

The middle row of the determinant are the operators d/dx d/dy and d/dz

and the last row are the x , y , and z components of F not the parial derivatives