T = (z-y)/r + (x^2 + y^2)^0.5/s ; dT/dy =