u = (1+xy)^z
∂u/∂x = yz(1+xy)^(z-1)
∂u/∂y = xz(1+xy)^(z-1)
∂u/∂z = [ln(1+xy)] (1+xy)^z