phi (x,y,z) = exp (-(x/x0)