-Mは-y+1ではなく、-y-1です。だから、2x-y+1は2x-Mになりません。
そこが計算が合わない理由です。

そこで、M = 2x+1とおけば(M-y)(M+y) = M^2-y^2とできます。