derive z=ln((x+1)/(y+1))