RD(x,y,z)=RD(0,y,z)-3xhz(1+O(xh)),