(I(x,y,z)=(Subscript[I, 0])exp(-(((x^2)+(y^2))/(2(Subscript[σ, x]^2)))-((z^2)/(2(Subscript[σ, z]^2))))(3))