Non mi trovo con il calcolo del potenziale. Andiamo ad integrare F3:
$\int [x e^{-y}-ye^{-z}+xy(1-z)e^{-x-z} ]dz $
$xze^{-y}+ye^{-z}+\int xye^{-x-z}dz-xy\int ze^{-x-z}dz = $
$ xze^{-y}+ye^{-z}-xye^{-x-z}-xy\int ze^{-x-z}dz$
Quest'ultimo lo faccio per parti:
$\int ze^{-x-z}dz = -ze^{-x-z}-\int -e^{-x-z} dz =-ze^{-x-z} -e^{-x-z} =-(z+1) e^{-x-z}$
Dunque otteniamo:
$U(x,y,z) = xze^{-y}+ye^{-z}-xye^{-x-z}+xy(z+1)e^{-x-z}$
che possiamo semplificare come:
$U(x,y,z) = xze^{-y}+ye^{-z}+xyz e^{-x-z} + h(x,y)$
Derivando rispetto a y e uguagliando a $F_2$ otteniamo:
$ -xze^{-y} + e^{-z} + xze^{-x-z} +h_y(x,y) = e^{-z}-xze^{-y}+xze^{-x-z}$
da cui semplificando:
$h_y(x,y) = 0$
e dunque h è costante rispetto a y cioè dipende solo da x: $h=h(x)$
Ora derivando invece rispetto a x e uguagliando a $F_1$:
$ze^{-y} + yze^{-x-z}-xyze^{-x-z}+h_x = ze^{-y} + yz(1-x)e^{-x-z}$
otteniamo che:
$ h_x = 0$
e dunque h è costante anche rispetto a x, cioé $h=c$. Dunque:
$U(x,y,z) = xze^{-y}+ye^{-z}+xyz e^{-x-z} + c$
Noemi