xy(y')^2+(x^2-y^2-u^2)y' -xy
= abcosθsinθ(-b/a cotθ)^2 + (a^2cos^2θ-b^2sin^2θ-u^2)(-b/a cotθ)-absinθcosθ
= abcosθsinθ(-b/a cotθ)^2 + (a^2sin^2θ-b^2cos^2θ)(-b/a cotθ)-absinθcosθ
= b/a cotθ( b^2 cos^2θ + a^2sin^2θ - b^2cos^2θ - a^2 sin^2θ)
=0
(3)も一緒