(i) div(F) = ∂x(x3z) + ∂y(-2yz2) + ∂z(xy3z4) = 3x2z - 2z2 + 4xy3z3
At (1, 1, 3), you just plug in x = 1, y = 1, z = 3 in the above.
div(F(1, 1, 3)) = 3•12•3 - 2•32 + 4•1•13•33 = 99
(ii) curl(F) = (∂y(xy3z4) - ∂z(-2yz2))i - (∂x(xy3z4) - ∂z(x3z))j + (∂x(-2yz2) - ∂y(x3z))k
= (3xy2z4 + 2z2)i - (y3z4 - x3)j + (0 - 0)k
At (1, 1, 3), you just plug in x = 1, y = 1, z = 3 in the above, you'll get what you want curlF(1, 1, 3).
(iii) It is basically asking for the divergence of the curl of F, which is zero (for any vector field with continuous partial derivatives) div(curl(F)) = 0.