D_y ln( (1-sqrt (x^2 + y^2 + z^3) / (1+sqrt (x^2 + y^2 + z^3) )