syms x y z;f=(1/2/pi)*exp(-(x^2+y^2)/2)*(1+x*y*exp(-(x^2+y^2-2)/2))-z;z=solve(f,'z'); for p=1:length(z)ezmesh(z(p));hold onend
