We can think of the two surfaces as "equi-potential" surfaces of [tex]f(x,y,z)= z- x^2- y^2[/tex] and [tex]g(x,y,z)= x+ y+ z- e^{xyz}[/tex]. Their gradients, [tex]\nabla f= -2x\vec{i}- 2y\vec{j}+ \vec{k}[/tex] and [tex]\nabla g= (1- yze^{xyz})\vec{i}+ (1- xze^{xyz})\vec{j}+ (1- xye^{xyz})\vec{k}[/tex] are perpendicular to the surfaces and perpendicular to the tangent planes at each point.
In particular, at (0, 0, 1), [tex]\nabla f(0, 0, 1)= \vec{k}][/tex] so the tangent plane is [tex]0(x- 0)+ 0(y- 0)+ 1(z- 1)= 0[/tex] or z= 1. And [tex]\nabla g(0, 0, 1)= \vec{i}+ \vec{j}+ \vec{k}[/tex] so the tangent plane is [tex]1(x- 0)+ 1(y- 0)+ 1(z- 1)= x+ y+ z- 1= 0[/tex] or [tex]x+ y+ z= 1[/tex]. The angle between those two planes is the same as the ange between the two perendicular vectors and you can use [tex]\vec{u}\cdot\vec{v}= |\vec{u}||\vec{v}|cos(\theta)[/tex] to find that angle.
The length of [tex]\vec{k}[/tex] is 1, the length of [tex]\vec{i}+ \vec{j}+ \vec{k}[/tex] is [tex]\sqrt{3}[/tex], and their dot product is 1. So [tex]1(\sqrt{3})cos(\theta)= 1[/tex].